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Abstract 

The non-relativistic energy levels of para-positronium are calculated in the 
quadrupole approximation of the interaction potential. This approximation tech- 
nique takes into account the anisotropy of the electrostatic electron-positron in- 
teraction in the lowest order. The states due to different values of the quantum 
number (l^) of angular momentum are found to be no longer degenerate as is the 
case in the conventional theory. The physical origin of this elimination of the con- 
ventional degeneracy may intuitively be attributed to the state-dependent inertial 
broadening of the rotating charge clouds; the corresponding anisotropic deformation 
(in the quadrupole approximation) lowers then the negative electrostatic interac- 
tion energy. The result of this influence of anisotropy is that the states with Iz = 
adopt smaller binding energy whereas the states with maximal value of | /j; | (for fixed 
principal quantum number n) have the largest binding energy within the angular 
momentum multiplet ( — |/z,max \ <lz < \lz. ,max|)- This yields a certain kind of elec- 
tric fine-structure splitting with the splitted RST levels being placed in a relatively 
narrow band around the (highly degenerated) conventional levels. 
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I Introduction and Survey of Results 

The present paper is intended to increase the accuracy of the RST predictions con- 
cerning the (non-relativistic) energy levels of para-positronium. In the lowest order of ap- 
proximation {spherically symmetric approximation), the corresponding RST predictions 
deviated from their conventional counterparts by some 5-10% for principal quantum num- 
bers up to n ~ 100 [1-4]. This provided us with sufficient motivation in order to consider 
now also the next higher order of approximation where the anisotropy of the interaction 
between positron and electron is taken into account, albeit only in the quadrupole approx- 
imation. But amazingly enough, this low degree of approximation is already sufficient in 
order to eliminate the energetic degeneracy which is present in the spherically symmetric 
approximation (and in the conventional theory). 

The latter kind of RST approximation is due to the assumption that the interaction 
potential between the positronium constituents (i. e. electron and positron) is spherically 
symmetric, i. e. it depends exclusively on the radial variable r of the spherical polar 
coordinates (r, 0) . Such a spherical symmetry of the interaction potential is adopted 
also in the conventional theory, namely in form of the Coulomb potential entering the 
conventional Hamiltonian H in the usual way: 

The corresponding energy eigenvalue problem can be solved exactly [5], and the conven- 
tional eigenvalues E^''' (n = 1, 2, 3, . . .) turn out as 

^<,.._^4._^^[eVl (1.2) 

qb = ^ ■ ■ ■ Bohr radius | . 

Me2 J 

Since the conventional quantum number of angular momentum £ ranges here from zero to 
n — 1, and its z-component £z from —i up to +i {—i < iz ^ +-^), the conventional degree 
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of degeneracy is thus given by 

n-l 

5^(2£+l)=n^ (1.3) 

This result is a mathematical consequence of the fact that the original three-dimensional 
eigenvalue problem can be separated by a product ansatz into two eigenvalue problems of 
their own, namely into the eigenvalue problem for the angular momentum operator and 
into a residual one-dimensional problem (in radial direction) for the energy. 

Such a separation of the original three-dimensional eigenvalue problem into a two- 
dimensional angular and a one-dimensional radial problem is also possible in RST, but 
only if one resorts to the spherically symmetric approximation for the interaction poten- 
tial. In this case, the angular momentum problem can also be solved exactly quite simi- 
larly as in the conventional theory (see Subsect.IV.2 ^^Angular Momentum Quantization 
in RST"), but the residual one- dimensional eigenvalue problem (in radial direction) for 
determining the energy eigenvalues is in RST much more complicated than in the conven- 
tional theory. And furthermore the conventional £^-degeneracy is eliminated in RST. The 
reason is here that the interaction potential (even in the spherically symmetric approxi- 
mation) depends on the considered quantum state and thus is not fixed once and forever 
as is the Coulomb potential in the conventional Hamiltonian H f ll.ip . Indeed, the RST 
interaction potential must be determined from a Poisson equation which has the charge 
density of the considered quantum state as its source, even in the spherically symmetric 
approximation (see ^^Spherically Symmetric Eigenvalue Equations^^ in Subsect.IV.3). 

The physical consequence of this state-dependence of the interaction potential is now 
that the states with different values of angular momentum are subjected to different 
interaction potentials and this entails of course the elimination of the otherwise present 
£^-degeneracy. Such a result seems not to be supported by the conventional theory [6]; 
but one must not forget the fact that the present elimination of the ^^-degeneracy in RST 
occurs in the quadrupole approximation; and therefore it may well be the case that the 
higher orders of approximation do weaken that effect of degeneracy elimination so that it 
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might escape the experimental verification by means of the present-day techniques. 

Indeed, there is some plausible reason why such an elimination should actually occur. 
Namely, the states due to different values of Iz niay intuitively be associated with differ- 
ent magnitudes of rotation (around the 2;-axis) which then would entail the emergence 
of different strengths of the centrifugal forces. As a consequene, the original spherically 
symmetric charge distributions would be flattened (or rolled out) in a different way by 
these different centrifugal forces, cf. Fig.s IV. A. 0— IV. A. 3. But such deformed electric 
charge distributions should then generate (via the Poisson equations) an electric interac- 
tion potential with a certain deviation from the purely spherical shape; i. e. the electric 
interaction energy (and therefore also the binding energy) should vary with the values of 
Iz elimination of the degeneracy). The fact, that this plausible effect is seemingly 
not observed experimentally, may be understood as a further hint at the "unreasonable" 
logic of the quantum world. 

Subsequently, these results are elaborated through the following arrangement: 
In Sect. II the fundamentals of RST are briefly put forward, with an emphasis on 
the two-particle Dirac equation for the fluid- dynamic quantum matter and on the (non- 
Abelian) Maxwell equations as the field equations for the matter interactions, cf. equa- 
tions (III.13P and (lII.lSp . This coupled system of Dirac and Maxwell equations is of 
such a logical structure that the standard conservation laws of charge (j^) and energy- 
momentum (T^jy) do emerge as an immediate and natural consequence, see equations 
(III.23P and f lll.32p below. The important point with such a logical structure is that a 
plausible definition of the total energy i?T, due to an RST field configuration, suggests 
itself: namely as the spatial integral of the time component Too of the energy-momentum 
density T^^, over whole three-space, cf. equation (1II.38P below. Such a definition of energy 
seems adequate for a fluid- dynamic theory, whereas in the conventional quantum theory 
(as a probabistic point-particle theory) the concept of energy is based rather upon the 
eigenvalues of some energy operator (H). 
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In any case, a definite energy can be expected only for a stationary system, sucli as 
positronium which is the object to be studied in the present paper. Here it is assumed that 
the corresponding gauge fields, mediating the interactions between the electron and the 
positron, are time- independent and the matter fields appear as the product of the usual 
exponential time factor (~ e~*^^) and a time-independent pre-factor ip{r), cf. flll.50|) . As 
usual, these assumptions convert the general RST field equations to the Poisson equations 
for the gauge potentials, cf. f lll.43p - flll.44p . and to the "mass eigenvalue equations" for 
the matter fields which we write down for the present purpose in Pauli form, cf. f lll.52ap - 
f]II.52bp . rather than in Dirac form. This coupled set of Pauli and gauge field equations 
is then our point of departure for working out the non-relativistic energy spectrum of 
para-positronium. 

Indeed, the RST treatment of positronium (i. e. the bound system of electron and 
positron) leads one to a certain observation which one encounters also in the conventional 
quantum theory: namely, this specific two-particle system occurs in two forms, i. e. ortho- 
and para-positronium, which we hereafter briefly term ortho/para dichotomy (Sect. III). 
However, this phenomenon of dichotomy is characterized in the conventional theory and in 
RST by means of quite different concepts, although both approaches do refer here to the 
particle spin. More concretely, the conventional theory resorts to the addition theorem of 
angular momenta for the spins of electron (sg) and positron (sp) which is realized by taking 
the tensor product of the single-particle spin spaces, such that the total spin 5* = Sg ± Sp 
is zero {S = 0) for para-positronium and S = 1 for ortho-positronium (sg = Sp = |). In 
contrast to this, RST is based upon the Whitney sum of single-particle vector (or spinor) 
bundles; and this means that each single particle does occupy a definite one-particle wave 
function which then carries individually the angular momentum of the composite system! 
Consequently, each positronium constituent (i. e. electron and positron) is described by 
a one-particle wave function due to either vanishing spin {s-p = Se — Sp = ^^para- 
positronium") or due to spin = 1 (= Sg + Sp), i.e. "ortho-positronium". In this 
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way, both manifestations of positronium appear to be equipped with integer spin (i. e. 
s-p = 0,so = 1) and thus are to be considered bosons. 

On the other hand, any single Dirac particle must be considered a fermion; and this 
entails of course the necessity to provide now in RST a clear definition of the bosonic and 
fermionic states for few-particle systems being composed of Dirac particles. 

Fermionic and Bosonic States in RST 

The resolution of this composition problem works as follows: the fermionic or bosonic 
character of a one-particle state is defined by the angular momentum eigenvalue of the 
selected spinor basis. Consequently, a fermionic basis (such as (]III.30a|) -( 1III.30d|) ) carries a 



half-integer value of angular momentum and a bosonic basis (such as (]III.35|) or ( IIII.36aP 



f lIII.36b|) ) carries an integer- valued angular momentum. Any Dirac spinor ip must now 
be decomposed with respect to such a spinor basis and thereby inherits its fermionic or 
bosonic character from just that selected spinor basis. By such a spinor decomposition 
there do arise as usual the corresponding spinor components ( ''amplitude fields" ) which are 
required to be always single-valued fields and therefore do always carry integer angular 
momentum, cf. f lIII.43a|) - (]III.43b|) . whereas the spinor basis will in general be multi- 



valued, cf. flIII.5naj) -f lITr50b|l . 



This multi-valuedness of the RST spinor fields does, however, not induce phathologi- 
cal elements in the logical structure because the observable objects of the theory (such as 
charge and energy-momentum densities) are always unique compositions of the non-unique 
spinor fields. For instance, the uniqueness of the charge density ^""^ko (r) as a bilinear con- 
struction of the spinor fields ipa{'>^ {ci = 1,2) ensures the uniqueness of the corresponding 
gauge potentials which feel the charge densities as their sources, see the Poisson equations 
f lIII.21|) -( ]III.22|) . In this way, one can arrive at a unique eigenvalue problem for the deter- 
mination of the positronium binding energy, although the underlying RST spinor fields 
are multi-valued, see the eigenvalue system (lIV.2ap -( IIV.2d|) below. 
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By means of all these rather technical preparations, one can now in Sect. IV tackle 
the proper problem of the positronium energy spectrum, where we are satisfied for the 
moment with the restriction to the para-case (-w Sp = 0). This problem is of course 
not exactly solvable so that we have to resort to various approximation assumptions. 
The first one of these refers to the non-relativistic situation; the next one concerns the 
electrostatic approximation which neglects the magnetic effects altogether. Finally we 
take the spherically symmetric approximation as our point of departure and furthermore 
treat the anisotropy effect only in the lowest order of approximation, i. e. the quadrupole 
approximation. Here, the assumption of "spherical symmetry" refers exclusively to the 
electrostatic gauge potential ■>^^Mo(r) (putting ^p^Aq^t) =^ [pUo(r); r = ||r||), not to the 
matter fields which are admitted to be anisotropic. Indeed, it is just this approximation 
assumption of spherical symmetry ^'P^AQ{f^ =^ ^^'^AQ{r) which allows us to separate exactly 
the original mass eigenvalue equation into an angular and a radial problem, where as usual 
the radial problem ultimately provides the wanted energy levels; but the angular problem 
is independent of the special interaction potential and thus may be solved exactly for all 
spherically symmetric interaction potentials. 

Angular Momentum Quantization in RST 

The angular problem is defined by the two coupled equations (lIV.8ap - (lIV.8bp below; 
and its solution is just what one may term ^^angular momentum quantization in RST" 
because it corresponds to the analogous phenomenon in the conventional quantum theory: 
there are two quantum numbers of angular momentum (termed i-p and iz) where i-p takes 
its values in the set of natural numbers {i-p = 0, 1, 2, 3, . . .) and iz is an integer whose range 
is —ip < iz ^ +ip. It is only the quantum number ip which explicitly enters the radial 
eigenvalue problem, cf. equations f lIV.9a[) - (]IV.9bp below, so that the second quantum 
number iz cannot immediately influence the energy spectrum! This entails the effect of 
degeneracy of the energy levels with respect to iz which occurs also in the conventional 
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theory (where the spherically symmetric Coulomb potential is adopted as the interaction 
potential between both positronium constituents). The various RST angular distributions 
of the current density, due to the possible values of i^, are displayed for £-p = 3 in the 
Fig.s IV. A on p.s A preliminary first analysis of this RST algebra of angular 

momentum is presented in App.D. 

Quadrupole Approximation 

This kind of /2-degeneracy would be also present in RST, but only if one is satisfied 
with the spherically symmetric approximation of the interaction potential! The corre- 
sponding result for the energy spectrum has been studied in a precedent paper [4], where 
it has been found that those "spherically symmetric" RST predictions undergo certain 
deviations from the conventional spectrum (\l.2\i ranging up to some 10% for principal 
quantum numbers n in the range 1 < n < 100. This appears to us as an acceptable 
result in view of the very rough approximation techniques being applied up to now within 
the RST framework. Especially, one becomes now curious whether perhaps the RST 
predictions for the positronium spectrum will come even closer to the corresponding con- 
ventional predictions (11.21 ) if those anisotropy effects are taken into account which are to 
be neglected in the very spirit of the spherically symmetric approximation. 

In order to clarify this question, one must first consider the anisotropy of the inter- 
action potential and then calculate the effect of this anisotropy on the binding energy. 
Since the gauge potential ^^^Ao{r,i!}) is tied up to the charge density ^^^ko^r,!!)) via a 
Poisson-type equation, cf. pV.10|) . and since furthermore the charge density ^^^ko is itself 
built up by the anisotropic wave functions, the latter objects do transfer their anisotropy 
directly to the interaction potential ^^^AQ(r,'&). Therefore it must be possible to express 
the anisotropy of the gauge potential ^''^^Ao{r,'d) in terms of the anisotropy of the wave 
functions (or amplitude fields, resp.); and indeed, this can be realized by expanding the 
anisotropic part of the gauge potential in terms of the angular pre-factors of the wave 
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functions, cf. equation (1IV.103|) below. The first non-trivial term of this perturbation 



expansion is the quadrupole term which appears as the product of an angular-dependent 
function ^^"'(^9) and a radial function ^^^'"(r), cf. fllV.1071) . 

The pleasant property of such type of perturbation expansion is now that the angular 
functions (such as ^^^A^^^d) for the quadrupole case) can be exactly traced back to the 
angular-dependent pre-factors of the wave amplitudes; and since the latter pre-factors are 
uniquely associated with the (discrete) values of angular momentum one thus obtains 
a unique link of the angular dependence of the gauge potentials to the various values of 
iz, see equation (lIV.104cP below. And it is now just this effect which brings about the 
elimination of the "spherically symmetric" ^^-degeneracy of the energy spectrum! Clearly, 
any specific angular dependence of the interaction potential (due to the various values of 
£z) must necessarily entail a different interaction energy and therefore also a different 
binding energy, see the estimate of quadrupole energies flIV.1291) due to the various values 
of iz for €-p = 3 (where the possible values of iz are iz = 0, ±1,±2, ±3). A certain 
peculiarity emerges here for i-p = 3, where it turns out that the quadrupole correction term 
for iz = ±2 is zero so that an improvement of the corresponding "spherically symmetric" 
energy level can arise not until the octupole approximation is considered, see equation 
fllV.104cp . 

Para-Positronium Spectrum 

But once the total energy (as sum of the kinetic and the electrostatic interaction en- 
ergy) is determined, one can (at least approximately) calculate the non-relativistic energy 
spectrum of para-positronium. As a nearby approximation procedure, one chooses some 
plausible trial configuration with included variational parameters for the concerned RST 
fields, cf. f lIV.131|) . and then one extremalizes the energy functional on that selected set 



of trial configurations. As an example, we discuss here explicitly the situation of eliminat- 
ing the ^2-degeneracy in the quadrupole approximation for £-p = 3, see Fig. IV. B, p. (HZl 
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The result is as expected; namely, the splitted RST levels become simultaneously shifted 
towards the conventional prediction (11. 21) . The RST level closest to the conventional 
prediction is due to 1^ = ±3. Its deviation from its conventional counterpart (11.20 is 
0,005 [eV] corresponding to a relative deviation of 1,2%. 

This motivates one to consider the whole spectrum for iz = ±£p in the presently 
discussed quadrupole approximation. The results are collected in Fig. IV. E on p. 11091 
The most striking effect of these results refers to the fact that the RST-splitted levels do 
form some kind of energy band around the (highly degenerated) conventional levels Eq 
( 11. 2p . For low principal quantum numbers n, the bands are well-separated; but it seems 
that for sufficiently large n the bands will overlap and thus constitute a rather intricate 
spectrum, rather different from the highly regular conventional spectrum (11.20 . It remains 
to be clarified whether such a result is merely an artefact being caused by all the applied 
approximations (especially quadrupole approximation) or whether the emerging band 
structure is also confirmed by the exact (hitherto unknown) solution of the RST eigenvalue 
problem. 
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II Positronium Eigenvalue Problem 



In order that the paper be sufficiently self-contained, it may appear useful to mention 
briefly some fundamental facts about RST. As its very notation says, the central idea is 
the Relativistic Schrodinger Equation which leads one in a rather straight-forward 

way to the Dirac equation ( ]II.13|) for few-particle (or many-particle) systems. In order to 
ultimately end up with a closed dynamical system for the fluid-dynamic quantum matter, 
one adds the (generally non-Abelian) Maxwell equations flII.15|) where this coupled system 
of matter and gauge field dynamics automatically entails certain conservation laws, such as 
those for charge (111.250 or energy-momentum (111.320 . This fundamental structure of RST 
is then subsequently specified down to the non-relativistic positronium system, especially 
to its para-form, with the main interest aiming at its energy spectrum. 

1. Relativistic Schrodinger Equations 

A subset of problems within the general framework of RST concerns the (stationary) 
bound systems. The simplest of those systems is positronium which consists of two 
oppositely charged particles of the same rest mass (M). The physical behaviour of its 
matter subsystem is assumed here to obey the Relativistic Schrodinger Equation 

thcv^^ = n^^ (ii.i) 

where the two-particle wave function \E'(a;) is the direct sum of the two one-particle wave 
functions ^ipaix) (a = 1,2) 

^(x) = V^i(x) ©V^2(x) . (II.2) 

The gauge-covariant derivative V on the left-hand side of the basic wave equation (III.ip 
is defined in terms of the u(2)-valued gauge potential as usual 

V^^ = d,^ + A,^ , {11.3) 
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or rewritten in component form 

D^^,=d^iP,-iAliJi (II.4a) 
D^^2 = ^,^^J2-^Al^2. (n.4b) 

Here, the electromagnetic four-potentials (a = 1, 2) are the components of the original 
gauge potential Afj, with respect to some suitable basis (a = 1, ... 4) of the u(2)-algebra 

A^{x) = A-(a;) = + B^{x) x - B;{x) x • (11.5) 

Here, the electromagnetic generators (a = 1, 2) do commute 

[ri,r2] = (II.6) 

and the exchange potential -B^ is put to zero B^{x) = 0) because the two positronium 
constituents (i. e. electron and positron) do count as non-identical particles. Recall 
that the exchange effects, being mediated by the exchange potential B^{x), do occur 
exclusively for identical particles so that B^{x) is inactive for the positronium constituents 
(-^ B^{x) = 0). Thus the bundle connection Afj,{x) flll.Sp becomes reduced to its u(l) © 
u(l) projection 

A,ix) A;:(a;)r,. (II.7) 

2. Dirac Equations 

For Dirac particles, which are to be described by four-spinors ipaix), the Hamiltonian 
in the Relativistic Schrodinger Equation (lll.ip obeys the relation 

T^'n^ = Mc\ (H.8) 

where is the total velocity operator and thus is the direct sum of the Dirac matrices 

r^ = (-7^)©7^. (H.9) 
16 



The mass operator A4 specifies the two particle masses (a = 1, 2) 

M = iM''Ta (11.10) 

and is required to be Hermitian {A4 = Ai) and covariantly constant 

V^M = 0. (11.11) 

This requirement is trivially satisfied for particles of identical rest masses (M^ = = M) 
since for such a situation the mass operator becomes proportional to the identity operator 

M = M 1 . (11.12) 



Thus the result is that, by virtue of the relation (1II.8|) . the Relativistic Schrodinger Equa- 



tion f lll.l|) becomes the two-particle Dirac equation 

ihcT'^V^'^ = Mc"^^ , (11.13) 

or in component form 

ihc 7^ z^i = -Mc^ iJi (II. 14a) 

ihc Y Df, %l)2 = Mc^ ^2 , (II. 14b) 

where the gauge-covariant derivatives (D) of the single-particle wave functions ipaix) are 
given by equations (ITL4a])-(ITL4b]). 



For identical particles, the Dirac equations flIL14ap - (III.14bp would couple both par- 



ticles much more directly since the exchange potential is generated cooperatively by 
both particles and simultaneously does act back on any individual particle which then en- 
tails the phenomenon of self-coupling. However, for the present situation of non-identical 
particles the coupling is more indirect: any particle does generate a Dirac four-current 
ka^i{x) {a = 1, 2) which is the source of the four-potential A'^{x) (see below). And then this 
four-potential of the a-th particle acts on the wave-function ipb{x) of the other particle 
{b 7^ a) as shown by equations flII.14a|) - flII.14b|) in connection with the gauge-covariant 



derivatives D f llL4al) -f llL4b|l 
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3. Maxwell Equations 

The bundle connection A^{x) fill. 51) is itself a dynamical object of the theory (just as 
is the wave function \E'(x)) and therefore must be required to obey some field equation. 
This is the (generally non-Abelian) Maxwell equation 



V^F^,, = -A'Kia,X (11.15) 

he 

Here, the bundle curvature is defined in terms of the bundle connection as usual, 
i.e. 

J'^u = V^A - + [A, A] . (11.16) 

For the present situation of non-identical particles, the connection becomes reduced 
to its (Abelian) u(l) © u(l) projection, cf. ( III.7|) . which then also holds for its curvature 



Decomposing here both the curvature J-'^^, and current operator J'^j^ with respect to the 
chosen basis of commuting generators (a = 1, 2) 

J't^u F;^Ta (11.18a) 
Ju i3\ra , (11.18b) 



one obtains the Maxwell equations flII.15p in component form as 



V'^F;, = 47ra,j,V (11.19) 

Since for the present Abelian situation the curvature components F'^^ (111.18ap are linked 
to the connection components (111.7^ as usual in Maxwellian electrodynamics (in its 



Abelian form) 

= - V^Al , (11.20) 
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the Maxwell equations f lII.19|) for the field strengths become converted to the 



d'Alembert equations for the four-potentials 

□ /l;: = 47^asJ^ (11.21) 
provided the gauge potentials do obey the Lorentz gauge condition 

W^A"^ = . (11.22) 
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4- Conservation Laws 



One of the most striking features in the description of physical systems is that both 
classical and quantum matter do obey certain conservation laws. For the presently con- 
sidered Relativistic Schrodinger Theory, as a fluid-dynamic theory, this means that there 
should exist certain local conservation laws, preferably concerning charge and energy- 
momentum. Moreover, these local laws should turn out as an immediate consequence of 
the basic dynamical equations, i. e. the Relativistic Schrodinger Equation flll.ip and the 
Maxwell equations fill. 150 . 

In this regard, a very satisfying feature of the Relativistic Schrodinger Theory is now 
that such conservation laws are automatically implied by the dynamical equations them- 
selves. In order to elaborate this briefly, consider first the conservation of total charge 
which as a local law reads 



But such a continuity equation for the total four-current can easily be deduced from 
both the matter equation (lll.ip and the gauge field equations (III.ISP : and this fact signals 
the internal consistency of the RST dynamics. First, consider the gauge field dynamics 
flll.lSp and observe here the identity 



which holds in any flat space-time. Obviously, the combination of this identity with the 
Maxwell equations (lII.lSp yields the following continuity equation in operator form 



V^j^ = . 



(11.23) 




(11.24) 



V'J. = . 



(11.25) 



Decomposing here the current operator in component form yields 








(11.26) 
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which furthermore simphfies to 



Wjl^Q (11.27) 



{a = 1,2 



under the Abehan reduction (jll.18ap -f lll.18bp . But when the individual Maxwell currents 



do obey such a continuity equation (]II.27p . then the total current 



2 



J, = E ("-28) 

a=l 

must also obey a continuity equation which is just the requirement (111.23p . 

On the other hand, we can start also from the matter dynamics flll.ip and can define 
the total current by 

= ^T^^ . (11.29) 

The divergence of this current is 

V^j^ = (V^) + ^T^, (V^) + ^ (P'T^) $ . (11.30) 

Here, one requires that the gauge- covariant derivative of the total velocity operator is 
covariantly constant 

V^T^ = , (11.31) 

and furthermore one evokes the Relativistic Schrodinger equation together with the Hamil- 
tonian condition f lll.8p which then ultimately yields again the desired continuity equation 
(III.23p . Thus, the conservation of total charge is actually deducible from both subdy- 
namics of the whole RST system; and this fact supports the mutual compatibility of both 
subdynamics (i.e. the matter dynamics (III.ip and the gauge field dynamics (III.lSp ). 

A further important conservation law does refer to the energy-momentum content of 
the considered physical system. Aiming again at a local law, one may think of a continuity 
equation of the following form 

W (^>r,,, = , (11.32) 
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where ^'^^^u is the total energy-momentum density, i. e. the sum of the Dirac matter part 
i^^^tMu) and the gauge field part {^'^^^u)'- 

(^>r,. = (^>r,. + (^>r,. . (11.33) 



Clearly, the validity of the total law (]II.32|) does not require an analogous law for the 



subdensities but merely requires the right balance of the energy-momentum exchange 
between the subsystems, i. e. 

(^>r^. = -v^ (^>r^. . (11.34) 

Indeed, the matter part has been identified as 

^""^M-^ = ^ [^T^{V,^) - {V,^)T^^ + ^T,{V^^) - {V^^)T,^] (11.35) 

and the gauge field part by 

= ^K^, (^F-^,Fi ' - \9,.Fl,F^^'^ , (11.36) 

where K^is is the fibre metric in the associated Lie algebra bundle. The (local) conserva- 
tion law 0II.32P comes now actually about through the mutual annihilation f lll.34p of the 



sources of both energy-momentum densities, i. e. 

(^>r^. = - (^>r^. = he f^, . (11.37) 

Obviously the sources of the partial densities ^^^^u and ^'^^^u are just the well-known 
Lorentz forces in non-Abelian form. 

It should now appear self-suggesting that the definition of the total energy (-Et) of 
an RST field configuration is to be based upon the time component '^"^^oo of the energy- 
momentum density '^'^^^^^ i.e. 



Et = j rfV(^>roo(r-). (11.38) 
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But since the total density '•"^•^T^y is the sum of a matter part and a gauge field part, cf. 



(III. 331) . the total energy (111.38^ naturally breaks up in an analogous way 



Eb + Eg , 



(11.39) 



with the self-evident definitions 



Eb = 



I 
I 



rf¥(^>roo(r-') 



(II.40b) 



(II.40a) 



But clearly, such a preference of the time component Tqo among all the other components 
Tfj^i, entails the selection of a special time axis for the space-time manifold. This then 
induces a similar space-time splitting of all the other objects in the theory, i. e. we have 
to consider now stationary field configurations which are generally thought to represent 
the basis of the energy spectra of the bound systems. 



In the present context, the notion of stationarity is coined with regard to the time- 
independence of the physical observables of the theory, i. e. the physical densities and 
the electromagnetic fields generated by them. In contrast to this, the wave functions 
do not count as observables and therefore are not required to be time- independent. But 
their time-dependence must be in such a way that the associated densities become truly 
time- independent . 



The simplest space-time splitting refers to the four-potentials A°'^, which for the sta- 
tionary states become time-independent and thus appear in the following form: 



5. Stationary Field Configurations 



Gauge-Field Subsystem 




(11.41) 



(a = 1,2). 
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A similar arrangement does apply also to the Maxwell four-currents j"^ 

J\ \^%ir^ ; -Ja{r)] , (11.42) 



so that the d'Alambert equations flll.2ip become split up into the Poisson equations, for 
both the scalar potentials '■"•*^o(r) 

A('^Uo(f) = -47ras (11.43) 

and the three- vector potentials Aa{r) 

AAa{f^ = -AwaJai^ . (11.44) 

Recall here that the standard solutions of these equations are formally given by 

WAo(rO = as f d^r' (II.45a) 

X(r) = a, I Sr' ^^^^ . (n.45b) 

Similar arguments would apply also to the exchange potential B^^ = j-Bo? ~-^| tiut since 
we are dealing with non-identical particles the exchange potential must be put to zero 

It is true, the particle interactions are organized here via the (electromagnetic and 
exchange) potentials which, according to the principle of minimal coupling, are entering 
the covariant derivatives D^ipa of the wave functions ipa as shown by equations (111.4ap - 



(]II.4bp . But nevertheless it is very instructive to glimpse also at the field strengths F"" . 



Their space-time splitting is given by 

E, = {('^)E^} = {nJ (n.46a) 
Ha = {^''^H^] = \)^e^\F\'^ (II.46b) 

and thus the linear Maxwell equations 0II.19P do split up in three-vector form (a = 1, 2) 
to the scalar equations for the electric fields 

V'K = ^T^a,''% (11.47) 
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and to the curl equations for the magnetic fields 



V X Ha = 47ras ja ■ 



(11.48) 



There is a pleasant consistency check for these linearized (but still relativistic) field equa- 
tions in three-vector form; namely one may first link the field-strengths to the potentials 
in three- vector notation (cf. flll.20p for the corresponding relativistic link): 



EJr) 



Hair) = V X AJr) 



(II.49a) 
(II.49b) 



and then one substitutes these three-vector field strengths into their source and curl 
equations (]II.47|) - (]II.48|) . In this way one actually recovers the Poisson equations (111.43^ - 
(1II.44I) for the electromagnetic potentials ^""^Aq, Aa- 



Matter Subsystem 

Concerning now the stationary form of the matter dynamics, one resorts of course to 
the generally used factorization of the wave functions ipa{r,t) into a time and a space 
factor 

M,c2 ^ 



ilJa{r,t) = exp 



-t 



(11.50) 



Here, the mass eigenvalues Ma {a = 1, 2) are the proper objects to be determined from the 
mass eigenvalue equations which we readily put forward now. For this purpose, observe 
first that the Dirac four-spinors tpa may be conceived as the direct sum of Pauli two-spinors 



(11.51) 



Consequently, the Dirac mass eigenvalue equations (to be deduced from the general Dirac 
equations flll.14ap - flll.14bj) by means of the factorization ansatz pi.50p ) are recast to their 
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equivalent Pauli form for the two-spinors ^°V±(^ 

(Observe here that we do neglect for the moment the magnetic effects by putting the 
three-vector potentials Aa{r) (III.4ip to zero: Aa{r) =^ 0). The mass eigenvalue for the 
positron (with rest mass Mp) is denoted by Mi and for the electron (with rest mass Mg) 
by M2. 

Summarizing, the RST eigenvalue system consists of the mass eigenvalue equations 
f lll.52ap - flll.52bp for the Pauli spinors ^"'^Lp±{f) in combination with the Poisson equations 
f lll.43p . Since the magnetic effects are neglected, the Poisson equations f lll.44p need not be 
considered here. But what is necessary in order to close the whole eigenvalue problem is 
the prescription for the link of the Pauli spinors ^""^Lp^^f) to the Maxwell charge densities 
^"'^jo{f), or more generally to the Maxwellian four-currents j""^ f lll.42p as the sources of 
the four-potentials A"-^, cf. the d'Alembert equations f lll.2ip . Surely, such a link between 
the wave functions ipa and the currents j"^ will have something to do with the Dirac 
four-currents which are usually defined by (a = 1, 2) 

kai, = iJaJulpa ■ (11.53) 

Indeed, a more profound scrutiny reveals the following link [7] 

J'V = = ^1 ^« ^1 (11.54a) 
j\ = -k2^. = -^2 7n . (11.54b) 

The change in sign of both Dirac currents reflects the positive and negative charge of both 
particles. In terms of the Pauli spinors ^""^Lp^^r) (a = 1,2) the space and time components 
of the Dirac currents read 

(r) = (-^Vl (r) (r) + V- (r) (r) (11.55a) 
4(f) = ^%\{r) a ^%_{r) + ^%l{r) a _{r) . (II. 55b) 
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Although the magnetic effects, which originate from the three-currents ka{r) via the mag- 
netic Poisson equations f lll.44p . are neglected in the present paper these current densities 
nevertheless play now an important part for identifying two essentially different kinds of 
positronium. 
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Ill Ortho/Para Dichotomy 



In the conventional theory, the manifestation of two principally different kinds of 
positronium is traced back to the two possibilities of combining the spins of the electron 
(e) and positron (p) |8]: if both spins Se and Sp add up to the total Spin S = 1 

So = Sp + Se = I (III.l) 

_ 1\ 
2 / ' 

one has ortho-positronium; and zero spin 

Sp = - Se = (III.2) 

yields para-positronium. This is the well known ortho/para dichotomy mentioned in any 
textbook on relativistic quantum mechanics. In contrast to this (generally valid) com- 
position rule for angular momenta, the ortho/para dichotomy in RST is based upon the 
(anti) parallelity of the Maxwellian three-currents jai'f^)- Here, it is assumed that both 
particles do occupy physically equivalent one-particle states ipa{r) (a = 1,2) in the sense 
that the Dirac currents (and therefore also the Maxwell currents) are either parallel or 
antiparallel. Thus we propose the following characterization of ortho-positronium [H]: 

A?i(f) = -hir) = h{r) (III.3a) 

Ji(^ = Mr) = jb{r) = h{r) (III.3b) 

A^{r) = A^i^r) = A,{r) (III.3c) 

Hi{r} = H2{f) = Mr) , (111.3d) 
and analogously for para-positronium: 

ki{r) = hir) = kpif) (III.4a) 

Ji(r) = -J2(r) = Jp(r-) = kp{r) (III.4b) 

A(r) = -A2(f) = Ap{r) (III.4c) 

Hi{r) = -Mf) = Hp{^ . (III.4d) 
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The interesting point with such a subdivision of positronium into two classes is the fact 
that this subdivision is based upon the magnetic effects which however are neglected for 
the present paper; but despite this neglection the subdivision is of great relevance also for 
the presently considered electrostatic approximation! Namely, even within the framework 
of the latter approximation scheme, there do emerge different distributions of electrostatic 
charge {^"'^ko{f^^ for ortho- and para-positronium with the corresponding quantum num- 
bers (-^ ortho/para dichotomy); and moreover there does arise also a certain ambiguity 
of the electric charge distribution even within the subclass of the ortho-configurations due 
to the same quantum number (-w ortho- dimorphism) . In order to elaborate these effects 
it is necessary to first specify the general eigenvalue problem down to the subcases and 
then to look for the corresponding solutions. 

1. Mass Eigenvalue Equations 

The hypothesis of physically equivalent states for both positronium constituents entails 
that both mass eigenvalues Ma are actually identical, i. e. 

Ml = -M2 = -M, . (III.5) 

Next, one considers the Maxwellian charge densities '•"-'jo(^) which must of course differ 
in sign for oppositely charged particles 

^')jo(r) = -(^)jo(r) (III.6) 

and this must be true for both ortho- and para-positronium. On the other hand, the link 
(lII.54a|) - (]II.54b|) of the Maxwellian currents j° to the Dirac currents ka^ shows that the 



requirement (1III.6P demands the identity of the Dirac densities ^"•^kQ{r), i.e. 

(%{f) ^ (%(f) , (III.7) 
or rewritten in terms of the Pauli spinors flll.55a|) 



%l{r) %+{f) + %l (f) %^ (f) = %l{f) %+{f) + (V- (r) %- if) . (III. 8) 
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However, the situation is different for the Dirac three-currents ka{r), since they may 
differ in sign, cf. ( 1III.3aP vs. (lIII.4a|) . Expressing the (anti) parallehty of both Dirac 



currents in terms of the Pauh spinors, cf. (]II.55bp . one requires 

(^VU^ + %l{r)a'^^^ip+{r) (III.9) 



where the upper/lower sign refers to the ortho/para case, resp. Both conditions OIILSp 
and fllll.Qp can be satisfied by putting 



(iV^(f) = %+{r) = ('/pV+(^) (III. 10a) 

(f) = i %_ if) = (''/^V- if) (III. 10b) 

where the upper/lower sign refers again to ortho/para-positronium, resp. 

It is true, the disposal flIII.10al) - flIII.10bl) satisfies both algebraic requirements fllll.Sp 



and fllll.Qp . but additionally there must be satisfied also a differential requirement: ac- 



tually, the spinor identifications f llll.lOap - fllll.lObp leave us with just one spinor field 



(i.e. ^^V±(^) for ortho-positronium and for para-positronium); and therefore 

both spinor equations (111.52ap - (111.52b|) must collapse without contradiction to only one 



spinor equation (either for '•''V±(^ or *'^V±(0)- In order to validate this requirement, 
we have to make a disposal also for the electrostatic gauge potentials *^"''ylo(r*). But this 
can easily be done by observing the link between the Dirac densities ^"'^^{f) and '■"''y4o(r) 
as it is implemented by the Poisson equations flll.43p . Indeed, this link entails that the 



potentials '^°''>AQ{f) must differ (or not) in sign when this is (or is not) the case also for 
the Maxwell densities ^"^jo(0- Therefore one concludes that for the positronium situation 
both electrostatic potentials must always differ in sign 

(lUo(r) = -('Uo(f) = (^/PUo(r) . (III.ll) 

But when this circumstance is duly respected, both mass eigenvalue equations ( ]II.52ap - 
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(]II.52bl) actually do collapse to a single one for ortho-positronium: 



ta-V ^%±{f) - (%(r) ■ ^%^ir) = c ■ %^ir) . (III. 12) 



The existence of one and the same eigenvalue equation (1III.12|) for both ortho-positronium 



constituents thus validates our original hypothesis that both the electron and the positron 
should occupy physically equivalent states. 

The relativistic pair (1111.12^ of Pauli equations has a single Schrodinger-like equation 
as its non-relativistic limit. Indeed, assuming that the "negative" Pauli-spinor ^^^(y9_(r) 
is always considerably smaller than its "positive" companion ^^^{p^{r) one can solve the 
upper one of the equations flIII.12p for ^^'>ip^ (r) approximately in the following form: 

^V(r1^^<?-V(V(r1, (III. 13) 

and if this is substituted into the lower equation (1III.12I) one finally ends up with the 
following non-relativistic eigenvalue equation of the Pauli form 



2M 



A ^%4f) - (%(f) • = ■ ^\+{r) . (III.14) 



Here, the non-relativistic eigenvalue emerges as the difference of the rest mass M and 
the relativistic mass eigenvalue M^,, i.e. 

E, = {M,-M)c\ (III. 15) 

Subsequently, we will be satisfied with clarifying the phenomenon of the ortho-dimorphism 
in the non-relativistic version (1111.14^ of the original relativistic eigenvalue equation 0111.12^ . 



But the case of para-positronium is a little bit more complicated. To begin with 
the positron equation flll.52ap . this becomes transcribed by the identifications fllll.lOap - 
(jlll.lObp . lower case, to the following form (para-positronium): 

i a • V (^V±(^ - ^''^^o(r) ■ ^^Vt(^ = ^*^^ c ■ ^''VtI^ • (HI- 16) 
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Obviously, this positron equation is just of the same form as the joint positron/electron 
equation (1111.12^ for ortho-positronium. However, the electron equation (]II.52b|) of para- 
positronium becomes transcribed by the identifications (lIII.10ap - (lIII.10bl) to a somewhat 
different form: 

ia-V (pV±(^) + ^"^Aoif) ■ (^Vt(^ = - ^*^^ c ■ ^P^if^ir) . (III.17) 

Since the sign of the potential term (~ *^^^y4o(r)) is reversed here in comparison to the 
positron equation (1111.16^ . the latter electron equation (1III.17P is the ''^charge conjugated'^ 
form of the first equation (lIII.lGp . 

The charge conjugation is defined here by the following replacements: 



'(^_(f) , ^P^ifi^ir) (111.18a) 



(pUo(r) -(^Uo(f) (111.18b) 
M, -M, . (III. 18c) 

Indeed, one is easily convinced that the two forms of eigenvalue equations flIII.16p and 
f llll.lTp for para-positronium are transcribed to one another by these replacements flIII.18ap - 
f lIII.18cp . This means that any solution of the positron equation flIII.16p can be interpreted 
also to be a solution of the electron equation (]III.17P : however, the charge conjugation 
does not leave invariant the Poisson equations, cf. f lIII.2ip - (]III.22p below. Therefore, we 
prefer here the use of solutions with the same non-relativistic limit! Indeed, it follows 
from both equations (IIII.16I) and (IIII.17I) that the "negative" Pauli spinor (r) can be 
approximately expressed in terms of the "positive" spinor (P^(y9_|_(r) through 



ih 



y^-[r) ~ ±— — — cr 



■2Mc" ^™-19) 
where the upper case refers to (1III.16P and the lower case to (lIII.lTp . This result may 
then be substituted in either residual equation (IIII.16P and (jlll.lTp which in both cases 
yields the same Schrodinger-like equation for the "positive" spinor ^^''ip^{r): 

- ^ A (^'V+(r) - /ic(PUo(f) ■ (^'V+(r) = ■ (PV+(f) . (III.20) 
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Thus it becomes again evident that the corresponding solutions of ( 1III.16I) and (]III.17|) 



do actually describe physically equivalent states. The notion of "physical equivalence" is 
meant here to refer in first line to the numerical identity of the energy being carried by 
anyone of the constituents of para-positronium, see below. 

2. Poisson Equations 

The mass eigenvalue equations do not yet represent a closed system and therefore can- 
not be solved before an equation for the interaction potential ^^/^''Ao(r) has been specified. 
On principle, this has already been done in form of equation (1II.43P so that we merely have 
to further specify that equation in agreement with the ortho/para dichotomy. Observing 
here the circumstance that the Maxwellian current of the first particle (a = 1, positron) 
agrees with the Dirac current, cf. flll.54a|l and flll.55a|l . the Poisson equation fill. 43 1) reads 



in terms of the Pauli spinors 

A (^/PUo(f) = -Ana, '^^'P^oi.r) (III.21) 

This Poisson equation closes the relativistic eigenvalue systems, both for ortho- 
positronium (1111.12^ and for para-positronium f lIII.16l) - flIII.17l) . For the non-relativistic 
limit, one merely suppresses the "negative" Pauli spinors (*/p)(y9_(r) so that the relativistic 
Poisson equations flIII.2ip simplify to 

A (^/^Uo(f) = -47ras ^''/^VU^ ^''^^V+(^ (III.22) 



which then closes both the non-relativistic eigenvalue equations ( 1III.14|) for ortho- 
positronium and (lIII.20p for para-positronium. 

Clearly, these coupled systems of eigenvalue and Poisson equations cannot be solved 
exactly (though exact solutions do surely exist), and consequently we have to resort to 
some adequate approximation procedure. But this suggests itself when we subsequently 
will establish the variational principle of minimal energy. 
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3. Non- Unique Spinor Fields 



Despite the fact that we originally subdivided the whole set of positronium configura- 
tions into two subclasses, i.e. ortho- and para-positronium (lIII.3ap - (lIII.4d|) . it may seem 



now that by the neglection of the magnetic forces we ended up with an eigenvalue problem, 
which does no longer offer any handle for sticking to that original subdivision into two 
peculiar subsets. Indeed, the adopted electrostatic approximation does admit exclusively 
an interaction force of the purely electric type (being described by the electric potential 
(''/p)Ao{r)), whereas the original ortho/para dichotomy was based upon the three- currents 
ka{r) as the curls of the magnetic fields, cf. f lll.48|) . As a result of this neglection of 
magnetism, the Poisson equation (lIII.2ip . or (lIII.22p . resp., holds equally well for both 
the ortho-configurations (b) and the para-configurations (p). But also the mass eigen- 
value equations, especially in their non-relativistic forms flIII.14p and flIII.20p are formally 



the same for ortho- and para-positronium. Does this mean that, through passing over 
to the electrostatic approximation, the difference between ortho- and para-positronium 
has gone lost? This is actually not the case because the different angular momentum (in 
combination with the hypothesis of the physical equivalence of both constituent states) 
leaves its footprint also on the electrostatic approximation. 

The crucial point here refers to the fluid-dynamic character of RST, as opposed to the 
probabilistic character of the conventional quantum theory. This entails that in RST the 
angular momenta of the subsystems cannot be combined (to the total angular momentum 
of the whole system) in such a way as it is the case in the conventional theory (-w 
addition theorem for angular momenta [ID])- More concretely: if we wish to insist on the 
viewpoint that the (observable) angular momentum of the considered two-particle system 
should emerge as the eigenvalue jz of the angular momentum operator J'z = Cz + S^, i. e. 

Jz ^6/p(^ = ^'/^^J. h ■ ^b/p{^ , (III.23) 
then the quantum number ^^^'''^jz due to the whole two-particle system must be carried 
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already by any individual particle! Namely, the wave function \l/;,/p(r) refers here to 
the two-particle system as a whole and thus, according to the RST philosophy, is to 
be conceived as the direct (Whitney) sum of the one-particle constituent wave functions 
^,(f) (a = 1,2) 

^,/p(f) = (^/fVi(r) © (''/^(r-) . (III.24) 

According to this sum structure, the total angular momentum is also the sum of 
the individual angular momenta 

Jz = ^'^1 © ^'^J. . (III.25) 

Therefore the eigenvalue equations (1III.23P for the total angular momentum jf^ can be 
decomposed as follows 

Consequently both one-particle spinors ?/'i(r) and 4j2{r) must obey the same eigenvalue 
equation as the total wave function flIII.24p . i. e. 



(^/^Vi (r) = ^^/f^j, h ■ (^/^Vi (r) (III. 27a) 

(2)J, (^/PV2(r) = h ■ (*/?'V2(r) . (III.27b) 

Furthermore, the Dirac four-spinors ipai'f^ can also be conceived as the direct sum of 
Pauli two-spinors ^"'^Lp_^_(r) and i-^- 

Mr) = ^"V+(0 © ^"V-(^) • (111.28) 

Therefore the same eigenvalue equation must also hold for the individual Pauli spinors, 
especially after those identifications Olll.lOap - piI.lObp : 



j{+) Wp)ip^{r) = h (III. 29a) 

Wp)^_(f) = iWp)j^ (^/pV_(f) . (III.29b) 
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This means mathematically that the eigenvalue ^^^^''jz of the two-particle state '^b/p ( 1111.23^ 



becomes transferred to any individual Pauli component of this two-particle state! In 
physical terms, the bosonic or fermionic character of the two-particle state \E' becomes 
thus incorporated in any individual constituent of the two-particle system. 

But now it is clear that positronium as a whole carries bosonic properties (-w is 
integer- valued) ; and this must therefore hold also for any Pauli constituent '^"•^ip±{f) of 
both particles (a = 1,2), cf. flIII.29ap -f lIII.29bl) . On the other hand, it is well known 
that the Pauli spinors do form a half- integer representation of the rotation group SO (3). 
This means that one can select in any two-dimensional spinor space a certain spinor basis 
{C'^e™} "with the following eigenvalue properties: 

/'C^r = j(j + i)^'C^T (ni.30a) 

Xc\'" = mhej^ (III.30b) 
L2(J;- = £(£+l)n2C^;™ (III.30c) 

§2 ^i, m _ ^ ^^^2 ^i, m _ (III.30d) 

Here the electron/positron spin is s = |; the orbital angular momentum is £ = 0, 1, 2, 3, . . . 



and thus the lowest possible value of j(= £ ± s) is j = 2 with £ = or £ = 1. Therefore 
we have two basis systems for j = \, namely \Q'^]Q' ^ \ and s Ci ; Ci^' ^ \. 

Fermionic States 

The existence of these two basis systems admits us to decompose now a fermionic 
state in the following way 



(4 



1 1 , s 11 



v+(ri = (''>;^+(r-) ■ cr + ^"^s+i^ ■ Cr~' (ni.3ia) 



('^V-(^) = {^''^-(r) ■ CP + ^''^5_(f) ■ Cf' ^} . (III.31b) 
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The action of the angular momentum operator on such a fermionic state is obviously 
= (4 ^">/^+(r1) ■ Cr + (4 ^"^'5+(f)) ■ Cr~' (ni.32a) 



4 + ^)(">^+(r) 



Co" 



2 



1_ 1 

■ 2 ' 2 



j(-) Wy,_(f) = (l, ('^>;^_(f)) ■ cr - ^ (4 (''^^^f)) • cf 



(III.32b) 



4 + ^]('^>;^-(r) 



1 1 

2 



1 _1 

2 ' 2 



The required results (IIII.29ap - flIII.29bp of the action of the angular momentum operator 
Jz on the Pauli spinors '^"•^ip±[f) (with a = 1,2 or a = h/p) are now deducible from the 
present equations f lIII.32ap - flIII.32bp by making the following arrangements: 

Lz (''^i {r)=izh- ('*>;^± (f) (III.33a) 
Lz (")5±(f) = (4 + l)h ■ (")5±(r-0 . (III.33b) 

Indeed, with these disposals the equations flIII.32ap - flIII.32bp adopt the required form 
of the eigenvalue equations f lIII.29ap - flIII.29bp with the eigenvalue of angular momentum 
being found as 



Jz 



+ 



1 



(III.34) 



Since the quantum number of orbital angular momentum is adopted as integer [i^ = 
0, ±1, ±2, ±3, . . .), we actually end up with half-integer quantum numbers ^"•'j^ flIII.34p 
for fermionic states! 

These fermionic states can obviously be realized by use of unique amplitude fields 
^"•^±(r) , ^'^^S±{f) and also unique spinor basis fields Co''^^, Ci^'^^ • Evidently, the latter 
fields work as the carriers of the spin, whereas the amplitude fields do contribute the 
orbital angular momentum. If this philosophy is tried also for the bosonic states we are 
forced to give up the uniqueness of the spinor basis! 
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Bosonic States 

Joining here the general conviction that bosonic states should have integer quantum 
numbers ^"^j, f lIII.29aj) - flIII.29bjl . i. e. ^"^j, = 0, ±1, ±2, ±3, . . ., it suggests itself to think 



that the amplitude fields should furthermore carry integer quantum numbers (f^) of orbital 
angular momentum; i. e. such equations as flIIL33ap - flIII.33bl) should persist also for the 
bosonic states. Thus the necessary modification must refer to the basis spinor fields (^'^"^ 
f lIII.30a|) - (]III.30dp . More concretely, we think of four basis spinor fields ^g^'', ^j^'', 



which for para-positronium obey the following eigenvalue equations 

Ji+) ej-) = Ji+) = Ji-' = Ji-' d"^ = ; (III.35) 
and similarly for ortho-positronium one wishes to work with a spinor basis 
rio^\ Vo'K vi^^ of the following kind: 

(III.36a) 
(III.36b) 
(III.36c) 
(III.36d) 

This says that for the para-case (1III.35I) the spins Sa (= |) of both constituent particles 
(a = 1, 2) do combine to zero spin quantum number sp of the para-type (sp = Si — S2 = 0, 
Si = S2 = |); whereas for the ortho-case (]III.36ap - (lIII.36dp the individual spins combine 















Ji- 




u {+) 









to unity {so = si + S2 = 1). Thus both basis systems f lIII.351) and flIII.36ap - flIII.36dp 



carry integer spin and therefore may be used for the corresponding decomposition of the 
Pauli spinors ^°'^Lp±{f) due to a bound two-particle system. 

For para-positronium {s-p = 0) one has now in place of the fermionic situation f]III.31ap - 
f ]III.31bp the following decomposition 



(^V+(r) = (^>7^+(rl ■ + ^"^'5+(f) • e^-^ (III.37a) 

-^|(^>;^_(f)■e^^^ + (^')5_(f)■e^^} . (iii.37b) 
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Here, the action of the angular momentum operator looks now as follows 



(-) 





(L,(^>;^+(r-))■e^+) + (L,(^)5+(f))■e 



(III.38a) 
(IIL38b) 



Consequently in order to satisfy again the para-form (p) of the eigenvalue equations 
(llll.29al) - (IIII.29bp one puts here 



(IIL39a) 
(IIL39b) 



so that the quantum number ^^^jz 0III.29ap - piI.29bp is solely due to orbital angular mo- 
mentum: 



(p). 



Jz 



(III.40) 



(4 = 0,±1,±2,±3,...) 



For ortho-positronium {sq = 1), the situation is somewhat different. First, the de- 
composition of the ortho-spinors ^^V±(^) looks quite similar to the para-case flIII.37ap - 
flIII.37bD : 



V.+ (r-) = (^>;^+(r-)■r^W + ('')5^.(f)■r^(-) 



ibyj^_ 



(f) ■ r^l^-' + 



(III.41a) 
(III.41b) 



But the action of the angular momentum operator Ji^"* on these ortho-states looks now 
as follows 



j(+) (b) 



■ 4"-' + 



(III.42a) 



J(-)(V(r) = -^|^f4 + /l)(^>7^_ 



■ 4^' + 



(L, - (^)5_(f) ■ r^S^^I . (III.42b) 
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For satisfying here again the eigenvalue requirement ( 1III.29a|) - (lIII.29b|) in its ortho-form 
(6) one puts in a self-evident way 

L, (^>;^± =i,h- ^^>R± (III.43a) 
L, ^^'^S± = {i, + 2)h- (^)5± . (IIL43b) 

This arrangement yields then the following eigenvalue equations for the ortho-spinors 

J± W^i(f) = (i^ + l)h- (^V±(r) . (III.44) 

Thus one finds the quantum numbers of the ortho-system: 

^'y. = 4 + 1 (III.45) 
= 0,±1,±2,±3,...). 

For a realization of the required basis spinors one takes the fermionic ba- 
sis < Co ' ^ , Ci ' M as the point of departure and introduces a general spinor basis 

through 



Cr (ni.46a) 

J-) ^ ,Q^--2 (III.46b) 

^(+) _ g-*6<^ .^il (III.46c) 

c^i") =e*^-Cf~^ . (IIL46d) 
Here it is easy to see emerging the following relations concerning angular momentum 

jWa;W = -(6-0 h.J^^ (III.47a) 

jW4-)= n-J,^ (III.47b) 

jWc.W = -(?)-l);i.c.W (111.47c) 

Ji-) c.;-) = (^?. - h-J-\ (III.47d) 

40 



Obviously, this w-basis contains some free parameter S (i. e. the boson number) and if 
this is chosen as S = | we obtain the desired ^-basis (1III.35I) for para-positronium; and 
if S is chosen as 5 = — | we obtain the r^-basis (IIII.36ap - (]III.36d|) for ortho-positronium. 



For S = we get back the purely fermionic C-basis (lIII.30ap - (]III.30dp . 

With the choice of the w-basis flIII.46a|) - flIII.46d|) the loss of uniqueness becomes now 



evident: since the original ^-basis flIII.30a|) - flIII.30dl) is unique (i.e. the C,^'J^{'d,(f)) consti- 
tute a "unique" spinor field on the 2-sphere S*^), the other two basis systems ^ flIII.35p and 
Tj ( 1III.36ap - (lIII.36dp are double-valued. More concretely, for both the ^- and the ?7-basis 



(^-w 5 = ±1) one finds by performing one revolution around the z-axis (0 < < 27r): 

eg) (0 + 27r) = e±- -ej J (0) = J (III.48a) 
r^S^ (0 + 27r) = e±- -r^g^ (0) = -r/g^ , (III.48b) 

and this says that we need two revolutions around the ^-axis (0 < < 47r) in order 
to return to the original basis configurations. Since we adopt all the amplitude fields 

unique scalar fields, the double-valuedness of the para- and 
ortho-basis becomes transferred to the para-spinors ^^V±(r) flIII.37ap - flIII.37bl) and ortho- 
spinors '•^V±(r) f|III.41ap - flIII.41bp and from here ultimately to the Dirac spinors ^v,o 



m^{r) = (pV+(0 © ^^V-(^) (III.49a) 
^^(f) = (^V+(^ © (^V-(^) • (III.49b) 

Summarizing, the Dirac wave functions \l/p(r') and \l/c)(r*) for ortho- and para-positronium 
must in RST be double-valued in the following sense: 

^p(r, ^9, + 2n) = -^p(r, ^, 0) (111.50a) 
^o(r, ^9, + 27r) = -*o(r, ^, 0) . (111.50b) 
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Uniqueness of the Physical Densities 

Naturally, in a fluid- dynamic theory (such as the present RST) the proper observable 
objects are the physical densities, such as those of charge, current, energy, linear and 
angular momentum, etc. A plausible condition on these densities is surely given by the 
demand that these physical densities should be single-valued tensor fields. But, as we will 
readily demonstrate, this condition can be satisfied also by non-unique wave functions; 
and this fact allows us to actually deal with such non-unique wave functions, as given 
for example by the double- valued positronium states flIII.50ap - flIII.50bl) . Therefore one 



wishes to have some condition on the wave functions which on the one hand admits their 
non-uniqueness but on the other hand ensures the uniqueness of the associated physical 
densities! 

Now according to the present RST philosophy, the non-uniqueness of the (Dirac) 
wave functions is to be traced back to the spinor basis, whereas the amplitude fields are 
furthermore required to be unique. Therefore the non-uniqueness of the wave functions 
is measured by the boson number S, cf. flIII.46ap - flIII.46dl) : and thus the condition of 
uniqueness of the densities is to be retraced to some condition for the fixation of the 
boson number 6. Such a fixation may be attained now by considering specifically the 
Dirac density ^""^Icq (f) and the three-current kai^^ which read in terms of the Pauli spinors 
^°V±(^ as shown by equations flll.55ap - fjll.55bp . Decomposing these Pauli spinors with 
respect to the w-basis flIII.46ap - flIII.46dp lets then appear the (Dirac) charge densities 



flll.55ap in the following form: 

4:71 

Evidently, these charge densities are unique in any case and therefore do not yet provide 
an immediate handle for fixing the parameter S. 

This situation changes now when one considers also the Dirac currents ka flll.55bp . 
which by their very definitions are always real-valued objects: 
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(^\^ = 1- . ^^yR_ + ^^)*s+ ■ ('^^^^ - (">^_ • (">;^+ - ^""^h^ ■ (111.52a) 

^""^k, = |e2^(^-^)<^ ■ C(,) - e-2^(^-^)* ■ C(.)| (IIL52b) 

W/t^ = i-z. ) i-yji^ . (-}jz_ + {-}jz_ . i-yjz^ _ . _ . I (111.52c) 

But here a nearby restriction upon the parameter S suggests itself, namely through the 
plausible demand that the Dirac currents ka ( ]II.55bp . with their components being spec- 
ified by 0III.52ap - (IIII.52c| . must be unique (!) vector fields over three-space (albeit only 
apart from the origin r = and the z axis = 0, vr). Evidently this demand of uniqueness 
reads in terms of the spherical polar coordinates {r, i?, 0} 

fc„(r,^?,0 + 27r) = 4(r,^9,0), (III.53) 

and thus the values of S become restricted to the range 

S = ^(n+1) (III.54) 

{n = 0,±1,±2,±3,...) 

which then entails also (half- jinteger quantum numbers for the z component of the angular 
momentum flIII.47ap - flIII.47dp of the spinor basis: 



= ± - 1 j = ±^ . (III.55) 

Here it will suffice to admit for the spinor basis (of the two-particle systems) exclusively 
the values 5 = ±i; other values of S come into play for bound systems of more than two 
fermions. 

Notice here that this (half- jintegrity arises as a consequence of the demand of unique- 
ness with respect to certain physical densities (i.e. Dirac current), whereas the corre- 
sponding integral quantum numbers of conventional non-relativistic quantum mechanics 
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are mostly traced back in the textbooks to the uniqueness requirement for the wave func- 
tions themselves (not the densities). The lowest values of Sz f llll.SSp are Sz = ±i for 6 = 
and Sz = 0, ±1 for S = ±i. Thus for the first case (5 = 0) we have a fermionic basis and 
for the second case (5 = ±|) one deals with a bosonic basis. In this sense, a Dirac particle 
is said to occupy a fermionic quantum state if the boson number'^ 5 of its spinor basis 
is zero {h = 0), and a bosonic quantum state if the boson number S equals ±|. Observe 
that through this arrangement the fermionic or bosonic character of the quantum state 
of a Dirac particle is defined by reference to the corresponding spinor basis. 
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IV Para-Positronium (5 = |) 



Though the best perspective on the ortho/para-dichotomy would surely result from a 
detailed synopsis of both configurations, it is nevertheless very instructive to first consider 
the case of para-positronium separately. The interesting point with the latter kind of 
positronium is here that the situation is much simpler than for the ortho-counterpart, so 
that a treatment of the para-case appears as an ideal preparation to the later discussion 
of ortho-positronium which itself decays into two subcases {^^ dimorphism"). Furthermore, 
the simpler case of para-positronium is also well-suited in order to present the quantization 
of angular momentum in RST from a more general viewpoint. 

1. Mass Eigenvalue Equations in Terms of Amplitude Fields 

Once the para-spinors *-''V±(^ are decomposed with respect to a bosonic basis, cf. 
(lIII.37ap - (]III.37b|) . one can substitute this decomposition into the two original eigenvalue 



equations in spinor form, cf. flIII.16|) - flIII.17p . and will then obtain the corresponding 
eigenvalue equations in terms of the amplitude fields ^^^±{r) and ^^^S±{f). But since 
both eigenvalue equations are assumed to yield physically equivalent states, there must 
exist a peculiar relationship between both eigenvalue equations. It is this relationship 
which first must be worked out. 

The identification (lIII.10ap - (]III.10b|) of the positron spinors ^^V±(^) with the para- 
spinors ^^V±(^) reads in terms of the amplitude fields 

(^yjz^(f) = (p>;^±(f) (IV. la) 

(i)5±(f) = (^')5±(f) (IV. lb) 

so that the spinor form 0111.16^ of the eigenvalue equation reappears in terms of the 
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amplitude fields as 



dr 



+ 



1 



r d(j) 2r 
1 d cot ^ 

r dd r 



M + M, 

h 



2r + u - 



_ 1 9 ^P>JZ+ _ cot 7? 
r r 



d(f) 



M + M, 



h 



dr r d(j) 
1 d cot ^ 

r r 



r 



2 ^ d<P 



M-M, 



h 



- c ■ (p>;^. 



ar r aa) r 



1 (9 (P>;^_ cot ?9 
r r 



M-M, 



h 



- c ■ ^p^s. 



(IV.2a) 



(IV.2b) 



(IV.2c) 



(IV.2d) 



On the other hand, those identifications Olll.lOap -f llll.lObl) of the electronic spinors *-^V±(^ 
with the para-spinors ^^V±(r) read in terms of the amplitude fields 



(2)7^^(r) = -^(p>;^±(f) 

(2)5^ (f) = -t^P^S±{r) 



(IV.3a) 
(IV.3b) 



so that the spinor form (1111.17^ of the electron equation reappears now from (]IV.2ap - 
flIV.2dD by means of the replacements flIII.18aj) - flIII.18cj) . 

In order to get some information about the peculiarities of this para-system (lIV.2ap - 
flIV.2dl) one first recalls the hypothesis that the spin Sp of the constituents should be zero, 
cf. flIII.351) . so that their total angular momentum is of purely orbital nature, cf. flIII.40p . 
Since, according to our hypothesis, the orbital part of total angular momentum is carried 
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by the amplitude fields, we may fix their general form as follows: 



a/ r sin {} 



(IV.4a) 
(IV.4b) 



a/ r sin 

If this form is substituted into the present eigenvalue equations flIV.2ap - flIV.2dp . this 
system reappears in terms of the real- valued amplitudes ^P^R±{r,'d) and ^P^S±{r,'d) 



d_ 

dr 



d_ 



(IV.5a) 



M + M, 

h 



■ ^P^R. 



^^'y 'y 'y 



pp)S_ + LJl . ip)s_ + (p)Ao ■ + - 

yy^y 'y 'y 



d_ 



(IV.5b) 



M + M, 

h 



(IV.5c) 



M-M, 

h 



■ ^P^R^ 



^{P)R 



-4 cot^9- 



(IV.5d) 



This eigenvalue problem may perhaps look somewhat complicated, but the assumption 
of a spherically symmetric potential ^P'^A^ir) (-w spherically symmetric approximation) 
admits us to separate it into an angular and a radial problem. For such a separation, one 
first resorts to the following product ansatz 



^P^R±{r,§) = fR{d)-'R^{r) 
^P^S±{r,^) = fs{^)-'S±{r) 



(IV.6a) 
(IV.6b) 
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and furthermore one makes the following identifications 



'R+ir) = 'S+{r) = <l+(r) (IV. 7a) 

'R^r) = 'S4r) = <l_(r) , (IV. 7b) 



which then ultimately separates the mass eigenvalue system flIV.5a|) - flIV.5dl) into an an- 
gular problem 

- 4 cot^ ■ fni^) = (/p + 4) ■ fsi^) (IV.8a) 
+ 4 cot ^ ■ fsi^) = (4 - Q ■ fni^) (IV.Sb) 

and a purely radial problem: 

- ^ • <l+(r) - (^Uo(r) ■ I._(r) = ^[±^1^. |,_(^) (iv.9a) 
dr r n 

%M + llln . ^_^r) + (PUo(r) ■ <l+(r) = ^^^^ c ■ l>+(r) . (IV.9b) 
dr r n 

Observe here that the eigenvalue Iz of orbital angular momentum does not enter the 
radial problem (lIV.9ap - (]IV.9b|) which ensures the occurence of energy degeneracy within 



the set of solutions belonging all to the same quantum number l-p but to different values 
of £2, provided the interaction potential '^^)Ao(r) is fixed from the outside. However, for 
our present two-body problem the interaction potential must obey the Poisson equation 
f lIII.2ip . which for the present case of para-positronium adopts the following form: 

A (PUo(rO = -47ras (^'A;o(f) (IV. 10) 

= -tts • 

sm v r 

Therefore the interaction potential *^^^Ao(r) depends upon the quantum state of positron- 
ium, i.e. more concretely: ^^^Ao{r) will in general depend also upon the quantum number 
iz of angular momentum; and this dependence will then be transferred also to the binding 
energy of the quantum state. Thus, it is just this effect which lets us suppose that the 
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angular momentum degeneracy of the conventional prediction (11. 2|) will be lifted in RST. 
Clearly, it remains to be clarified whether this lifting is (or is not) a negligibly small effect. 

But in any case we have arrived now at the definite form of the energy eigenvalue 
problem for para-positronium which (in the electrostatic and spherically symmetric ap- 
proximation) evidently consists of the angular equations flIV.8a|) - flIV.8b|) . the mass eigen- 
value equations flIV.9a|) -f lIV.9b|) and the Poisson equation fllV.lOl) . Here it should be clear 



that exact solutions of this eigenvalue problem cannot be found so that we will have to 
be content with the elaboration of approximate solutions which, however, are expected 
to display some of the qualitative features of the unknown exact solutions. The central 
point of our approximation procedure refers now to the solution ^P^AQ{r) of the Poisson 
equation (llV.lOp which may formally be written as 

^'%ir-) = a^ [ rfV'I^Mfl (IV.ll) 



47r J sin {}' J \\f—f'\\ 

From here it is easily seen that the interaction potential ^^^Ao{r) adopts the Coulomb form 
in the asymptotic region (r — > oo) 

lim(PUo(f) = — , (IV. 12) 



r 



provided the Dirac charge density ^P^ko{r) is normalized to unity 

d^f^P^oir) = 1 . (IV. 13) 

Since the density '^P^fco(r) is factorized, cf. fllV.lOp 

(p)fco(f) = ^''^koi^) ■ ^P^fco(r) , (IV. 14) 

the normalization condition ( IIV.ISP can be passed on separately to the longitudinal and 
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radial factors: 



J Att sm 

j dr r^^P^koir) = J dr r l>^(r) + $2 (r) =1. (IV.lSb) 




2. Angular Momentum Quantization in RST 



Obviously, both angular equations flIV.8ap - flIV.8bp represent a closed eigenvalue prob- 



lem for the quantum numbers i-p and iz, i-e. the quantization problem of angular mo- 
mentum in RST. The solution of this problem must yield the possible values of i-p and 
£z which then does apply quite generally for all bound 0(3)-symmetric systems, since 
the special form of the interaction potential ^^^AQ{r) does not enter the angular problem 
at all! Clearly, this is the RST counterpart of the corresponding eigenvalue problem in 
conventional quantum mechanics [10] 



L'^ \i,m) = fi^iii+l) \i,m) (IV.16a) 
Lz \i,m) =mh \i,m) . (IV. 16b) 



For a closer inspection of the coupled angular problem f 1IV.8a|) -f ]IV.8b|) it is very 



instructive to decouple those equations by differentiating once more which then yields the 
following decoupled second-order equations 

+ fp ■ fsi^) - 4(1 + 4) ■ ^ = . (IV.17b) 
dv'' sm v 

Here, the first striking item refers to the fact that the substitution iz =^ —iz entails the 
replacements /^(■i?) =^ fs{^), fs{^) =^ ~fR{^)- Of course, this pleasant property of the 
solutions fni^), fs{^) could have been deduced also from the original equations f lIV.8ap - 
f lIV.8bp . We will take advantage of this symmetry by solving the angular system only for 
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a restricted number of possible values of £2 (e. g. > 0) and will deliver subsequently 
the solutions for the other ^z by the above mentioned replacements. 

Concerning now the non-singular solutions of the decoupled system (IIV. 17ap -( 1IV. 17b"|) . 
it suggests itself to try an ansatz in form of a finite power series expansion with respect 
to sint?, i. e. we put (x = sin??) 

/k(^) ^ Fr{x) (IV.lSa) 
Fs{x) (IV.lSb) 

with 

Fn{x) = J2pnX'' (IV.19a) 

n 

Fs{x) = Y,^nX^ . (IV.19b) 



n 



The range of the powers (n) remains to be determined by substituting this ansatz back into 
the second-order system (IIV. 17ap - (]IV. 17bp which reappears now through the alterations 
fllV.18al) - (IIV.18bp in the following form 



,2- 



<PFr{x) dFnix) , ,2 ^ / N . .X Fb.{x 



X — X 



+ 4 ■ Fr{x) - 4(4 - 1) • = (IV-20a 



dx"^ dx ^'^ 



x^ 



^.d^Fsix) dFsix) , ^2 ^ / X 4(1 + *., 



(1 - x') ' - X — + t'p ■ Fs{x) - , ■ Fs{x) = . (IV.20b) 

ctx^ ax x"^ 

The solutions hereof are then obtained via the recurrence formulae 

^"^^ = (n + 2)(n + l")- 4(4-1) • ('^-2^") 

^2 _ £2 

= („ + 2)(n + l)-\(4 + l) ■ ■ 
In order that both series (IIV. 19ap - (]IV. 19bp remain finite one obviously has to impose the 
halt condition 

r^max = (IV.22) 

{ir > 0) 
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which fixes i-p as an integer [i-p = 0, 1, 2, 3, 4, . . .). Furthermore, the starting power n„ 
is found as 



and 



1 + 



1 + 



for flIV.21aD 



(IV.23) 



for flIV.21bD 



(IV.24) 



Therefore 4 is also found to be an integer (— £p < < ^p)- Thus the desired angu- 
lar functions fni^) and fs{'&) are ultimately found in terms of powers of trigonometric 
functions. 

Perhaps it is most instructive to consider a simple example {i-p = 3, say). First, taking 
iz = simplifies both equations f lIV.17ap - flIV.17bp to 



2 +9fs{^) = 



(IV.25a) 
(IV.25b) 



with the obvious solutions 



/«(^) = W- cos(3i9) 



— sin(3'i9) . 



(IV.26a) 
(IV.26b) 



This solution satisfies the first-order system flIV.8ap -f ]IV.8bp for ip = 3,- 



dd 
d^ 



-3/ij(^) 



(IV.27a) 
(IV.27b) 



and moreover it satisfies also the angular normalization condition f llV.lSap . 

Because of its untypical simplicity, the present solution flIV.26ap -f lIV.26bp for 4 = 
0,ip = 3 could be guessed directly from the second-order system f lIV.25ap -f lIV.25bp or 
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also from the first-order system (]IV.27a|) - (]IV.27b|) . But of course it must be possible 
to find it also by following our systematic procedure being based upon the formal solu- 
tions (IIV. 19ap -( ]IV. 19b"|l ■ These may be written down beyond the present example more 
generally as 



( h 



Fr{x) = { 



(i^ - i^) =^ even; 4 > 



X' 



Iv - 



odd; 



< 1 



Fsix) = I 



(Jyi X , 



CT^ X 



{ip - 4) ^ even; 4 < 



> -1 



(IV.28a) 
(IV.28b) 

(IV.29a) 
(IV.29b) 



(i^ - 4) ^ odd; 

Because of i-p = 3 and 4 = 0, we have to start here with either the second-order solution 
Fji{x) ( 1IV.28b|) and then have to construct fs{^) from the first-order equation (llV.Sap . 
or one can start also from Fs{x) ( ]IV.29bp and then one determines fni'^) from ( jIV.Sbp 

For this situation [l-p = 3,£z = 0) one finds from the recurrence 
starting, say, with Fs{x) (IIV.29bp 

4 



a3 = -3^i 



formula f lIV.21bj) by 
(IV.30) 



so 



that the desired second-order solution Fs{x) flIV.29bp becomes specified to 



Fsix) 



cri (x-^x-^ 



(IV.31) 



(£p = 3,4 = 0) . 
This then yields for the corresponding angular function fs{'&) f llV.lSbp 

/ A \ 



= io-isin(3?9) . 



(IV.32) 
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Of course, this function fsi^^) must satisfy both the second-order equation (lIV.17b[l and 



the first-order equation (]IV.8b|) from which we deduce the first angular function fni^'d) as 

fj^l'd) =^ |cos?9 - ^ cos^??! (IV.33) 
= — - (Ti cos(3't?) . 



Finally, by means of the angular normalization condition f llV.lSaP the constant (Xi 
becomes (up to sign) fixed to 

^1 = -3^1 , (IV.34) 
so that our method of decoupled second-order equations with the results (lIV.32p - (lIV.33p 
actually reproduces the immediate first-order result (lIV.26ap - (]IV.26bl) ! (Would we have 
started with Fr{x) fllV.28bp ) in place of Fs{x) fllV.29bl) . we had found both functions 
/r(^^)5/s(^) being interchanged in that result flIV.26ap - fllV.26bp ). Clearly, this equips 



us with sufficient confidence into the second-order method fllV.20ap -f llV.20b"|) in order to 
construct by means of it the higher angular momentum states due to iz = ±1, ±2, ±3. 

For instance, for = 1 (and still l-p = 3) we have to start now from the solution 
f llV.28ap of the second-order equation fllV.20ap . i.e. more concretely 



Fr{x) = pix + p3X^ (IV.35) 

with the coefficient pa to be taken from the recurrence formula flIV.21aP as 

P3 = -^Pi. (IV.36) 

Thus the second-order solution -Fr(x) becomes explicitly 

4 



Fnix) = p,\^x-^x^] (1V.37) 
(£-p = 3,£z = 1) 



and consequently the associated angular function fn^i^) fllV.18aP must look as follows 

fni^) = pi (^in^ - ^ sin^ ^ . (1V.38) 
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As a check, one is easily convinced that this function actually obeys the second-order 
equation (lIV.17a|) . Furthermore, one substitutes this angular function (1IV.38I) in 

the first-order equation (llV.Sap and thus finds from here the associated angular function 

fs(^) = -lp^cos^ sin^ ^ (IV.39) 

which, of course, must then obey the second-order equation (lIV.17bp . Finally, one deter- 
mines again the residual coefficient pi by means of the angular normalization condition 
f llV.lSap 

[24 

Pi = \-, (IV.40) 

V IT 

which puts the desired solution into its final form 

fj^{^) = sin ^9 - ^ sin^ {IV AU) 

2 /24 

/5(^) = _ W — sin^t^cos?? (IV.41b) 

3 V vr 

{£r = 3,4 = 1) . 

The residual cases for = — 1, ±2, ±3 and i-p = 3 can be dealt with in the same way and 
the results are collected in the following table. Observe here also the realization of the 
symmetries fni^) =^ fsi''^), fs{^) =^ —fR{^) being included by =^ —^z, as was 
already remarked below equation flIV.17bl) ! (For a more systematic computation of the 
angular eigenfunctions see App.D). 
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4 









cos(3?9) = (-3 COS ^9 + 4 cos^ ?9) 


- sin(3i9) = - y^ (3 sin ?9 - 4 sin^ ?9) 


+ 1 


■ (sin — 1 sin^ -i?) 


— 1 /l^ ■ sin^ -(9 cos 


- 1 


— ■ sin^ -i? cos ■& 


— y^ ■ (sin — 1 sin^ -i?) 


+ 2 


a/I^- sin^i^cos^? 


-sin^T? 


-2 


-sin^^? 


^/l^- sin^'i^cos^? 

Y 3"" 


+ 3 







-3 








As pleasant as this RST method of angular momentum quantization may appear, 
there is also a somewhat critical point which refers to the case = 0. Indeed for this 
case of vanishing iz, there arises an additional singularity on the z-axis! This is readily 
realized by considering the Dirac current kp{r) (llll.4bp which may be decomposed with 
respect to the spherical polar coordinates {r, {}, 0} as usual 

kp = kj. • Bj. -\- k^ • e^ -|- kfp • e^ , 

where the components ^^^kj., ^^^k./), ^'^^k^ have already been specified in terms of the ampli- 
tude fields by equations ( 1III.52a|) - (lllI.52cp . Substituting therein the present parametriza- 



tion of the amphtude fields 'R-±,S± by means of the angular functions fni^), fs{^) (IIV-6ap - 
f lIV.6bp and the radial fields $±(r) (lIV.7ap -( ]IV.7bp . one ultimately finds that only the 



azimuthal component is non-trivial and factorizes to a product form with respect to the 



56 



spherical polar coordinates, i. e. one finds 



^P^krir,^)=0 (IV.42a) 
^P^k4r,^)=0 (IV.42b) 
^P^k^{r, ^) = ^f^k^i^) ■ ^''^k^ir) , (IV.42c) 



with the azimuthal factors being given by 



r 



Now substituting here in the angular part (lIV.43aP the solution ( 1IV.26aP - (]IV.26b|) for 
4 = yields 

= (IV.44) 

and this result remains finite on the z-axis {i} = 0,7r). Clearly, this signals the presence 
of a magnetic singularity on the 2;-axis, cf. the curl equation ( 1II.48p . Thus the conclusion 



is that for the states with £^ = the presently adopted electrostatic approximation is 
possibly of limited use, and magnetism should be taken into account. 

From this reason, it seems worthwile to consider also the other states with non- 
vanishing iz- Intuitively, one expects that the centrifugal forces will press matter off 
the axis of rotation (i. e. the 2;-axis; = 0, vr) so that the (azimuthal) current density 
^^^kfj,{'d) becomes zero on the z-axis and adopts its maximal value the farer away from this 
axis the larger is the value \iz\ of angular momentum. For instance, consider the largest 
possible angular momentum (i. e. £z = ±3 for £-p = 3) and deduce from the above table 
that the corresponding azimuthal part of the current density ^^^k^p^i!}) (lIV.43aP becomes 

Obviously, this current density vanishes like ~ when approaching the z-axis ("i? — >■ 0). 
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Next, turn to the case £z = ±2 and find again from tlie corresponding results in tlie 
above table for the (angular) current strength ^^^k^f.^'d) (]IV.43al) 

= ±^ sin^ ^ 1 1 - ^ ^i^' ^ f • (1^-46) 

Comparing this to the precedent result (]IV.45|) for = ±3 we see that now for £z = ±2 
the current density vanishes somewhat slower (~ -i?^) for approaching the 2;-axis (-i? — t- 0). 
Clearly, this rule is continued also to the case £z = il: 

4 r 56 1 

^^^k^{^)\^^^^^ = ±— sin^ ^y-Y si^^ ^ + 12 sin^ ^ > , (IV.47) 

which meets with the expectation that for £z = ±1 the current density vanishes now 
slowest (i.e. ~ t?^) for approaching the 2;-axis (?? — )■ 0), see Fig. IV. A. 
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^ sin(7a) 
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71 sin(d) 



Fig. IV.A.O Azimuthal Current Density ^P^fc^(i9) ^IV.44\) 

For = 0, the current density ^^^k^ is non-zero on tlie z-axis (i!) = 0, tt) wliicli signals 
the existence of a magnetic singularity just on the axis of rotation. For the ground- 
state {Ip = 0), there occurs an additional singularity of the electric type (see Fig.4b in 
ref [Ij). This fact might eventually be the origin of the relatively large deviation of the 
RST prediction for the groundstate energy (-7,2305 [eV])) from the conventional predic- 
tion of —6, 8029 . . . [eV] see equation (]I.2p for n = 1. Thus the expectation is that the RST 
groundstate prediction will be shifted towards its conventional counterpart by a more ad- 
equate treatment of both groundstate singularities. Observe also the enlargement of the 
current density close to the axis of rotation (z-axis) which keeps the angular momentum 
as small as possible (^ 1^ = 0). 
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Fig. IV.A.l Azimuthal Current Density ^P^k^{'d) (TV^J^ 



For Iz = ±1, the current density ^'^^k^{'d) takes its maximum aside of the axis of 
rotation (in contrast to 1^ = 0, precedent figure). This may be intuitively associated 
with the larger amount of angular momentum (i.e. Iz = ±1 in place of 1^ = 0). Since 
the current density vanishes now on the z-axis (?? = 0, vr) one expects that no magnetic 
singularity can occur here. 
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Fig. IV.A.2 Azimuthal Current Density ^P^k(p{'d) ^IV.46\) 



For Iz = ±2, the maximum of tlie current density ^'^^k(j,('d) becomes shifted farther 
away from the axis of rotation {-t} = 0, tt) which may intuitively considered again to be 
due to the even larger amount of angular momentum (i.e. \lz\ = 2 in place of = 1, 
precedent figure). The peculiarity of this specific angular distribution of charge is that the 
associated quadrupole correction (II V. 104c]) of the spherically symmetric potential tPly4o(r) 
is zero; see the discussion below equation (lIV.lOSp . The small region around d = ti /2 has 
opposite direction of rotation. 
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Fig. IV.A.3 Azimuthal Current Density {P^fc0(i?) ^IV.45\ ) 



For = ±3, the maximum of tlie current density ^'^^k(f){i!}) is now farthest away from 
the axis of rotation which meets with the intuitive expectation that such a configuration 
should possess the largest value of angular momentum due to a given value of /p 
(i.e. Iz = ±/p)- Simultaneously, this specific geometry of the flow pattern generates 
also the largest amount of electrostatic(!) anisotropy energy ^111 ^ flIV.128j) within an 1^- 
multiplet, see the corresponding estimate flIV.129p . Thus, this effect of dislocation of 
charge and current to the outside regions with increasing angular momentum 1^ turns 
out to be the origin of the elimination of the /^-degeneracy occurring in the spherically 
symmetric approximation. The whole energy spectrum due to the present subset defined 
by Iz = ±/p is studied in Subsect.IV.4. 
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3. Energy Functional and Principle of Minimal Energy 



In order to decide whether the states due to different values of ^z (but the same 
value of l-p) do carry the same energy, one needs an energy functional (-E{t}; say)- This 
object will then be equipped with two intimately related meanings: First, its value on 
any solution of the positronium eigenvalue problem (precedent subsections) yields the 
binding energy of the corresponding positronium state; and second, it provides us with 
the basis of the variational principle of minimal energy. This fact may then be exploited 
in order to construct approximate variational solutions so that the values of the energy 
functional -E{t} on these solutions gives to us a first rough estimate of the energy spectrum 
of para-positronium. 

Of course, in view of its importance the construction of such an energy functional -E{t} 
has been the subject of several precedent investigations [Zitate] so that we can restrict 
ourselves here to merely quoting the relevant results. In principle, the desired functional 
is composed of two different kinds of contributions, i. e. the physical part Ej^r^j and the 
constraint terms which are neccessary in order to reproduce the right extremal equations. 
Of course, the latter must then be identical with the system of eigenvalue equations. 
Thus, the choice looks as follows [1] 

^[T] = e[JJ) + 2\d ■ iVo + Ag) ■ iV^^) (IV.48) 

with the physical part being given by 

Ef^l^ = 2Mc^ -Zl + A (P>rkin + . (IV.49) 

No doubt, the more interesting part is the latter one, because for the practical purpose 
of extremalizing the total functional (IIV.48P one will select those trial configurations 
which per se annihilate the constraints N-q and N^^ (see below) so that one can deal 
exclusively with the physical part flIV.49|) . which is the sum of renormalized rest mass 



energy, kinetic energy and gauge field energy. Naturally, all these different contributions 
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are to be understood as functionals of the basis RST fields (i. e. matter field \l/ and bundle 
connection A); but in the course of our investigations we did reparametrize these fields 
ultimately in terms of the angular functions fR{^),fs{^) (HV.6ap -( 1IV.6bp . spherically 
symmetric amplitude fields $±(r) (]IV.7ap -( ]IV.7bp . and electric part ^^^Ao^r) of the gauge 
field where we first may resort to its spherically symmetric approximation ^PAo{r) {r = \ f\). 
Consequently, we first have to specify all the energy contributions in terms of these latter 
fields before we can determine the energy spectrum of para-positronium, namely through 
deriving great benefit from that energy functional ( ]IV.48p . 

Relativistic Mass Renormalization 

The first term of the physical part ^^^'^ (]IV.49P is not simply the rest mass energy 
(2Mc^) of both particles (i. e. electron plus positron) but contains also the renormalization 
factor which is defined by 



= y d-^r i)'p{r)i)'p{r) , (IV. 50) 

or in terms of the Pauli spinors "-^ V± (0 iii place of the Dirac spinors ilj-p (r*) flIII.49aP 



0-2 



Inserting here the former decomposition flIII.37ap - flIII.37bp of the Pauli spinors with re- 



spect to the bosonic spinor basis yields the renormalization factor in terms of the ampli- 
tude fields (P>;^±,(P)5± as 

^ J (p)7^^ . (pyR+ + ■ (IV.52) 
_(p)^_ . {pyji_ _ {v)*s_ . {p)s_ 

and since furthermore the original amplitude fields '■^•*7^± , *^^^iS± have been ultimately trans- 
formed to the angular functions f r{'&) , f si'^) and radial amplitudes $±(r) via the equa- 
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tions firvriill -f irObll . firTOl) -f Hv:6b|l . and (IIV2al)-(|IVlb]) one ends up with the follow- 
ing form of the relativistic normalization factor 

oo 



2-2 



0-2 



47r sin d 



$^ (r) 



(IV.53) 



r=0 



Thus the final form of the (renormalized) rest mass energy of the two particles is found 
to appear as 



2Mc2 • 



2Mc^ 



dr r 



Kir) 



$2 



(IV.54) 



provided one agrees to apply here the separate angular normalization condition (]IV.15aP . 

Since (for the present prelimary estimate of the positronium spectrum) we will be sat- 
isfied with the non-relativistic approximation, we have to reveal now the non-relativistic 
limit of the renormalization factor Zj, flIV.54p . For this purpose, one might perhaps be 
tempted to simply neglect completely the "negative" constituent <l>_(r) of the amplitude 
field and consider the "positive" constituent $+(r) as the true non-relativistic approxi- 
mation. However, such an argument would be somewhat too superficial. Actually, we 
have to observe here that for the relativistic treatment the "negative" constituent ^-(r) 
is needed in order to obey the radial normalization condition f lIV.15b|) . Even if we take 



the viewpoint that 



$_fr) 



, it cannot be neglected 



is always much smaller than 
for the transition to the non-relativistic limit! Thus we have to keep in view $_(r) at least 
approximately; and we do this by solving (approximately) the first eigenvalue equation 
(lIV.Qap for $_(r) in the following form 

h 



$_(r) 



2Mc 



^ 1 

dr r 



(IV.55) 



where the mass eigenvalue is (approximately) identified with the rest mass M of the 
particles. Next, substituting this result for $_(r) in the radial normalization condition 
f lIV.15bp yields 



dr r $i 



— ) 

2Mc/ 



r=0 




(IV.56) 
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and thus the radial part of the renormahzation factor Z-p (IIV.SSP becomes 

oo 



dr r 








2Mc / \\ dr 



(IV.57) 

So we see that the non-relativistic approximation of the mass renormahzation term flIV.53p 
is not buih up only of the particle rest mass energy (Mc^) but also of the non-relativistic 
kinetic energy "^P^i^kin, i. e. we find 

Mc^ ^ - / (/«(^) + ■ ^"^^kin[$] = Mc' - (^')Ekin[$] (IV.58a) 

. (IV.58b) 

This result is interesting in so far as it says that the rest mass contribution (first 




term) to the physical energy -E^^T^ (IIV.49P actually is the proper rest mass energy 



IT] 



minus 



the non-relativistic kinetic energy (P^Ekin- This is the reason why the physical energy 
-Ep,^^ (IIV.49P (for two particles!) contains the four- fold kinetic energy '-^■^Tkin of a single 
particle, namely in order that the two-fold negative kinetic energy in the rest mass term 
2Mc^ ■ 2p f lIV.58aP can be compensated for. Clearly, for the non-relativistic form E[t] of 
the relativistic -E[t] f lIV-48p we can omit that ordinary rest mass term (2Mc^) because it 
does not take part in the dynamics at all. 

Kinetic Energy 



The second physical contribution to the energy functional i?[T] (IIV.48p - (IIV.49p is four 



times the one-particle kinetic energy ^^Tkin- Analogously to the precedent case of the rest 
mass energy, this contribution must now also be expressed in terms of the fields ^± and 
/i{('i9), fs{i^)- To this end, it is very convenient and also instructive to subdivide the total 
kinetic energy '■^Tkin into three parts which are associated with the motion along the three 
spatial directions being specified by the spherical polar coordinates r, {}, cj): 

ipyr^.^ = ip% + (p>r^ + (p>r<^ . (iv.59) 
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Turning here first to the radial part ^^^r one is easily convinced that this energy contri- 
bution is itself the sum of two terms 



(p>r, = (p>^,[7^] + ^""Tris] , (iv.eo) 

where the first one is built up solely by the amplitudes ^^^R±{r,'d) (lIV.4aP 



~ he r 8 ~ 

(p>^,[7^] (pyTr[n]=-^ / rf¥{(f)i?_(r,^9) ■ ^(^')i?+(r,^?) (iv.ei) 

or r ^ 

{(ff= r dr (I'd) 

and the second term ^''^^[iS] looks identical apart from the replacement ^P^R±{r,'d) =^ 
(P)^±(r,t9). Thus we are merely left here with the task to apply the last transformations 
(IIV.6a|) -( 1IV.7b|) which then ultimately brings the total kinetic energy ^'^Tr of the radial 



type (1IV.60|) to the following form 



71" OO _ 

(^'>f, = Y / + fsi^)] I dr r|<l_(r) ■ (IV.62) 





dr r 

Properly speaking, the product ansatz flIV.6ap - flIV.7bl) represents a kind of pre- 
selection among all the possible matter field configurations so that the radial kinetic 
energy f]IV.62p becomes a functional of only four very restricted functions of r and 
i.e. $±(r) and fsi''^)- But it is also possible to go one step further by demanding 

that the angular functions obey the normalization condition fllV.lSap which may also be 
written as 

TV 

Nf = \jdd {fm + fm) -1 = 0. (IV.63) 
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In this case, the radial kinetic energy ^"^^r becomes reduced to a function of $±(r), i. e. 



{pyf 



he 



, , . , , rf$+(r) -- ci$_(r) 



dr 



dr 



(IV.64) 



so that the extremahzation procedure (due to the principle of minimal energy) can refer 
only to the amplitudes $±(r). On the other hand, one may prefer to subject the energy 
functional also to the extremahzation with respect to the angular functions fni'^), fs{^) 
which then has the consequence that we have to work with the original functional ^^Tr 
(1IV.62P but have to add in this case the angular normalization condition (]IV.63p as a con- 
straint. Observe also that the radial normalization condition ( 1IV.15b|) has to be satisfied 
in any case. 

Next, the energy ^^T^ due to the longitudinal motion is found to appear in terms of 
the amplitude fields ^^^R±, ^^S± in the following form: 



he f d'^r 



d 



d 



^P^R^(r, m ■ TT77^^^'5+(r, ^) - ^^^S+ir, ^) ■ —^P^R^(r, ^) (IV.65) 
ov ov 

ov ov 



If the former product ansatz 0IV.6ap - flIV.7bl) is again introduced herein, one finally ends 
up with the longitudinal analogue ^^T^^ of the radial energy ^^^T^. ( ]IV.62p . i.e. 



he 



dd 



dd 



d§ 



dr $ I fr) • (^_{r) . 



(IV.66) 



Finally, the azimuthal part '^^'^tf, of the kinetic energy ^^^kin f irvMj) is found as 



he 



d'^f 



(pyt^ =—i^l—^ (p)5^(r, ^) . (P)5_(r, ^) - ^P^R+{r, ^) ■ (r, ^) (IV.67) 



+ cot?9 
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or by applying again those transformations (lIV.6a|) - (]IV.7b[) 



OO 



= ^4 / M { fl{d) - fl{d) +2coidfR{d)fs{d) \ I dr |.+ (r) • <l_(r) . (IV.68) 







Electric Interaction Energy 

According to the principle of minimal coupling, the interaction of both particles occurs 
via the gauge potential Afj_, cf. flII.3p - flII.5P : but since we are dealing with non-identical 
particles and do also neglect the magnetic forces we have to observe merely the electric 
part ^^^Ao(f) of the non-Abelian bundle connection A^. Moreover, recall that we are 
presently satisfied with the spherically symmetric approximation where the interaction 
potential (^Uq (r*) is assumed to depend only upon the radial variable r of the spherical 
polar coordinates (i.e. we put '^^^y4o(r') =^ ^^AQ{r)). The interaction energy can now be 
expressed by the interaction potential ^^Aq^t) in two different ways. The first possibility 
consists in identifying the interaction energy with the gauge field energy {E^, say) due 
to the electro-static potential [2] 

he 



^ Ana 



s 



j d^f V(*'Uo(f) ^ (IV.69) 



which becomes simplified in our electrostatic approximation to 



OO 





The second possibility refers to the ''mass equivalent" of the gauge field energy E^'^ 
which is defined in the following general way 

M^^^c' = -he J d¥(PUo(f) ■ (P)/to(r) , (IV.71) 

i. e. in the special case of our spherically symmetric approximation 

OO 

M(")c^ Mt^lc^ = -he j dQ ^P^A;o('^) j dr r^ ^P%{r) ■ ^P^A;o(r) , (IV. 72) 
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where we have already made use of the factorization (11 V. 141) of the charge density ^^^ko(r). 
Here it is also reasonable to apply the angular normalization condition (llV.lSaP so that 
the mass equivalent Mt^lc^ (HVTD appears in its final form as 

oo 

MHc2 = J drr ^^%{r) ■ <l^(r) + ^i{r) . (IV. 73) 



By this procedure, the angular functions fs{^) disappear from the mass equivalent 

quite similarly to the situation with the kinetic energy ^^Tr. Therefore we have to include 
the angular normalization condition flIV.63p as an extra constraint for the subsequent 
principle of minimal energy. But the crucial point with the mass equivalent Mt'^'c^ is now 
that it must be identical to the gauge field energy (]IV.70|) which may be reformulated 
as a further constraint for the variational principle: 

iVgl = E^^^ - M^c^ = . (IV. 74) 

Relativistic Principle of Minimal Energy 

Collecting now all the results, we find that our original energy functional -Ejt] (1IV.48P 
reappears in its final form (but restricted to the electrostatic and spherically symmetric 
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approximation) as follows: 



E 



[T] 



eL =2Mc^ ■ I drr 



<l>i(r) - $^ (r) 



+2hc I drr {^^{r) P-^-(^Ar)- 



(IV. 75) 





7r 

+2hc j d'd 





dr 



dr 

d^ 



dr ^+(r) ■ ^^(r) 



+24^c / d^ {fsi^) - fU^) + 2 cot ^ fR{^)fs{^)} / rfr<i>+(r) ■<!>_( 



dr 





oo 

/ (ir 





+ Ag^- iVgl + 2Xf ■ Nf + 2Xd ■ • 



In this result, one recognizes the original form (IIV.48p - (IIV.49p of the energy functional Ep^] 
in a more elaborate form: the last three terms are the constraints, where the normalization 
fHVl3|) of the density ^P^ko{r) in the constraint form 



Nd = j d^f^P^koir) -1 = 0, (IV.76) 

with the Langrangean multiplier 2Xd, is split up into the separate angular normalization 
constraint ( 1IV.63P and its radial analogue 



dr r 



<ll(r) + <l^(r) 



1 = 0. 



(IV.77) 



Of course, the latter constraint is nothing else than the radial normalization condition 
f llV.lSbp . The associated Langrangean multipliers are denoted by A/ and A^, resp. All 
other terms in the final form (1IV.75P of the energy functional E^^^^ do refer to its physical 
part Ef^l^ (HYani), i.e. the sum of the (renormalized) rest mass energy (lIV.54p . three 



kinds of kinetic energy (IIV.64p . (IIV.66p . fllV.GSp . and the gauge field energy (IIV.70p . 

Once the right energy functional is written down, one can look now for the extremal- 
izing configurations in order to check whether the system of eigenvalue equations is really 
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reproduced. First, turn to the extremal configurations of (1IV.75P witli respect to 
tlie angular functions fni^^) and fs{i^) and find just those first-order equations ( llV.Sap - 
f lIV.8b|) . provided the Langrangean multiplier A/ is identified with the quantum number 
ip in the following way 

^ ^ = . (IV.78) 



2hc J dr $+(r) ■ $_(r) 



Next, one wishes to look for the extremal equations for the amplitudes $±(r) and for 
the potential t^Uo(r). But for this purpose it is very convenient to first eliminate the 
angular functions fR{^),fs{^) from the energy functional -Ej^j flIV.75p . namely just by 
use of their extremal equations f lIV.8ap - (IIV.8bp . In this way, one is left with the reduced 



functional which then acts only over the configuration space of the amplitudes ^± (r) 
and the potential ^P^Ao{r): 



=2Mc' I dr r 



oo oo 

2 



$i (r) - (r) 



he 7, 2 fd^^'^ir]' ' 
dr r ' 




ocs J \ dr 



oo ^ _ _ 

+2hc [ drr I ^(r) ■ ^1±M. _ $^(^) . _ (i + 2£p) ■ 

J \ dr dr 



Observe here again, that the constraints appear in the third line; all other terms are 
of truly physical nature, i. e. the rest mass energy and the electrostatic interaction energy 
(first line) and the kinetic energy (second line). Together with the angular functions 
fR{^),fs{^) there has also disappeared the quantum number of angular momentum 
which says that the energy does not explicitly depend on iz', however, there is of course 
an indirect dependence on via the interaction potential ^^'^Aq^t). The reason for this is 
that the potential [^Uo(r) depends on the electrostatic charge distribution ^P^kQ{r,'d), see 
the Poisson equation (lIV.lOp . which takes account of the angular functions fR{'&), fs{'&) 
and therefore also of the quantum number i^. Consequently we have to expect a certain 
lifting of the energy degeneracy with respect to angular momentum ij 
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But concerning our reduced functional -Ej^j (]IV.79p . it is now a standard variational 



procedure to reproduce from it the system of eigenvalue equations: the extremaliza- 
tion with respect to the amplitude fields $±(r) generates the mass eigenvalue equations 
(IIV.9ap -( ]IV.9bp . with the Langrangean multiplier A^j being given by 

= -M.c^ . (IV.80) 

However, since we are satisfied for the moment with the spherically symmetric approx- 



imation, the extremalization of i?r$] (IIV.79|) with respect to the gauge potential t^^Ao 



does of course not yield the Poisson equation f llV.lOp but rather the spherically symmetric 
approximation thereof: 

Ar [^ko(r) = -as ■ +^ ^ ^ . (IV.81) 

r 

This may be understood as a kind of angular average (over the two-sphere) of the original 
Poisson equation (lIV.lOp . 

The spherically symmetric Poisson equation fllV.Sip closes now the present RST eigen- 



value problem in the spherically symmetric and electrostatic approximation. Its formal 
solution is given by 

bU,(,) = ^/:^!i:*M±|^ (,V.82) 

^ ^ Air J r' \\f-f'\\ ^ ' 

which may be used in order to calculate the potential [Ply4o(r) in terms of the amplitude 
fields $±(r). For the latter fields one may try some ansatz, as realistic as possible, and may 
then substitute the result of integration [pUo(r) in the energy functional (]IV.79p in 



order to determine the ansatz parameters in the trial amplitudes through extremalization 
of The result of this procedure would then be the desired energy spectrum (see 

below). However, a somewhat better picture of this spectrum will surely be attained if we 
succeed to go (at least partially) beyond the spherically symmetric approximation ^P^AQ{r). 
For this purpose, one may reconsider the original Poisson equation fllV.lOp whose formal 
solution is given by equation (lIV.lip . Clearly, one expects that if we will be able to use 
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that anisotropic potential '•^■'Ao(r*) (llV.llj) in place of its isotropic aproximation [^Uo(r) 
(1IV.82|) . then we will get an improved energy spectrum (see below). The improvement 
must especially help to clarify the existence of an energetic degeneracy of the states due 
to different quantum numbers but the same ip. 

Non-Relativistic Energy Functional 

Up to now, our line of arguments referred always to the relativistic situation; but the 
main purpose of the present paper aims at the deviations of the RST predictions from the 
conventional non-relativistic spectrum fll.2p (a comparison of the relativistic corrections 
between both theoretical frameworks must be deferred to a later study). Therefore one 
is now first interested in the non-relativistic approximation to both the energy functional 
-E[<i>] (1IV.79P and its extremal equations (i. e. the non-relativistic form of the eigenvalue 
equations (lIV.9ap - (]IV.9bp in combination with the Poisson equation (llV.lOp ). Of course, 
one will demand from reasons of internal consistency that the wanted non-relativistic 
field equations should represent just the extremal equations due to the non-relativistic 
form of the energy functional. Thus we will first look now for the non-relativistic energy 
functional (]E[$], say) by working out the non-relativistic form of its individual energy 
contributions. 

The non-relativistic form of the first contribution to the relativistic functional 
(lIV.79p . i. e. the renormalized rest mass energy, has already been determined by equation 
pV.58aP and ultimately appeared in the following form: 



2Mc^Z^ 2Mc2 -2 (%kin[<i.] , (IV.83) 



namely by making use of the angular normalization condition fllV.lSap . Here, the non- 



relativistic kinetic one-particle energy *^^^ii^kin is specified by equation fllV.SSbp for which 
the "negative" Pauli amplitude (r) has simply been neglected. Indeed, this amplitude 
can be eliminated completely for passing over to the non-relativistic approximation so 
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that we can deal here exclusively with the "positive" Pauli component $+(r), being 
henceforth simply rewritten as <i>(r). Thus the non-relativistic one-particle energy '•^■'-Ekin 
appears simply as 

oo ( / ~ \ 2 








The second energy contribution to the functional i?[<j,] 0IV.79P (first line) is the elec- 
trostatic gauge field energy E^^ ( lIV.TOp . which is built up exclusively by the electrostatic 



potential [^Uo('")- Therefore this term does not at all change its formal appearance for 
the transcription to the non-relativistic situation. But, of course, the relativistic *^^Mo(r) 
cannot preserve its functional form for the non-relativistic approach. Indeed, it is easy to 
deduce this desired non-relativistic form from its relativistic predecessor flIV.82p : namely 
by simply omitting again the "negative" Pauli component <l>_(r) ^<l>+(r) ^ <l>(r)^ 

.U„„.^/^ *!(!;L. (,v.85) 

47r / r r — r' 



see the corresponding remarks above equation flIV.84p . However, one must not forget here 



the fact that both potentials [PUo(r) f lIV.82p and f lIV.85p are a consequence of the spher- 
ically symmetric assumption ((*'Mo(r) =^ [^ly4o(r)). Subsequently we will go one step 
further beyond this spherically symmetric approximation, namely be expanding the orig- 
inal anisotropic potential '^'^^Aoif) f llV.lip around the isotropic potential [^Uo(r) (lIV.85p . 



But for the moment we are satisfied with merely deducing the non-relativistic approxi- 
mation of the relativistic functional (1IV.79P within the framework of the spherically 
symmetric approximation. 

The last contribution of physical nature to the energy E^^ flIV.79P is the kinetic energy 
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^^^kin[4>] (second line) 

^^^kinf*i = 2hc f drr h(r)- ^1±M. _ $^(^) . (iv.86) 
J \ dr dr 



The non-relativistic version hereof can easily be obtained by simply solving (approxi- 
mately) the first eigenvalue equation ( lIV.QaP with respect to $_(r) as 

^-ir)^4rl'^---^^ir)\ (IV.87) 
^ ^ 2Mc ] dr r +^ M ^ ^ 



and then substituting this into the second eigenvalue equation fllV.Qbp which yields the 
four-fold non-relativistic kinetic energy ^"^^ Ei^i^^^^^ flIV.58bp : 

^^^kin[$] 4(P)Ekin[$] (IV.88) 

oo ( / ~ \ 2 

2^2 /• (d^{ry 



, dr r . , 
M J \ \ dr 





V 




Finally, it remains to consider also the non-relativistic forms of the constraints for the 
principle of minimal energy (see last line of equation ( ]IV.79|) ). The relativistic Poisson 



constraint n]^ (1IV.74P is nothing else than the difference of the electrostatic field energy 
f lIV.7Qp and its mass equivalent Mt^'c^ fHV73|l . The non-relativistic version (M'^'^c , 



say) of the mass equivalent is obtained again by simply omitting the "negative" Pauli 
amplitude $_(r) and putting $+(r) ^ <l>(r): 

oo 

MI^Ic^ M''^c^ = -fm j rMAo(r)l>2(r) , (IV.89) 



so that the non-relativistic Poisson constraint (Nq ) reads 

iVgl ^ = E^^ - M'^c" = . (IV.90) 
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And finally, the non-relativistic normalization term (N$) due to its relativistic predecessor 
N<s> (1IV.77|) must evidently look as follows 



2 



=^ N$ = dr r^\r)- 1 = . {W.91] 







Thus, collecting all the non-relativistic contributions, one ultimately finds for the non- 
relativistic approximation E[$] of the relativistic case (]IV.79|) : 

E[<,] ^ = 2 (^)Ekin[$] + + Ag^ ■ Ng + 2 A, • N$ (IV.92) 

(here the Langrangean multiplier A^ is to be understood as the non-relativistic version of 
Xd). Observe also again that the proper rest mass energy 2Mc^ is omitted (because of 
its dynamical ineffectiveness) and that the relativistic renormalized rest energy 2 MZj, 
combines with the four-fold relativistic kinetic energy (P^fkin, cf. f irTOj) . just to the 
two-fold non-relativistic single-particle energy (^'^-Ekinf*], as it must be expected for a two- 
particle system. Clearly, this supports the intrinsic consistency of the relativistic and 
non-relativistic versions of our energy functional. 

Spherically Symmetric Eigenvalue Equations 

A further consistency test refers to the extremal equations of the non-relativistic func- 
tional E[$] (lIV.92p . Indeed, one expects here that those extremal equations will coin- 
cide with the non-relativistic approximation of the relativistic eigenvalue system f llV.Dap - 



(lIV.Qbp to be complemented by the non-relativistic version of the relativistic Poisson 



equation f llV.Sip . i.e. 



r 1 $^(r) 
AWAo(r) = -a,^^, (IV.93) 
r 



which has the formal solution displayed by equation (lIV.SSp . But indeed, it is easy to see 
that this consistency test is positive: the extremalization of the present non-relativistic 
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functional E[$] (1IV.92P with respect to the amphtude field $(r) generates the following 
Schrodinger-like equation 

1 \ ^ _ h'' 

2M \dr^ 

and furthermore the extremalization with respect to the gauge potential [^Uo(r) lets 
emerge just the spherically symmetric Poisson equation (lIV.93p . On the other hand, 
we can deduce also this Schrodinger-equation flIV.941) as well directly from the rela- 
tivistic eigenvalue system (IIV.9ap - (IIV.9bp : for this purpose, resort again to the former 



equation flIV.87p and substitute this into the second equation pV.9bp in order to find 



2M V dr2 



= (M, - M) c^ ■ l>(r) . 



Comparing this result to the former one (1IV.94P yields the identification of the non- 
relativistic Lagrangean multiplier with the difference of mass energies 

A, = (M - M,) c^ = -E, , (IV.96) 

i. e. the non-relativistic multiplier As plays the part of the Schrodinger eigenvalue £"* in 
the Schrodinger-like equation flIV.95p . 

The present spherically symmetric approach to the positronium eigenvalue problem, 
consisting of the equations f lIV.92p - flIV.96p has been used in order to get a first rough esti- 
mate of the energy spectrum |3]; and it has been found that the corresponding predictions 
deviate from the conventional predictions fll.2p by some 5-10%. This result is not too 
bad, especially when one recalls that it came about by neglecting the anisotropy of the 
interaction potential ^^^Aq^t). But in the strict sense, this potential cannot be spherically 
symmetric, not even for the product ansatz described by the equations f lIV.4ap -f lIV.4bp . 
f lIV.6ap -f lIV.6bp and f lIV.7ap -f lIV.7bp . see the integral representation fllV.lip in terms of 
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the charge density ^P^ko(r). This anisotropy of ^^^Ao(r) does exist also for all the excited 
states even if they have vanishing ^-component of angular momentum (-w = 0, > 0). 
Therefore, if we wish to get a more accurate RST prediction of the positronium spectrum 
(especially for ^ 0), we have to regard also the anisotropy of that interaction potential 
(pUo(r), at least in some approximative sense. 

Anisotropic Interaction Potential 

The point of departure for studying the anisotropy in question is, of course, the integral 
representation fllV.lip of ^^^AQ{r), which in the non-relativistic approximation becomes 
simplified to 



OO 



dr' r' t—^-4tj ■ (IV.97) 



4:71 sint?' J \\r — f 



For the subsequent expansion of this anisotropic potential with respect to the magni- 
tude of anisotropy it is very helpful to first clarify its relation to the spherically symmetric 
approximation ^p^Aq^t) (jIV.SSp . To this end, consider the average of the anisotropic po- 
tential ^P^Aoir,^) over the 2-sphere; indeed, for this one finds by explicitly integrating 
over the solid angle fl: 

4tt J Att smv' 

where one has made use also of the angular normalization condition (llV.lSap . Thus, 
restating this result because of its importance for the subsequent perturbation expansion, 
we note down for our spherically symmetric approximation 

^P%{r)=[ — ^PUo{r,^). (IV.99) 
J 47r 

This namely suggests to define the anisotropic part A^^{r,'d) of ■>^^Mo(r, -(9) as just its 
deviation from the spherically symmetric approximation, i. e. we put pQ 

{r,^) = [^ko (r ) + (r, t9) . (IV. 100) 
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Or conversely, combining both integral representations (1IV.85|) and (]IV.97p one gets the 



following integral representation for the anisotropic part: 

A"-(r,79) = ^P%{r,^) - [PUo(r) (IV.lOl^ 



as 



47r 



sin ■}}' 



dr r — — — 







Would we neglect here the angular dependence of the radial (i. e. the last) integral, 
then we would find the anisotropic part vanishing (A^'^(r, -(9) =^ 0) and thus the potential 
^^Mo(r, ^9) (IIV.97P would coincide with its spherically symmetric approximation [^lylo(r). 
Therefore the angular dependence of ^'P^AQ{r, t?) flIV.97p comes about through the regard of 
the angular dependence of the radial integral! This fact may now be exploited in order to 
establish a perturbation expansion of the anisotropic part A^'^lr,-}}) f llV.lOip with respect 
to the magnitude of the anisotropy, namely just by constructing such an expansion for 
that radial integral. 

The point of departure for the wanted perturbation series is the following expansion 
of the denominator in the radial part of (lIV.lOip 

- + 3ir'r')+- F [r-r') (IV.102) 



2 



5 r^-r'^ .n 3 35 r'^ . r'^ 
+ o [^'^ I + s 9" [''""'" I + 



2 ^^2 _|_ ^/2^ 2 V / 8 ^^2 _j_ ^/2^ 2 



Substituting this series expansion into the integral representation fllV.lOip for the anisotropic 
part A^'^lr, ^) yields an expansion of the following kind: 

A'^-(r,^9) = ■ ^^'U'(r) + • ^PU"(r) + ^^^"'(^9) ■ ^^'U"'(r) (IV.103) 

Here, the angular pre-factors in front of the radial auxiliary potentials ^P^A\r), ^^^A^\r), 
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^^^yl"'(r) depend on the angular functions fsii^) in the following way 



i.e. for = 3: 



dQ' 

4:71 

3 r dn' 

2./ 17 



1 



sin '(9' 

sin d' 
sin -(9' 



= 
1 



(IV. 104a) 
(IV.104b) 
(IV. 104c) 



i[cos2^-|] ,4 = 
|[cos2^-|] , 4 = ±1 
, 4 = ±2 

[cos^ ^9 - |] , 4 = ±3 



15 
32 



(the result for = does hold for any /p, see App.B). 

As a check of the results for these angular pre-factors, one deduces from the combina- 
tion of equations f lIV.99|) and (]IV.100|) that their angular average must vanish, i. e. 



^^^]^pU"'(^?) = dn ^pU^(^?) 



. 



(IV.105) 



This is trivially satisfied for ■^^^'"('i?) (]IV.104cP with 4 = ±2; indeed the vanishing of 
{p}^iii(^^~) for 4 = =t2 says that in this case the spherically symmetric approximation 
[ply4o(r) is a very good approximation and corrections can arise only in the higher pertur- 
bation orders. This circumstance equips the spherically symmetric approximation with 
additional relevance, since its predictions remain partially true even in the first order of 
perturbation. Observe also that the first non-trivial order of anisotropic perturbation is 
of the third kind flIV.104cjl : the first two correction terms ^^^'(^9) flIV.104aj) and ^pU"(^9) 
(lIV.104bp must always vanish. 

For the radial correction factors ■>^^M'(r), ^'^^A^\r)^ ... of the anisotropic part yl™(r, -t?) 
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(1IV.103P one obtains the following results [T] 



{p}^iii( 



dr' 



■ r I dr' 



(r2 + 



5" 

^2 _j_ j^/2j 2 



(r2 + 



(r2 + r'^) 



(IV.106a) 
(IV. 106b) 
(IV.106c) 

(iv.ioed) 

(IV.106e) 



After having selected some trial amplitude $(r), one could in principle explicitly cal- 
culate these radial potentials in terms of the ansatz parameters contained in the selected 
trial amplitude. However, our method of approximation just consists in leaving unspeci- 
fied these (non-trivial) radial correction potentials ■>^^^y4"'(r), ^^^A^ {r) etc. and determining 
them through solving their extremal equations due to the principle of minimal energy. 
Thus, the next step must concern the reformulation of the energy functional E[<j,] flIV.92p 
in terms of the chosen set of field variables $(r), ^^'^AQ{r)^ ^^^A^^\r). For a first rough 
estimate of the para-spectrum we restrict ourselves to the first non-trivial anisotropic cor- 
rection (~ •fpU"'(r)), i. e. we let the expansion fllV.inai) of ^PU^°(r,^9) stop at the third 
term: 



III/ 



(IV.107) 
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Energy of the Anisotropic Gauge Fields 



In order to get a first impression of the anisotropic effects we still retain the spherical 
symmetry for the non-relativistic amplitude field $(r) but we regard now the anisotropic 
effect on the electrostatic gauge field energy E^^ (^=^ ^r^^ its mass equivalent 
M^^'^c^ (^^ M^c^y The reason is that we think those effects to be the leading ones 
concerning the expansion of all fields with respect to the magnitude of anisotropy. 

The starting point for the calculation of the "anisotropic" gauge field energy E^^ is 
(IIV.69p . Inserting therein the present decomposition (lIV.lOOp of the anisotropic potential 
^P^Ao{r, with the anisotropic part A^^{r, t?) fllV.lOSp yields an analogous decomposition 
of the gauge field energy E^^ flIV.69|) into a "spherically symmetric" part e]^' flIV.70|) 
and an "anisotropic" part Ein , i.e. 

= J d'r 1 1 V WAo(r, ^) 1 1' = 41 + Ei^ , (IV.108) 

where the "anisotropic" part Ei^ is given in terms of the radial correction potential 
frU"'(r) by 

{00 2 °° 

(^^63 ■ jdrr^ (^^^^) + ^'^3 ■ / dr (i^U'"(r))^ 


(IV. 109) 

since we are satisfied for the moment with the third-order approximation. Here, the 
coefficients ^^•'63 and '^^•'/3 are determined by the angular correction potential ^^^'"(-(9) 
through 

We3 = j 1^ (Maiii(^9))' (IV.llOa) 



_ 47r V d^ 

Since the angular correction potential ^P^A^^\'d) is explicitly known, cf. equation flIV.104cp 
both coefficients are also known numerically, see the following table for i-p = 3: 
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4 





±1 


±2 


±3 




1 

80 


9 

1280 


n 


5 
256 


(p)f„ 


3 

40 


27 
640 


n 


15 
128 




1 

16 


9 
256 





25 
256 


(p)e3 


6 


6 




6 


(p)m3 
(p)e3 


5 


5 




5 



Besides the gauge field energy E^^ there is a further gauge field contribution to the en- 
ergy functional E[$] ( 1IV.92|) whose presence is needed in order that the extremal equations 
do coincide with the eigenvalue system, i. e. the Poisson constraint Nq ^ 



I^^>c2 







(IV.lll) 



which is the anisotropic generalization of the spherically symmetric constraint Nq as 
shown by equation OIV.QOp . The "anisotropic" gauge field energy E^^ is given by equa- 
tions (lIV.108p - (lIV.109p . and the anisotropic generalization M^'^j'c^ of the isotropic M^'^^c^ 
(1IV.89P must look as follows: 



01^2 



he / d^f^P%{r,^)^Pho{r,^) . 



(IV.112) 



Substituting herein again the decomposition of the anisotropic potential ^'P^Ao{r, fllV.lOOp 
and the factorization flIV.14p of the charge density ^'P^kQ(r,i!}) lets this mass equiva- 
lent M^'^^c^ fllV.112p reappear also as a sum of the "isotropic" part Mt'^'c^ flIV.89p and 



"anisotropic" part 



I.e. 



I{e}c2 



ie],2 



(IV.113) 
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where the "anisotropic" part consists here solely of the third-order correction potential 
WA"'(r,79) (= ■ WA"'(r)), i.e. 

oo 

Mf^c^ Mfif c^ = -hc (fWg ■ J drr ^^^'"(r) ^^(r) . (IV.114) 







Here, the angular part '^^^"'('(9) of ^P^A^^{r, -(9) f lIV.107p builds up the numerical coefficient 
^^Vns in the following way 



'ms = [dn ■ ^P^koii}) = [dn ■ ^^^f + ^'}^^ (IV.115) 

/ / 47r sm V 



ip), 

(for the numerical values of ^^Vna see the above table on p. 1841) . 

With respect to the principle of minimal energy there is now an important circum- 
stance; namely the splitting of the "anisotropic" Poisson constraint fllV.llip into 
certain subconstraints. First observe here that, due to the splitting of E^^ (IIV.IOSI) and 
its mass equivalent M'^^'^^'c^ (lIV.llSp . the original Poisson constraint splits naturally 
up into the "spherically symmetric" part Nq (IIV.QOP and an "anisotropic" part Nln^: 

= + , (IV. 116) 

with the "anisotropic" part being given in an self-evident way as 

Nl;^ = ^in^ - M^c^ , (IV.117) 

i. e. more concretely since we are presently satisfied with the first non-trivial approxima- 
tion order (~ 

Nf,^^ = E,^^ - Mf^c^ ^ . (IV.118) 

This means that, because of considering both gauge field constituents tp]Ao(r) and ^^^'"(r) 
as independent field degrees of freedom, we can impose both Poisson constraints fllV.QOp 
and fllV.llSI) independently from each other. Thus, the one constraint Nq^ = (lIV.llip 



is broken up into the two independent constraints Nq = f llV.QOp and N^^I = OIV.llSp 



85 



Summarizing the results for the "anisotropic" gauge field energy (i. e. more precisely: 
the energy of the anisotropic gauge field configurations), we find it convenient to split up 
the present energy functional ]E[$] (]IV.92p into a (Dirac) matter part '•'^'*lE[<i>] and a gauge 
field part *-'^'*IE[$] which, however, enters the present electrostatic approximation of the 
energy functional only via its electric component '■^^Ej^j, i.e. we put 

E[$] = (^>E[$] + (^>E[$] (IV. 119) 

with the matter part being given by 

(^>E[$] = 2 (p)Ekin [$] + 2 A, ■ N$ (IV.120) 

and the electrostatic gauge field part by 

= + Ag) ■ Ng^ . (IV.121) 

Obviously, each of both contributions consists of a physical term and a constraint; but 



the proper meaning of the splitting (1IV.116P refers to the principle of minimal energy 



S E[^] = (IV.122) 

which is required to reproduce the complete system of eigenvalue equations (including 
the gauge field part thereof). In this respect, the observation is now that the gauge 
potential ^P^AQ{r,'d) is contained solely in the electric part *^'^^E[$] flIV.12ip . Consequently, 
it is sufficient to consider only the latter part of the energy functional for the purpose 
of deducing the gauge field equations. More concretely, since the gauge field has been 
decomposed into the two independent field degrees of freedom WAo(r) and ^PU"'(r), the 



corresponding field equations must be deducible from flIV.1221) through the two separate 
extremalization processes 

<5o ^%'^] = (IV. 123a) 

Su\ ^%^] = , (IV.123b) 
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where ((^iii) denotes the extremahzation with respect to ^^^AQ{r) (''^^^yl"'(r)) , resp., and 
the matter part (1IV.120|) plays no part at all. 



Non-Relativistic Gauge Field Equations 

Now that all the details of the gauge field part 0IV.12ip of the energy functional 

have been clarified, one can proceed in determining the corresponding extremal equations. 
Since the spherically symmetric approximation [^ly4o(r) is conceived to be an independent 
field, the associated extremahzation process ( 1IV.123a]) appears in a more detailed form as 

^0 + =0 (IV. 124) 

and yields of course the former "spherically symmetric" Poisson equation (]IV.93p with 
Xq^ = —2. Clearly, this is not a new result; but the "anisotropic" counterpart f lIV.123bj) 



is just what we need in order to discuss the elimination of the energetic degeneracy of the 
states due to different (but the same i-p). 

That second extremahzation process flIV.123bl) reads in a more detailed form 

5„, {ES;^ + Ag)-Nf,?}=0 (IV.125) 



where the "anisotropic" energy E^^^ is specified by equation (]IV.109|) and its mass equiv- 
alent Mlfi'f-'^c^ by equation (lIV.114p . Thus the extremahzation process (]IV.125p generates 



the following quadrupole equation for {piyl"'(r) [Zitat] 



I (ir^ r dr J 

r 

Amazingly enough, this equation adopts the same form for all values of angular momentum 
iz (besides the trivial case = ±2 for £p = 3): 

c/2{pUiii(r) 2d^PU^^Ur) MA"'(r) ^\r) 

1-5 \ 1 6 = -5as , (IV.127) 

ctr^ r dr r 



see the table of numerical coefficients on p. [HH By comparing this quadrupole equation for 
the radial correction potential ^^^'"(r) to the Poisson equation (]IV.93p for the spherically 
symmetric approximation ^^AQ{r) shows that the coupling of ^^^'"(r) to the spherically 
symmetric amplitude field $(r) is five times stronger than the corresponding coupling of 
[^ly4o(r). This, however, does not entail that the energy E}^f fllV.inQj) located in the field 
(r) is greater than that being concentrated in the isotropic field ^^^Ao(r). Indeed, we 
will readily see that the anisotropy energy due to ^^^'"(r) is only a small fraction of the 
"isotropic" energy due to t^Uo(r). 

Though all the states with the same value of i-p (= 3) but different iz (= 0, ±1, ±3) 
share the same correction potential ■>^^M"'(r), this does not entail that they all possess the 
same anisotropy energy E^^^ (lIV.lOQp . In order to elucidate this fact, rewrite the latter 
quantity in the following form 

2 / r„i , III / \ \ 2 



^111 



so that these anisotropy energies for iz = 0, ±1, ±3 must obey the following interrela- 
tionship: 



^111 





^111 




^111 
























4=0 








lz=±3 



(IV.129) 



1 : 0,5625 : 1,5625 



(provided the source term ^~ <l>^(r)j in the quadrupole equation (IIV.127P is kept fixed). 
Furthermore, since the anisotropy energy Ej^^^ OIV.lOQp is always negative and the spheri- 
cally symmetric approximation alone predicts somewhat too high energies, the corrective 
property of the present approximation method is optimal for iz = =t3. Since in this case 
the value of iz is itself extremal (i. e. \iz\ = iv), such a state would correspond to circular 
motion in the Bohr-Sommerfeld picture of the atoms; and thus it would be interesting to 
calculate the whole spectrum (for i-p = \iz\ = 0, 1, 2, 3, 4, . . .) within the framework of the 
present approximation method (see below). 



Generalized Eigenvalue Equation for $(r) 

The final act of our principle of minimal energy must of course concern tlie eigenvalue 
equation for the non-relativistic amplitude field $(r). Obviously, the wanted equation 
will represent some generalization of the "spherically symmetric" eigenvalue equation 
(1IV.94|) since the present "anisotropic" energy-functional E[$] (]IV.119|) - (]IV.121|) differs 
from its "isotropic" counterpart (1IV.92P just by the inclusion of the anisotropy energy 
E^^^ (IIV.109p . Or in other words, the extremalization of the present "anisotropic" en- 
ergy functional ]E[$] (1IV.119P with respect to the amplitude field $(r) must reproduce 
that "spherically symmetric" eigenvalue equation flIV.95|) with additional inclusion of an 
"anisotropic" term (~ ■^^^^'"(r)). Indeed, the extremalization process with respect to 
$(r) yields now an eigenvalue equation of the expected form: 

-JL_ + ^r) + — ^ il ■ $(r) (IV.130) 

-he ([PUo(r) + ■ ^^U"'(r)) = E, ■ l>(r) 

where the non-relativistic energy eigenvalue E^, agrees again with the Lagrangean mult- 
plier Xg, see equation (IIV.96p . Here it is evident that the numerical coefficient '^^Wa 
(11 V. lisp plays the role of a coupling constant for the coupling of the amplitude field $(r) 
to the anisotropic part My4"'(r) of the gauge potential ■i^*'Mo(r, t?), see the perturbation 
expansion f llV.lOap of ^PU^'^'ir.^'}). Concerning now the coupling strengths of the ampli- 
tude field $(r) to the isotropic and anisotropic parts of the gauge potential, we see from 
the present eigenvalue equation (lIV.lSOp that the ratio of couplings 1 : ^P\n3 is in the 
order of magnitude 1 : 3^, or even smaller. This says that the coupling to the isotropic 
part [pUo( r) is dominant so that the anisotropy energy E}^f f lIV.128p can actually be 
expected to represent only a small correction to the result of the spherically symmetric 
approximation (see below). 
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Selecting a Plausible Trial Amplitude $(r) 

Being motivated by the acceptable success of the precedent results [I]- [2], our choice of 
an appropriate trial amplitude $(r) for extremalizing the energy functional E[$] is again 
the following: 

|.(r) = <^>,r'^e-'^^ (IV.131) 
where the normalization constant must adopt the value 

in order to satisfy the non-relativistic normalization condition fllV.Qip . Consequently, the 
constraint term N$ can be omitted for our energy functional E[$] (^M> or (^1$ fllV.120p . 
resp. Thus our trial ansatz (1IV.131|) contains two variational parameters (i. e. u and /3) 



for extremalizing the corresponding energy function E*^^^^(/3, u), see below. Surely, such a 
trial amplitude (lIV.lSip with only two variational parameters will be too rough in order 
to predict the RST energy levels with a satisfying accuracy; but it may be sufficient for 
the present purpose of demonstrating the level splitting for different values of iz (but due 
to the same value of i-p). 

As mentioned above, the suggestion for a first rough estimate of the para-spectrum 
consists now in taking the value of the energy functional E[$] (IIV.IIQP on the selected trial 
amplitude $(r) (IIV.13ip and to look for the extremal points of the corresponding energy 
function K^^^\(3, v). First, the value of the kinetic energy functional (P)^kin[#] fliV.84p on 
$(r) (lIV.lSip becomes by straightforward integration 

^kin(^^, P) = 7^ (2/3aB)' ■ £kin(^^) (IV.133) 
ZaB 

with the kinetic function £:km('^) being given by 
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Observe here that the kinetic energy does not depend on the quantum number ^z (but 
exclusively £p)! 

Next, consider the isotropic gauge field energy (1IV.70P which must equal its mass 
equivalent M'^^^c^ (lIV.89p . provided we can use for the calculation of these two objects 
the exact solution of the "isotropic" Poisson equation (IIV.93|) . In this case, we then 



can dispense with the "isotropic" Poisson constraint term in our energy functional E[<j,] 
fllV.119IWIlV.121l) . i.e. U§ = 0. However, such an exact solution (r) is easily obtain- 
able for our selected trial amplitude ^(r) (]lV.13ip - (]lV.132|) . Here it is convenient to pass 
over to dimensionless objects by putting 

tPUo(r) = 2/3«, -5,(2/) (IV.135a) 
'y = 2(3r) , (IV. 135b) 



which recasts that "isotropic" Poisson equation fllV.93|) in dimensionless form 



The exact solution hereof is given by [1], [1] 

~a,{y) = li-e-y.y — I . (1V.137) 

2z/+l| ^r(2z/ + 2 + n)^ J ^ ^ 

This exact solution can now be used in order to calculate explicitly both the "isotropic" 
gauge field energy ( I1V.70P and its mass equivalent M'^^c^ (]1V.89|) . and the result is 



the following: 



= M^c^ = (2 Pas) ■ £pot(^^) (IV.138) 

as 



with the potential function Epoti^) given by [T], |1] 



+ 1 y 24-+2 £^2" r (2z/ + 2) • r (2z/ + 2 + n) y ■ y ^- '^^) 
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From reasons of computational simplicity it is somewhat more convenient to determine 
the potential function epot(^) (1IV.138I) from the mass equivalent which yields 



-pot I 



r(2z/ + 2) 



dy y^''^^ My) \ (IV-MO) 



whereas the determination from the "isotropic" gauge field energy would necessitate 
to calculate a somewhat more complicated integral, i. e. 

^potM= ldyy^(^-^\ , (IV.141) 



dy 



which then would produce a double sum in place of the simple sum f lIV.139p . see the 
discussion of this numerical form of the Poisson identity in ref.s [I] -[3]. 

Finally, we have to face the problem of determining the non-relativistic anisotropy 
energy E}^f flIV.109jl and its mass equivalent Mfi^-'^c^ ( 1IV.114p . The situation is here the 
same as for the precedent isotropic case: if we succeed in finding the exact solution of the 
quadrupole equation (lIV.127p . the quadrupole constraint N|^|'f^ (IIV.IISP will vanish so that 
we can omit it in the wanted energy function ¥,^^"^^((3, u) which arises as the value of the 
functional f lIV.119j) - flIV.121j) on our trial configuration |<l>(r), [^'Uo(r), ^PU"'(r)|. But 
fortunately, the desired solution can be exactly constructed by some simple mathematics 
which will be presented in greater detail in App. A. The result is the following general 
form of the anisotropy energy 

Ejf = -(^)e3 ■ - (2 Pas) ■ ell{i^) (IV.142) 
as 



where the "anisotropic" potential function £pot(z^) is obviously the anisotropic counterpart 
of the "isotropic" Spoti^-') (1IV.139p - (]IV.14ip and thus is a relatively complicated (but well- 
determined) function of the variational parameter i/, see App. A. 



92 



Partial Extremalization 



Now that all physical contributions to the energy functional E[$] (1IV.119|) on our 



subset of trial configurations |<l>(r), *^pMo('")) ^^"^^"'('")| have been fixed, we arrive at the 
following result: 

E(^v) u) = 2 EUP, ^) + 4' + ^m^ , (IV.143) 

i. e. we get by use of the results f lIV.133[) - (]IV.134l) for the kinetic energy E\an and the 
results (I1V.138I) - (I1V.139I) for the "isotropic" gauge field energy E^^ plus fllV.1421) for its 
anisotropic part E^^^ : 



p2 p2 



E(iv)(/3,z/) = — (2/3aB)'-£kinM (2/3aB)-£potM (IV.144) 

ttB dB 

-^^W- — (2/3aB)-/if,f(z/) 
as 

= — {{2/3aB) ■ ei,i,,{u) - {2 l3aB) ■ etot{j^)} , 
as 

where the total potential function £:tot(^) is evidently defined through 

^rtot(//) = e^otiiy) + ■ /ifii^(^^) . (IV.145) 

It is important to observe here the fact that the total energy function E*^^^^(/3, i/) 
(1IV.144P depends on the two continuous variational parameters (3 and u but also on the 
two discrete quantum numbers i-p and The first one (i-p) is contained in the kinetic 

energy function ekinl'^) flIV-1341) and the second one C-^Ws) is associated to the discrete 
values of the quantum number iz, cf. the table on p. |8H Therefore the extremalization of 
the energy functional transcribes now to the search for the extremal points of the energy 
function E,^^^\(3, u) flIV.144p with respect to the continuous variables (3 and v whereby the 
discrete parameters i-p and ^P^n^, must be kept fixed. Thus for any allowed pair of these 
fixed parameters one obtains the corresponding binding energy of para-positronium, and 
this provides us with the possibility of studying the question of degeneracy (concerning 
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the same energy for different pairs of discrete parameters), albeit only in the present 
quadrupole approximation. 

The wanted extremal points of the energy function E'-^^\/3,u) are determined now as 
usual by the vanishing of its first derivatives, i. e. 

9E(iv)(/3,z/) 



d/3 

gE(i^)(/3,;/) 
du 



(IV. 146a) 

. (IV.146b) 



These two conditions are of a rather different degree of complication so that it is advanta- 
geous to first deal exactly with the simpler one (i. e. (]IV.146aP ) and afterwards solve the 
residual condition ( II V. 146b]) by means of an adequate numerical program. Indeed, the 
energy function E(^v)(/3,z/) flIV.144j) is a simple polynomial function of the first variational 
parameter /3 so that the first extremalization condition flIV.146a|l fixes /3 to /3* being given 
by 

2/3.a^= /"^^';\ , (IV.147) 

and if this is substituted back in the energy function E^^^\/3, v) flIV.144j) one is left with 
the problem of determining the minmal values of the reduced energy function (Ep(z/), say) 

E,M^Eav)(^..)|^^^,._jL.£|M. (,V.148) 

It is true, the denominator ekinl'^) is here a relatively simple function, cf. App. A, but 
the numerator etQt{v) f llV.MSp is given only in form of a rather complicated power series. 
Therefore one will have to use some suited numerical program for determining the local 
minima (Ep"-*^, say) of the reduced energy function Ep(z/) (IIV.148P for any possible value 
of the discrete parameters l-p and ^z "-^Vna). 
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4- Para- Spectrum 



For a brief demonstration of how the energy degeneracy (with respect to £z) of the 
spherically symmetric approximation becomes now eliminated in the quadrupole approxi- 
mation we will restrict ourselves here to the precedently discussed quantum number ip = 3 
—3 < ^2 < 3). Clearly, similar degeneracy effects must emerge for all quantum num- 
bers i'p > 1 but this will be deferred to a separate paper. For the moment, it is perhaps 
more instructive to consider the case = ±ip for all values ip (i. e. i-p = 1, 2, 3, 4, . . .). 
The point here is that the preliminary estimate of the quadrupole level splitting flIV.129p 
demonstrates that for i-p = 3 the lowest-lying energy level of the corresponding 
multiplet" is that for = ±£p. Therefore the hypothesis suggests itself that this perhaps 
might be true for the £p-multiplet due to any given value of ip. Consequently, the energy 
spectrum due to iz = for all values of ip is the adequate one in order to be compared 
to its conventional counterpart f ll.2p . We will readily plunge into this matter of question 
after first having glimpsed at the £^-degeneracy for ip = 3. 

Elimination of the iz- Degeneracy, Ip = 3 

For the search of the local minima of the reduced energy function Ep(z/) f lIV.148p one 
inserts there the kinetic energy function £:kin(^) (1IV.134|) for = 3 and the total potential 
function Etoti^) (1IV.145P which yields 

e' (£po.(i')+ '-'™3-A'f|f(i'))= 
^-'"'l'.- = + x " ■ (IV-WS) 



Here, the potential function e^otii^) is given by equation (IIV.139P and /in? by flA.19p or 



flA.20p of App.A. The coefficient *^P^m3 is linked to the possible values of 1^ according to 
the table on p. [HH and thus the four functions ¥.p{v) for lp = 3 and Iz = 0, ±1, ±2, ±3 are 
well-defined in order to determine the four corresponding binding energies lEp^|i^ as the 
local minima (with respect to the variational parameter u) of the energy function Ep(z/). 
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But observe here also that this minimahzation process is based on the simple trial 
function $(r) (lIV.13ip which has no zero. Therefore we can obtain by our extremalization 
procedure only the "groundstate" in the effective potential which itself appears as the 
sum of the centrifugal potential (~ l"^) and the electrostatic potential (due to [pUo(r) 
and ipUlll(r) in the Schrodinger like equation flIV.130|) ). This means that the principal 
quantum number [rip, say) must be linked here to the angular momentum quantum 
number through 

np = /p + l. (IV.150) 

Thus, the (non-relativistic) configurations under consideration are those being character- 
ized by the quantum numbers = 4,1^ = 3, Iz = 0, ±1, ±2, ±3. The corresponding 
results for the para-positronium energies Ep ^ are collected in the following table: 








±1 


±2 


±3 


Ej,'^ [eV] 


-0,40298. . . 


-0,39001... 


-0,37369... 


-0, 42003 . . . 


E^' [eV] 


-0,42518... 


-0,42518... 


-0,42518... 


-0,42518... 




5,2% 


8,2% 


12, 1% 


1,2% 



For an evaluation of these results one introduces the relative deviation Ap''L from the 



conventional predictions through 



A(4)| 



E, 



(4) 



e: 



(4): 



E. 



(4) 
C 



(IV. 151) 



where E^ is the conventional prediction fll.2|) for n = A 

1 



E, 



(4) 



0, 42518 . . . [eV] . 



(IV. 152) 



4aB 42 

Obviously, the deviation from the conventional result is maximal for 1^ = ±2. Indeed 
for this special value of Iz the anisotropy correction is zero, (cf. (IIV.104c]) ). and therefore 
we obtain a deviation of 12, 1% which is due to the spherically symmetric approxima- 
tion, see ref. However, for all other values of Iz one gets a smaller deviation A^p^ 
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which adopts its minimal value of 1,2% just for the maximally possible value of Iz, i.e. 
Iz = ±3 = ±/p. In this way, the conventional n^— fold degeneracy (11. 3p is eliminated as 
far as the 2;— component of angular momentum is concerned. Whether there is also an Z^- 
degeneracy for the same value of with < < ri^ — 1 necessitates an extra discussion. 
The elimination of the /2-degeneracy is illustrated also by the subsequent Fig. IV. B 



-0.3 



-0.4 



-0.5 



-0.6 



-0.7 



« = 5 



« =4 



« = 3 



I =±2 



l=±\ 

z 

1=0 

z 

I =±3 

z 



Fig. IV.B RST Elimination of the I z-Degeneracy for = 3 

The conventional energy levels Eq^ fll.2p (broken lines) are n^— fold degenerated, cf. fll.Sp . 
The RST quadrupole approximation eliminates this /^-degeneracy for the states with prin- 
cipal quantum numbers rip, 1^ = rip — 1; here = 4,1^ = 3. Presumably, the higher-order 
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approximations will bring the RST predictions even closer to their conventional counter- 
parts p.2|) (broken lines), at least for small values of the principal quantum number n^. 
For the magnitude of splitting for larger n^, see the table on p. llOll and Fig. IV. C below. 

Spectrum Due to iz = ±^7? 

The restriction = ±£p selects just one energy level from an £p-multiplet {—ip < 
< ^p); and therefore we have to carry through all the necessary intermediate steps for 
just this one multiplet term, however for all possible values of i-p (= 1, 2, 3, 4, . . .). 

The first step concerns the solution of the first-order angular system (IIV.8ap - (IIV.8bp . 
Here we can restrict ourselves to the case iz = since the commutation symmetry 
Iz =^ —^z simply entails the interchange /i?(^9) -v^ fs{^)- Thus, for £z = i-p that system 
(IIV.8ap -( ]IV.8bp becomes specialized to the following form: 



(I'd 



- tp cot ^ ■ fn{:&) = 2ip- ism (IV.153a) 
+ £pcot^9-/s(^9) = 0. (IV.153b) 



But here one is easily convinced that the second equation (lIV.153bp does admit exclusively 
the trivial solution js[^^ = which is non-singular on the 2;-axis = 0, vr). This, however, 
recasts the first equation pV.153aP to the homogeneous form 



^M!) _ £p cot ^ ■ fn{^) = , (IV.154) 
and this admits the non-trivial solution (obeying of course also flIV.17aP ) 



where the separate angular normalization condition (llV.lSaP has already been regarded. 
As a consistency test, specialize this result (1IV.155P to the case ip = 3 and find just the 
result of the table on p. |56l last lines. 
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The next step must refer to the calculation of the angular factor ^^^^'"(-(9) (lIV.104cp . 
Inserting therein the trivial function fs{'&) = together with the result (1IV.155P for fni'^) 
and carrying out explicitly the integration yields 



+ 1 



' 2 1 
COS V — - 



(IV.156) 



Indeed, for i-p = 3 one is immediately led back to the former result (lIV.104cp . last line. 
But once the angular pre-factor ^^^y4"'(7?) is known explicitly, one can proceed at once 
to the calculation of the numerical coefficients ^^h^ pV.llOap . ^^^f^ pV.llObj) . and ^^Vn^ 
OIV.llSP : the result is the following: 



ip)r 



1 


(i 


£v - 1 


80 




V + 1 


1 


(i 


iv - 1 


40 




V + 1 






iv - 1 






V + 1 



(IV.157a) 
(IV. 157b) 
(IV.157c) 



Obviously, these results do generate just those numbers of the table on p. (HU last column, 
if the quantum number i-p is again specified down to i-p = 3. But more generally, we have 
here for any value of ip (and even for any allowed value of Iz, see App.C): 



6 



ip), 



'63 



(IV. 158a) 
(IV. 158b) 



This ensures again that the general quadrupole equation flIV.126p adopts its specific form 
(IIV.127j) : and furthermore the anisotropy energy ^Si^ f lIV.128|) reads now 



{e} 
III 




(IV.159) 



The quadrupole solution ^P^A^^\r), to be used herein, is elaborated in App. A and is of 
course the same as previously used for the special case ip = 3. 
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Since that solution ^^^^'"(r) is an exact one, we can evoke again the "anisotropic" Pois- 
son identity in coefficient form (lA.lSp and thus we ultimately obtain the wanted energy 



function E^'v}(/3, z/) in the same form as previously, cf. (lIV.143p - (]IV.145p . with merely 
the coefficient of mass equivalence ^P\n^ to be replaced now by its value given by equation 
flIV.157cp . Furthermore, because of this formal equivalence of the total energy functions 
u), their partial extremalization with respect to the variational parameter works 
in the same way; and thus we end up with the residual problem of looking for the min- 
ima of the reduced energy function Ep(z/) (]IV.148|) where now not only the denominator 



£]im{^) depends upon the discrete quantum number i-p, cf. (]IV.134[) . but also the denom- 
inator Etoti'^)- In cosideration of the numerical calculations, it is convenient to express 
this dependence on ip explicitly, i. e. we finally get for the reduced energy function 



^2 



Ep{u) = Sp{u) (IV.160) 

4 as 

with the spectral function S-piv) being specified by 

The following table presents the energy spectrum E^jy' ^ rap = £p -|- 1 = 1, 2, 3, . . .) 
to be determined by the local maxima of the spectral function Sp{v) for any fixed value 
of Ip: 



E^,"^) = E^(z/,) (IV.162a) 
dEJu) 



dv 



(IV.162b) 



100 



rip 




7T.{n) r„-\n 


Ap /o 




111 V . ioza 1 


111 V . iOi|l 


111 V .iOZD 1 








2 


—1.58097 


—1.70072 


7.0 


1.71724 


3 


—0.71761 


—0.75588 


5.1 


3.86851 


4 


—0.42003 


—0.42518 


1.2 


6.21902 


5 


-0.27878 


—0.27212 


—2.4 


8.83210 


6 


-0.19969 


-0.18897 


-5.7 


11.66043 


7 


—0.15062 


r\ "1 o o o o 

—0.13883 


—8.5 


14.68005 


8 


-0.11793 


-0.10630 


-10.9 


17.87403 


9 


-0.09499 


-0.08399 


-13.1 


21.22077 


10 


-0.07824 


-0.06803 


-15.0 


24.71323 


15 


-0.03690 


-0.03024 


-22.1 


43.99860 



Obviously, the RST quadrupole approximation for the lowest-lying level (i.e. Iz = ±/p) 
yields predictions of smaller binding energy up to rip = 4 and of larger binding energy 
for Up > 4, as compared to the conventional levels E^(^^ fll.2p . The other members of 
an /p-multiplet (i.e. \lz\ < Ip) have smaller binding energy (than those due to Iz = ±/p) 
and thus do form a kind of band structure around (or in the vicinity of) the conventional 
levels, see Fig. IV. C. 
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Fig. IV.C Lowest-Lying RST Levels eS"^ 



Iz — l^. 



(TVA6M) 



For increasing principal quantum number > 4 the lowest energy levels of an /^-multiplet 



-L <lz<+l 



with Ip = Up — 1, have ever larger binding energy in comparison to 
their conventional analogues Ejf^ (IL2]). Therefore it seems that, for sufficiently large n„, 
there may occur a crossing of the conventional levels (broken lines) and RST levels (solid 
lines) . 
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Magnitude of Fine- Structure Splitting 

The elimination of the degeneracy, as being displayed by Fig. IV. B on p. [97] for 
i-p = 3, may be conceived also as a certain kind of (electrostatic) fine-sctructure splitting 
of the energy levels. For such an interpretation one would like to have an estimate for 
the splitting width for any i-p. For instance, for the special case of £j> = 3 (see the table 
on p. [96]) . one observes a maximal energy difference in the quadrupole approximation for 
those levels which have iz = 2 and = 3, namely an energy difference of (0,42003 — 
0,37369) [eV] = 0,04634 [eV]. This corresponds to 12,4% of the "spherically symmetric" 
energy of 0,37369 [eV]. Such a magnitude of the splitting seems not to be supported by 
the experimental results [6]; and this necessitates to consider the RST splitting width also 
in the higher-order approximations (i.e. octupole, etc.). But as a first rough estimate of 
that questionable splitting phenomenon it may suffice to discuss the splitting magnitude 
already in the present quadrupole approximation. For this purpose, we take iz = C-v and 
£^ = £p — 1 for a preliminary estimate of the range of the splitting. 

More concretely, we have to transcribe the former result f lIV.160p - flIV.16ip for iz = i-v 
to the present choice £^ = £p — 1. The first step for this is to look for the solution of the 
angular system f lIV.8ap -f lIV.8bp for our special choice of = £p — 1. For this choice, the 



angular system appears as 



dfsW 



- {ir - 1) cot^- fn{^) = {2ip-l)- fs{^) (IV.163a) 
+ {iv - 1) cot^- fs{^) = -fni^) . (IV.163b) 



d^ 

The solution of this coupled first-order system can be obtained by the method being 
already described in Subsect. IV. 2 and is then found to look as follows 



fnW = (2tj, - 1) ll^^ll • (sintf)''-' costf (IV.164a) 
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For a brief test of this result, one puts £-p = 3 and thus recovers both functions and 
fs{^) for £p = 3, = 2 as displayed by the table on p. [561 

Next, the angular factor ^^^"'(^9) (llV.104cl) due to the present solution fllV.164ap - 
pV.164bp must be calculated. The straight-forward integration yields the result 



{pUiiif^) = _3(2£p-lH£p-3)_ 



cos d 



(IV. 165) 



Observe here that the angular dependence of the pre-factor ^'p^J^^\'§) is always the same 
in the quadrupole approximation, namely [11] 



cos^ § — - 



(IV.166) 



with the constant -f^jp} being given by 



(IV. 167) 



For the present solution (]IV.164a|) -( 1IV.164b|) the explicit integration yields 



(IV. 168) 



2£p (£p + l) ' 

and if this is inserted in the general expression flIV.166p one just recovers the claimed 
result flIV.165p . Observe also that this latter result just agrees with the former result 
f lIV.104cjl for £p = 3,4 = 2 where ^^"'(^9) must vanish! 

Once the angular factor ■^^^'"('i?) is known, cf. (]IV.165p . one can calculate the corre- 
sponding quadrupole parameters ^^^^63 (llV.llOap . ^'^^f^ (llV.llObp and ^^Vn^ (]IV.115p . The 
straight-forward integration yields the following result 



(pl 



'1713 



1 

80 
3 

40 
1 

16 



(2 




-1) 


{iv- 


3) 






(£p 


+ 1) 




(2 


iv- 


-1) 


{h- 


3) 






{iv 


+ 1) 




(2 




-1) 


{h- 


3) 



(IV. 169a) 
(IV. 169b) 
(IV. 169c) 
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The striking feature of this result refers again to the ratios 



ip), 



6 
5 . 



'63 



(IV. 170a) 

(iv.irob) 



Obviously, this agrees with the corresponding ratios found also for the set of quantum 
numbers {i-p = 3; iz = 0, ±1, ±3}, see the table on p. IHH and it agrees with the ratios 
found for the quantum numbers {i-p; = i-v} for all ip. This suggests that these ratios 
are the same for all allowed conbinations of ip and iz in the quadrupole approximation. 

In any case, the observed ratios f lIV.170ajl -f lTV.170bj) do ensure that the general quadrupole 
equation (1IV.126P adopts again its specific form f lIV.127|) whose solution is elaborated in 
App. A. The associated quadrupole energy e}^^ (IIV.128P adopts now by use of ^^V^, 
(IIV.169a| the following form 



E^ 



III 



1 

80 



(2£p-l)(£p-3) 



n 2 



1 

16 




(IV.171) 



i2ip-l){ip-?>) 



ip yip 



+ 1) 



as 



(2aB/3)-W 



where again the "anisotropic" Poisson identity is used in its coefficient form flA.lSp . Since 
furthermore the "isotropic" gauge field energy E^ f llV.lSSp is the same as before, the 
total gauge field energy E^^ fllV.lOSp emerges also again in the same form as before, i. e. 



E- 



R 



as 



(2 as/?) ■ Etotii^) 



where £tot(i^) is now given by 



Spotiiy) + — 



{2ip-l) {ip-?,) 



(IV. 172) 



(IV.173) 



iv {iv + 1) 

Moreover, the kinetic energy can also be taken over for the present situation so that 
we ultimately are left with the problem of determining the minima of the reduced en- 
ergy function Ep(z/) flIV.160p where, however, the spectral function Spiv) is now of the 
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following form: 



|£pot(z^) + 1^ 


"(2£p 


-l)(£p-3)" 






1 








2u+l 


4 ^ 2u 







(IV. 174) 



in place of the precedent situation (lIV.16ip . 

This formal similarity with the precedent situation admits us to determine the minima 
of the reduced energy function E-p(z/) (1IV.160P by use of the same numerical method, with 
merely the spectral function S'p{v) pV.161j) to be replaced by the present one flIV.1741) . 
The corresponding results are collected in the following table: 





eJ,") [eV] 


4^^) [eV] 


A[r) % 






flIV.162aP 


m 


(IIV.151P 


(IIV.162b) 


2 


-1.65471 


-1.70072 


2.7 


1.72555 


3 


-0.68113 


-0.75588 


9.9 


3.76782 


4 


-0.37369 


-0.42518 


12.1 


5.87475 


5 


-0.24146 


-0.27212 


11.3 


8.17493 


6 


-0.17148 


-0.18897 


9.3 


10.69620 


7 


-0.12931 


-0.13883 


6.9 


13.43108 


8 


-0.10161 


-0.10630 


4.4 


16.35354 


9 


-0.08228 


-0.08399 


2.0 


19.45641 


10 


-0.06818 


-0.06803 


-0.2 


22.71757 


15 


-0.03306 


-0.03024 


-9.3 


41.08929 
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Fig. IV.D RST Levels 
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For quantum numbers = — 1) the RST levels E^'' (]IV.162a|) (solid lines) have 
lower bounding energy then their conventional counterparts E^^'^ (11.2p (broken lines) in 
the range 2 < /p < 9, see the table on precedent p. 11061 However for 1^ > 9, the RST 
binding energy E^^ becomes larger than the conventional prediction E^'' fll.2p . Thus the 



RST predictions Ep ^ have lower degree of degeneracy but are placed within some band of 
relatively narrow width around the conventional (highly degenerated) levels E'^f'' (II. 2p . 
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Summarizing all the numerical results obtained up to now, it may be sufficient to 
collect only the three levels for 1^ = 0,lp — 1, Ip (see the tables on pp. I120fl06fl01l) into 
one figure in order to attain a first impression of the RST band structure, see Fig. IV. E 
below. For low principle quantum numbers (n < 5, say) the RST energy bands are 
neatly separated and are situated in the neighborhood of the conventional predictions 
fll.2|) (broken lines). But for increasing quantum numbers (n > 5, say) the splitting of the 
degenerated conventional levels E^^ fll.2p becomes larger and larger; and it seems that 
finally the RST spectrum may adopt a chaotic character. Whether this is a consequence 
of the roughness of our approximation procedure or is a true feature of the exact RST 
spectrum must be studied separately. 
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Fig. IV.E Band Structure of the RST Levels 
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Appendix A 



Quadrupole Energy 

In order to determine the quadrupole energy fllV.1281) we first have to inquire 
into the solution ■tpM"'(r) of the respective quadrupole equation (lIV.127p . Passing here 
over again to dimensionless objects and referring to the selected trial amplitude $(r) 
fllV.131l) - fllV.132|) . that quadrupole equation reappears now in the following form: 



fmyl^l^^.e^^.,,^-^-K^r'^ ,A.la) 

dy y dy y^ 

(y = 2Pr, ^r(l/) = ^^f^-^^^^"'w) . (A.lb) 

If it were possible to find the exact solution of this equation, one could substitute it in 
the first-order anisotropy energy E^^^ (1IV.128P which reads in terms of the dimensionless 
potential correction Vlj\y) (lA.lbp as follows 

e};;^ = -^^h ■ — (2 Pan) e^U^) . (A.2) 

Here, the potential function ^potl^^) is the anisotropic counterpart of the "isotropic" £:pot(^) 
flIV.139p and reads in terms of the (dimensionless) anisotropic potential Vl^\y) flA.lbp 





Naturally, one would prefer to find the exact solution of the (dimensionless) quadrupole 
equation ( lA.laP in order to get the corresponding exact value of the quadrupole energy 
eI^^ ( IA.2p . Fortunately, this quadrupole problem allows us to work out the desired exact 



solution, albeit only in form of a power series. 
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The Exact Quadrupole Solution 



The first observation for finding the desired solution refers to the fact that the quadrupole 
equation (lA.lap can be recast to the following form: 



1 d 

dy 



y ■ -r 



dy \ y^ 



-5y^^-^-e-y . (A.4) 



Integrating this from arbitrary y up to infinity {y ^ oo) yields for arbitrary z/ (> |) 



For further integrating this intermediate result we have to observe the boundary conditions 

limP'"(?/) = const. ■ y^ (A.6a) 
limP:"(l/) = ^ (A.6b) 

3/— s>oo U'-' 

which may easily be deduced from the homogeneous version of the quadrupole equation 
flA.4p . The point here is that the solutions V^J\y) of the inhomogeneous equation flA.4p 



(or flA.lap . resp.) must adopt the form of a solution of the homogeneous equation in those 
regions of three-space where the inhomogenity (~ y'^'^~^ ■ e~^) is practically zero, i.e. for 
y ^ oo and also for ?/ — )■ for sufficiently large values of the power u. 

With regard to this circumstance concerning the boundary conditions, one may further 
integrate (by parts) the intermediate result (lA.Sp in order to find by means of some simple 
manipulations 



°o 2u-l+n 

T{2u-l).y'e-yYl 



r(2z/ + n) 



Here, it remains to determine the integration constant which emerges in connection with 
the solution of the homogeneous equation (~ y^). This problem can be dealt with very 
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conveniently by specializing that general result (]A.7|) to integer values of 2i>. Indeed, it is 
easily seen by means of some simple manipulations that for this case that general result 
(1A.7P adopts the following form 



V'J\y) const. ■y2 + (2z/ + 3)! 1""=° (A.8) 



2v-2 

y 



n=0 



But here it is now evident that the required boundary condition (lA.6bp is violated, ex- 
cept the integration constant in equation (1A.7P adopts the special value const. ^ {2u — 2)! 
so that the "integer-valued" solution (lA.Sp appears now as 

1 _ p-y . y^2;/+3 

piii(y) = {2u + 3)! ^3 (A.9) 



y 



2*^-2 n 

y 



n=0 



or, by generalization of this result, the solution for general (but real- valued) v is 

p-y °° 2v+i+n 

Prfe) = r(2. + 4)^.g^^A__ (A.10) 

f °o 2u-l+n 



T{2u + n) 

n=0 ^ ' 



Obviously, there appears to be a certain peculiarity with this result; namely concerning 
the behaviour of V]^\y) at the origin (?/—)■ 0). This can be most easily seen by inspecting 
a little bit more closely the "integer- valued" result f lA.Op . Indeed, this form of V^J\y) can 
be recast to the following shape which more clearly displays the behaviour for y — )■ 0: 

(2z/ + 3)! (2z/-2)! 



oo 



V^J\y) = {2u-2)\y' + e-y.Y^ 



n=0 



(2z/ + 4 + n)! (2z/-l + n)! 



_y2u+L+n_ (A.ll) 



For instance, for z/ = 3 this looks as follows in the vicinity of the origin [y — > 0) 

V^^\y)^2Ay'-y- + l^y' + ... (A.12) 
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The point here is that for such a differential equation as (lA.4p . where the inhomogeneous 
term (source term) behaves as (~ y^^'^) near the origin, one expects that the solution 
V]^\y) should then behave there as ~ y^^'^^ because the differential equation is of second 
order. However, this expectation is not validated by our result (lA.lip since (for z/ > 1) 
there always dominates the "homogeneous" solution (~ y^) in the neighbourhood of the 
origin. Actually, it is the second (i. e. exponential) term in flA.ll|) which goes like ~ y"^""^^ 
and thus meets with the original expectation. 

Observe here also that one cannot simply subtract off the "homogeneous" solution 
(~ y'^) since in this case one would spoil the required boundary condition ( IA.6b|) at 
infinity! This is easily seen by returning to the equivalent form (1A.9P which is numerically 
identical to the latter form (lA.lip of V]^\y). Finally, let us also mention that the latter 
integer- valued form flA.lip of the solution V]^\y) admits the obvious generalization to 
arbitrary (but real- valued) u, i. e. 

oo 

V]j\y) = T{2u -l)-y' + e"^ ■ J^pM v'"^''"' (A.13a) 

n=0 

, , r(2z/ + 4) r(2z/- 1) , , 

^ r(2z/ + 5 + n) T{2u + n) ^ ^ 

This is a complementary form of the precedent result flA.lOp or flA.Qp . resp. Indeed, 
whereas flA.Qp clearly demonstrates the required boundary condition flA.6bp at infinity, 
the present result (]A.13aP validates the condition ( IA.6aP near the origin. 

As a check of this result, one considers the function Vlj\y) (lA.lSaP not the result but 
rather an ansatz with unknown coefficients p„ which depend of course on the variational 
parameter z/. Substituting this ansatz into the quadrupole equation (lA.lap . one gets the 
following recurrence formula for the determination of the coefficients pri: 

2{2u + 3 + n)- pn+i - Pn 
(2z/ -1-4-1- n){2h' + 3 + n) — 6 
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with the two lowest-order coefficients being given by 

Pq = —-, ^ (A. 15a) 

^" (2z/ + 2)(2z/ + 1) - 6 ^ ' 

^^=(2. + 3)(2. + 2)-6-^°- ^^-^'^^ 

From consistency reasons, the former coefficients Pn{i^) flA.lSbl) must then turn out as the 
desired solutions of the present recurrence problem flA.14p - flA.15bp . Indeed, this claim 



may easily be verified by a few simple algebraic manipulations. 

Quadrupole Energy 



For calculating now explicitly the quadrupole energy 

dD we simplify again the 



numerical computation by resorting to the quadrupole identity fllV.llSp . This means that 
we do not substitute our solution V]^\y) in the "anisotropic" potential function ^potl^^) 
flA.Sp in order to obtain the anisotropy energy E^^^ flA.2p . but rather we substitute the 
solution V^J\y) into the "anisotropic" mass equivalent M|^|'|"'^c^ flIV.114p which reads in the 
dimensionless notation 



d^^c^ = -(P)m3 — {2PaB) ■ ^S?(z^) (A.16) 



2 



^111 



with the mass-equivalent function yU^n (z/) being defined through 







This procedure is legitimate because our general quadrupole solution flA.lOp is an exact 
solution of the quadrupole equation f lA.lap so that the "anisotropic" Poisson identity 
fllV.llSp is actually satisfied. Of course, this is the same effect as mentioned already in 
connection with the spherically symmetric approximation, cf. the remarks below equation 
flIV.139p . Consequently, substituting both results flA.2p and f lA.16p into the "anisotropic" 
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Poisson identity (1IV.118|) we obtain the numerical identity of both energy coefficients 



-III / 

'pot I 



flA3|) and /if||^(z/) flATT]) : 



(A.18) 



which may be used as a consistency check for the numerical correctness of the computer 
programs. 

The explicit computation of ^\\^ (]A.17|) may be performed now by means of either 
the precedent form (lA.lOl) of the solution V]^\y) or also by means of its second form 
flA.13ap - (IA.13bp . In the ffist case ([OOj) we find 



/ifi^(^) 



T{2v + A) 
r(2i/ + 2)2 



n=0 



+ i^.r(2. + 4) 



1 r(4z/ + 3 + n) 
24^+3+" ' r(2z/ + 5 + n)_ 

Vi2v - 1 



(A.19) 



r(2z/ + 2)2 
and in the second case 



,{e}. ,_ r(2z. + 4)-r(2z.-i) 



r(2z/ + 2) 



n=0 



1 r(4i/ + 3 + ra) 
2iv+i+n ' r(2z/ + n) 



r(2z/ + 2)2 



Of course, both results for /i[|*|^(z/) must be identical (thus providing a confidence test for 
the corresponding computer programs). For half-integer values of z/ (-w = 1, |, 2, |, . . .) 
the mass-equivalent function /i|-|'|^(z/) flA.20p can also be represented by a finite sum: 

(2z/ + 3)(2z/ + 2) 



r(2z/ + 2)2 



n=0 



r(4z/ + 3 + n) 

24v-|-3+n 



(A.20) 



(2Z/+1)! 



with 



E 



{(2z/-2)!-U-5-K^} 



(2z/ + 3 + - m)! 



(n + Q)\ ^ 2^v+^^n-m 
n=0 ^ ' m=0 



m — m ! 



E 

n=0 



1 (2z/-2 + n)! 

n! 22^^-1+" ' (2z/-2)! 



(A.21) 

(A.22a) 
(A.22b) 



see Appendix E of ref. [Ij. Inserting now this result flA.20p via the Poisson identity 
( lA.lSp in the expression f lA.2p for the anisotropy energy yields the anisotropic gauge field 
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contribution to the energy functional E[$] ( 1IV.119|) which thus adopts the simple form 
(II V. 144|) because all constraints are automatically satisfied. 
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Appendix B 



Special Case iz = 0, Vi-p 

The situation for vanishing quantum number iz and arbitrary i-p (£p = 0, 1, 2, 3, . . .) 
represents a special case insofar as the relevant quadrupole objects ^^^ko^'d), ■>^^M"'('i9), ^^h^, 
(p)f^^ 'yP\n^^ Ei^i'f-'^ are all independent of the quantum number (.-p. 

In order to validate this assertion, one first turns to the angular density ^'P^k^^'d) 

47r sin v 



which for the product ansatz may be separately normalized to unity, cf. fllV.lSap . Now 



differentiate this object and find by use of the first-order angular system (IIV.8ap - flIV.8bp 



' + cot ^ ■ WA:o(^) = {2/«(^)/5(^) + cot ^ [fU^) - fli^)] } . (B.2) 

For the presently considered case 4 = 0, this becomes simplified to 

^^^ + cot^.W;toW = (B.3) 
with the solution (normalized according to the prescription OlV.lSap ) 

ZTT'^ sm tf 

Observe here that our demand iz = eliminates completely the quantum number i-p from 
the equation (1B.2P so that the resulting solution flB.4p is actually independent of ip\ Of 
course, this independency must then be transferred to all the other objects which derive 
from the angular density ^^^kQ^i!)). 

The first one of these objects is the anisotropic potential correction ■^^^'"(■(9) (lIV.104cp . 
being defined quite generally by 



(B.5) 
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or equivalently by 



cos — 



(B.6) 



with the numerical constant -f^jp} being defined through 



= / dn cos2i?-^?'>A:o(t9) . 



(B.7) 



This constant is easily calculated for the present angular density 'f^^/co(i?) flB.4p and yields 
the result 

= I ■ (B-8) 
Thus the anisotropy correction ^^^^'"(■i?) flB.6p becomes 



cos -i? — - 



(B.9) 



The interesting point with this result is now that it does hold for any value of the quantum 



number ip, see for instance the case L_ 



in equation flIV.104cp . 



The next quadrupole objects are the numerical constants (^^63 f llV.llOap . (P^/g f llV.llObp 
and ^^Vn^ (lIV.llSp . Inserting in these definitions the present potential correction 
f lB.9|) and the angular density ^^^ko^i!)) ( 1B.4P yields for these constants 

1 

80 
3 



''=80 



•' 40 
1 



16 



(B.lOa) 
(B.lOb) 
(B.lOc) 



Observe here again that this result does hold for all values of the quantum number i-p, 
not only for i-p = 3 being displayed by the table on p. [HH This conclusion must then be 
true also for the anisotropy energy eI^^ f lIV.109|) or ( 1IV.142p . resp. 



{e} 
III 



16 qb 



(B.ii: 
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Accordingly, the total electrostatic energy E^^ f llV.lOSp becomes 



{c} 



R 



(2/3aB) ■ etot{j^) 

as 

which finally lets appear the energy function Ep(z/) flIV.148p as 

„2 



(B.12a) 
(B.12b) 



Ep(z/) 



4a 



B 



■ S'p{v 



(B.13) 



with the specific spectral function S-p^u) being given now by 



1 + £i 

2v+l V 4 2i/ 



(B.14) 



Here, the peculiarity due to the special case ^2 = becomes now manifest: in contrast to 
the general situation (f^ ^ 0), e.g. 0IV.16ip or f lIV.174p . the numerator is independent 
of the quantum number £p which thus infiuences exclusively the denominator being rep- 
resented by the kinetic energy i?kin (1IV.133p -( 1IV.134p . Or, equivalently, the electrostatic 
interaction E^^ (IB.12ap - (lB.12bp does not depend explicitly upon the quantum number 
i-p of angular momentum. From this one concludes that possibly the corresponding en- 
ergy spectrum, due to £z = 0, may show some characteristic features missing for the other 
values {^z 7^ 0) of angular momentum. The following table represents the corresponding 
energy spectrum to be calculated again by the method of partial extremalization described 
below equation (1IV.148P : 
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h = 0, VI 





[ev] 


[eV] 

m 


A[r) % 




(lIV.162a|) 


flIV.151) 


flIV.162b) 








2 


-1.65471 


-1.70072 


2.7 


1.72555 


3 


-0.71761 


-0.75588 


5.1 


3.86851 


4 


-0.40299 


-0.42518 


5.2 


6.10402 


5 


-0.25889 


-0.27212 


4.9 


8.50058 


6 


-0.18075 


-0.18897 


4.3 


11.03747 


7 


-0.13356 


-0.13883 


3.8 


13.69725 


8 


-0.10282 


-0.10630 


3.3 


16.47923 


9 


-0.08167 


-0.08399 


2.8 


19.36094 


10 


-0.06648 


-0.06803 


2.3 


22.34116 


15 


-0.03013 


-0.03024 


0.3 


38.50906 


20 


-0.01719 


-0.01701 


-1.0 


56.43818 



The numerical evidence says here that the RST energy predictions (second row) 
are placed above the conventional results E^^^ (third row) up to principal quantum num- 
bers rip ~ 15. For > 15 the RST binding energy (— E^'') becomes larger than its 
conventional counterpart {—Eq'^). As far as the present numerical results are concerned, 
the RST energies due to = do mark the upper edge of the band of non-degenerate 
RST levels due to the same quantum numbers n^, 1^ (= — 1), see Fig.s B.I and IV. E. 
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n = l5 
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■5 
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Fig. B.I RST Levels due to 1^ = 0, V/^, 

The energy predictions flRTSjl - flRlil) due to = 0, V/^ (solid lines) appear as 
being situated closest to their conventional counterparts E^^"* f ll.2p within the bands of 
non-degenerate RST levels; and furthermore these levels (due to Iz = 0) do form the 
upper boundary of the /p-multiplets defined by — < < +lp- For > 15, this upper 
boundary slides below the conventional levels E^^ (II. 2p . 
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Appendix C 



Universality of the Quadrupole Ratios 



and ^ 



'63 



63 



At various occasions one observes the emergence of the following ratios: 



6 



(C.la) 
(C.lb) 



e.g. for £p = 3 and = 0, 1,3 (see the table on p. |H1]) or for = (see equations 
f llV.158ap - (IIV.158bp ) or also for 4 = ± {^v - 1), see equations (llV.170ap - (IIV.170bp : and 
finally for = 0, V£p (see equations fIB.lOaP -f lB.lOcP ). Thus, one is led to the supposition 
that these ratios fIC.lap - flC.lbp must hold for any allowed combination of the quantum 
numbers iz and i-p. 

Indeed, it is a simple matter to convince oneself of the correctness of that supposition. 
The point of departure is the general form of the anisotropic correction potential ^^^'"(i?) 
(IIV.166P which we rewrite here as 



■ 2 1 
COS V — - 



(C.2) 



Wyl"'(t9) = am ■ 
with the constant am being defined by 

a,M = -^ (l-3ir|y . (C.3) 
Now insert this general form of ^"p^A^^^i!)) flC.2l) in the definition fllV.llOaP of "^^^63 and find 

(C.4) 

Next, insert the general form flC.2p also in the definition of '■^^/s fllV.llObI) and find 



(P)p, - _ . fl2 



15 
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''III 



(C.5) 



Thus, both general results (lC.4|) -( !cr5|) validate the general character of the first ratio 

The case of the second ratio (]C.lb|) is a little bit more complicated. But inserting 
again the general definition of ■"^^^'"(^9), see equation (lIV.104cp or (IB.Sp . resp., in the 
definition of ^^Vns flIV.1151) yields by some simple manipulations and the normalization 
condition fllV.lSal) for ^P^kQ^i!)) the following form of ^P\ns: 



'1713 



(C.6) 



After having carried out here the azimuthal integrations, this result reappears as 



- ( (^Vna + - 
3\ 2 



dn cos2i?^%(^?) 



1 2 



+ 



dn smH^%o{^) 



(C.7) 



Now observe here the definition of -^'p} flIV.167p . together with the normalization condi- 
tion (llV.lSaP so that this result reads in terms of K^^^j 



I.e. 



3^^m _1 



or in terms of am flC.3l) 



ip), 



m3 = \ -am 



(C. 



(C.9) 



(C.IO) 



Now combining this with the precedent result for ^^^^63 flC.4p yields just the claimed ratio 

dcibD. 



123 



Appendix D 

Angular- Momentum Algebra for Para-Positronium 



There are certain similarities between the conventional eigenvalue problem (]IV.16ap - 
(lIV.lGbp and its RST analogue (lIV.8ap - (]IV.8b|) . The quantum numbers m do 



occur in boath approaches in the same sense (namely as eigenvalues of L^), but the other 
quantum number i-p seems to emerge in a somewhat different way than is the case with 
the conventional i which is due to the square f llV.lGap . Therefore one may raise the 
question whether in the RST formalism there is also an operator whose eigenvalues are 
just given by the values of that second quantum numer i-p. 

Subsequently, we will inquire into this question and we will identify the wanted op- 
erator which has ip as its eigenvalue (see equation (ID.14p below.) In the course of this 
inquiry it should become clear that the RST algebra of angular momentum deserves fur- 
ther elaboration beyond the present extent. 

Eigenvalue Problem for Angular Momentum 

In order to deal with the angular momentum algebra as generally as possible, we 
restart with the original eigenvalue system (lIV.2ap - (lIV.2dp . But in place of applying 



the somewhat special product ansatz f lIV.4ap -( ]IV.4bp . we prefer now the following more 
general transformation: 

ir^,(r> = (D.la) 

V r sin 

(..5,(0 = . (D.lb) 

V r sin 



This recasts the first half flIV.2ap - flIV.2bp of the eigenvalue system f lIV.2ap -f irV2dl) to the 
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new shape 



dr r d(j) 

— cott} ■ 
r 0(h 



+ (D.2a) 



M+M, 



dr 



r 



90 



1 9(P>^+(r) 



r 



9^9 



and similarly for the second half flIV.2cp - flIV.2dp 



d(pyjz_{r) ^ 1 

dr 



r r 0(p 

1 (9(p)5-(f) ^ ^ coU9 (9(P)5- _ M-M, ^ (p)^ ^ 
r dd r dd) h 



(D.2b) 



(D.3a) 



dr r r dd) 



+ + ' ^ = -P, ^ ■ '^+ ^ 

r dv r dd) n 



(D.3b) 



The next substitution is again a certain generalization of the former case f llV.Gap - 
f lIV.6bp . namely we try now a product ansatz for the new amplitude fields ^P^±{r),^P^S±{r) 
of the following form 



■ <l_(r) 



(D.4a) 
(D.4b) 
(D.4c) 
(D.4d) 



When this product ansatz is inserted in the modified eigenvalue system flD.2ap - flD.3bp . 
there occurs a separation of this three-dimensional system in a two-dimensional angular 
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problem in terms of the angular coordinates "i?, and in a one- dimensional radial problem 
in terms of the variable r. The angular subsystem looks as follows: 

^+^cot^-^) fM<P)+^^^^^ = i^.fsi^,<P) (D.5a) 

^ - ^ cot ^ ■ ^) fsi^, 0) + ^ ^^^^ = -iv ■ fni^, 0) • (D.5b) 

Here it is easy to see that this angular system represents the immediate generalization 
of the former, more special system f lIV.8ap - (llV.8bp . since the latter system emerges from 
the present one (ID.5ap - (ID.5bp by simply putting 

/«(79,0) = e^^^'^-/«(^) (D.6a) 
/5(^,0) = e^'^*-/5(^). (D.6b) 

Obviously, the present more general system flD.5ap - flD.5bp is nothing else than the eigen- 
value problem of angular momentum for spin-zero RST systems {sp = 0). This may be 
recast also in a more concise form, i. e. 

e-"^L+ Mi^, 0) - 4 fs{i^, 4>)=iv fs{^, 4>) (D.7a) 
e'^L. fs{^, 0) + /^(^, 0) = £p ■ fni^, 0) , (D.7b) 

with the usual definitions (up to h) of the angular momentum operators 

= (4 + ^Ly) = { ^ + ^ cot } (D-8a) 

L_ = (4-^4) = -e-^^ -^^°t^^} (D.8b) 
1 d 

^ - — . (D.Sc) 

Concerning now the remaining radial problem one finds oneself being left with {spher- 
ically symmetic approximation) 

±^ - . r - [^Uo r ■ $_ r = ^ — ^c ■ $_ r p.9a 

dr r n 

+ $_(r) + ^^Uo(r) ■ $+(r) = c ■ $+(r) , (D.9b) 



dr r h 
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but this is identical to the former system (]IV.9a|) - (]IV.9b[) and thus gives no new informa- 
tion. 



Second-Order Form of the Angular-Momentum Eigenvalue Problem 



The special second-order differential equations (IIV.17ap - (lIV.17bp do hint at the possi- 
bility that also the generalized form (ID.5ap -( lD.5b|) could be recast in second-order form. 
In order to attain this goal, it is favourable to transcribe the present eigenvalue problem 
of angular momentum flD.5al) - flD.5b|) to the two-dimensional spin space. This works as 
follows: First, define three operators Qx,Qy,Qz as follows 



d 

Qx = i cot i!} — = sin (f) ■ Ly + cos (j) ■ 



dcj) 



d 



Qv 4= —i 7777 = COS ■ — sin ■ 
ov 



Qz = -i 



d 



(D.lOa) 

(D.lOb) 
(D.lOc) 



and find the corresponding commutation relations as 



Qxi Qy 
Qyi Qz 



sin^ i) 

Qz) Q 



Qz 







(D.lla) 
(D.llb) 



Next, define the operator Q, to act in the suitably defined Pauli spinor bundle over the 
two-sphere S"^ by the following self-evident construction: 



Q. = QxC^x + Qy<yy + Qz<yz = (lx + Q4J + kz 



(D.12) 



and furthermore collect both functions //{("i?, 0) and 0) to a Pauli two-spinor f (i? 
i. e. in short-hand notation 



f(^9,0) 



/5(^,0) 



(D.13) 
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And now convince yourself that the eigenvalue system (lD.5a|) - (lD.5b[) transcribes to the 
two-spinor form as 

(if(^,0)=£p-f(^9,0). (D.14) 

This reveals the quantum number i-p as the eigenvalue of the operator Q,. Finally multiply 
this eigenvalue equation (lD.14p from the left by the "quaternion" Q, (1D.12P and thus find 
its desired second-order form as 

Q,2f(^,0)=4.f(^,0). (D.15) 
Complementing finally this eigenvalue equation by 

2,f(^,0)=4-f(^,0) (D.16) 

it becomes evident that both eigenvalue euquations (ID.lSp and (1D.16P are just the RST 
analogues of the conventional eigenvalue problem fllV. 16al) - flIV. 16bl) . Of course, the si- 
multaneous presence of two eigenvalues (i.e. d-p and is a consequence of the fact that 
the corresponding operators do commute 

[i.,Q,]=0 (D.17) 

so that both eigenvalue equations (1D.14I) and (ID. 161) can simultaneously be valid. Observe 
here that the first eigenvalue equation (1D.14p plays the part of a kind of "square root" of 
the equation (ID.lSp . 

The desired second-order form of the RST eigenvalue problem for angular momentum 
emerges now simply by explicitly writing down the squared operator Q,^: 

1^ ov^ sm d(j)^ J sm -d ocj) 

so that the equation flD.lSp reads in explicit component form 

+ —T-, in-. fs{^, 0) - TTT /5(^, 0) = 4 ■ fs{^, 0) . (D.19b) 



d-t}"^ sin d dcj)"^ J sin d 
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Observe here that the squared operator Ci^ (ID.ISP is diagonal (provided one uses the 
standard form of the Pauh matrices <7x,<Jy,crz), and this ensures the result that the second- 
order equations (lD.19ap -f lD.19b|) appear decoupled. Moreover, trying the product form 
(ID.6ap -( fD.6bp we are left with the simpler system 



sm V 



aw^ sm v 



which is nothing else than the former second-order system (lIV.17ap - (lIV.17b|) . 



(D.20a) 
(D.20b) 



Ladder Operators 



In order to find the complete set of states due to an £p-multiplet (i. e. for —ip < iz ^ 
+i'p) it is very instructive to first collect the former results in a table (next page). 

On principle, one could find all the missing members (...) of such a "ladder" of 
angular states (due to the same value of ip) by solving the coupled first-order system 
f lD.5ap -( lD.5bp separately for any member and then applying the normalization condition 
f llV.lSap . But it would be surely more efficient to possess an algebraic algorithm which is 
able to generate the neighbouring states ^^'^^fR,e^±i, ^^'^^fs,ez±i from some point of departure 
'■^^^fR/z, ^^^^fs,e.- Observe here that the symmetry fji{i3) fs{i^), fs{^) -fni'^) 
mentioned below equation fllV.lTbp induces for = a certain peculiarity, namely the 
occurence of two angular states: the upper half of the line due to £2 = refers formally 
to iz = 0+ and the lower half to £2 = 0_. Indeed it is easy to see that both solutions do 
obey the coupled first-order system f lIV.8ap -( IIV.8bp for iz = 0. This double solution for 
iz = is needed below when we generate the angular states due to iz ^ from the states 
with i, = 0. 
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4=; 






Ct, 

I1IV.155D 







- 1 

(IIV. 164a )-(IIV. 164b) 


. / ' {2^v 1) ■ So^^^li; ■ e^^'^-^^-^fsin ^Y^~^ COS ^9 










(App. B) 


. /I . cosf^-o-i?) 


-^/^sinff-Di?) 


■sin(£p^9) 


■cos(£p^) 








-{h - 1) 


/ 


1 

^2 (2£^ 1) . . e*{i-^p)</'(sin^)^p-icost9 










General Ladder Arrangement of the angular eigenfunctions 



The starting point for obtaining the states with £^ > is then the solution due to 
= 0+ (upper half) and for obtaining the states with £z < ^ one starts from the angular 
state with = 0_ (lower half): 

4 = 0+: /^(^) = y|.cos(£p^9), fs{^) = -^sm{ir^) (D.21a) 

4 = 0_: /«(^) = y|-sin(£p^9), fs{^) = -^cos{ir^) . (D.21b) 

In order to elaborate this program, it is now self- suggesting to define two kinds of 
ladder operators, and 5*+, by 

R^ = e"^\ — + i cot ^— } (D.22a) 



d'd d(j) 

^+ = <( A + ^ cot \ e'^ = e-'^R+e''^ . (D.22b) 



di!) d(j) 

Furthermore, it is convenient to introduce an economical notation for the angular eigen- 
functions /ij('(9, 0) and fs{'&, 0) so that the associated quantum numbers i-p and 4 become 
immediately evident: ^^'^^^fn/^i'd, (p) and ^^''^^fs/^{'d, (J)). With this arrangement, we easily 
realize for the upper end of the above table that the highest-order state (i. e. iz = ^v) 
becomes "annihilated" under the action of the ladder operator flD.22aP : 

Since the highest-order state ^^^^fsz-pii!^, <P) of the ^-ladder is already zero, the correspond- 
ing relation for this state ^^'^^fs,e-p{^, 0) is trivially evident 

S+^'^%er = 0- (D-24) 



Thus, the combined action (1D.23|) plus (1D.24|) of both operators -R+ and 5*+ annihilates 



the highest-order state of any £p-multiplet. 
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Turning now to the next lower states (i. e. iz = iv — 1) of the table, i. e. 



we find again by subjecting those states to the action of the ascendive operators -R+ and 

K ''-^inu-^ = -^-^^'-^fR,e^ (D.26a) 

5+(^-)/5,^,_i = = (^-)/v, . (D.26b) 

From these results one becomes tempted to suppose that the operators R+ and 5+ do 
act as "ladder operators" which generate the next higher state (due to + 1) from the 
precedent one (due to i^) according to the recipe 

= C^iv, 4) ■ (D.27a) 

with the normalization constants C+(£p,4) and D^{i-p,iz) being still to be determined. 
Their special values for = £p — 1 may be already read off from the equations (lD.26ap - 
f lD.26bj) as 

C+{irJv-l) = ^ (D.28a) 

D+{i'p,iv-l) = 0- (D.28b) 



Let us test our general construction flD.27al) - flD.27bl) for those special states which 
have iz = 0,\/i-p (see the table on p.(?)) 

^'^^fR,oA^, 0) = y|cos(£p^9) (D.29a) 

^'"Vs,o+(^, 0) = -APsin(£p^9) . (D.29b) 
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For this situation, one finds by straight forward calculation for the action of the operators 

R^^'^^fR,oA^A) = -^e"t'sm{ev^) (D.30a) 

S+ ^^^^/5,o+(^?, 0) = -^e'^ ■ cos(£p7?) - cot ^9 ■ sin(£pt9)} . (D.SOb) 

Comparing this to the general action of both ladder operators /2+ and 5*+ flD.27ap - (ID.27bp 
we find that the angular functions due to = 1 must look as follows: 



^'"-^fR^i^A) = -]l- c_^^^^ 0) ■ e^^sin(£p^) (D.31a) 
^'^V5,i(^, 4>) = -\h n fl n^ ■ ■ cos(£p^) - cot t9 ■ sin(£pt9)} . (D.Slb) 

V 77 L)+[t'p, (Jj 

But these two angular functions ^^^^/ra and ^^^''fs,! flD.31ap - flD.31bp must now obey the 
coupled angular system flD.5ap - flD.5bp for = 1 without producing any contradiction; 
and this demand can actually be satisfied provided the constants C+(£p, 0) and D^{£-p, 0) 
are related by 

D^iir, 0+) = ■ C^ir, 0+) . (D.32) 

The symmetry operation z ^ —z {i^ ^ —^z) results in the substitution (p) =^ 
fs{^,4>), fs{'d,(f>) =^ —fR{d,(j)) (take the complex conjugate of the system flD.5a[) - f]D.5bp 
or recall the remark made below equation fllV.lTbp ): and therefore one may introduce 
the descendive operators -R_ and S_ quite analogously to the precedent ascendive case 
flD.22aD - flD.22bD through 

d_ 

Here it is again obvious that these descendive operators do act on the lowest-order states 
of an £p-multiplet quite analogously to the ascendive cases flD.23p - flD.24p for the highest- 



S_^- e-'^ \^-icot^^} = - e-'^ R- e'^ . (D.33b) 
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order states, i. e. 

R. ^'^^fR,-e^ = (D.34a) 

5_ = . (D.34b) 

And the lowest-order states of an £p-multiplet are again generated by 

Rj'^^fR,i.i^ = = (D.35a) 

SJ'-%i-e. = ■ , (D.35b) 



quite analogously to the ascending situation (lD.26a[) - (lD.26b|) . 



Clearly this analogy suggests to define the action of the descendive operators -R_ , 
quite generally through 

Rj'^^i^ = C^ivJz) ■ ^'^W-i (D-36a) 
SJ'^%e^ = D4iT., 4) ■ ^'"^5^-1 , (D-36b) 

with the constants C_ and D_ to be determined later on. The special values for 4 = l—iv 
may be deduced again already from the equations (lD.35ap - (]D.35bp : 

C_(^p,l = (D.37a) 

np 1 

D.{ir,l-iv) = -^=^. (D.37b) 

Recall here the analogous relations f lD.28ap - flD.28bp which apply to the angular states 
with positive values of (i- e. the upper end of the ladders). Furthermore, one deduces 
also from the survey table on p. 11301 by straightforward application of R- and S- the 
following descendive relation for the i?-ladder 

C_(£p,£p) = V^ (D-38) 
and similarly the ascendive relation for the S-ladder 
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From reasons of symmetry {iz ^ —^z) one wishes again to test the consistency of the 
ladder ansatz (]D.36ap - (]D.36b|) for the special case iz = 0, V^p. This says that we take 
from the table on p. (?) the two angular functions 

Vi?,o_ 4>) = \^- sm(£pt?) (D.40a) 

V TT 

^'^^5,0. = vf ■ cos(£p^9) (D.40b) 



and let now act the descendive operators R_ and S_ f lD.33ap - (ID.33bp on these functions 



R^^^fRfl. = e-*"^ {iv ■ cos{iv^) - cot ■(9 ■ sm{iv^)} (D.41a) 

Sy^-^o. = ^Jl^~"^ {^v ■ sin(£p^9)} . (D.41b) 

On the other hand, the ladder construction (lD.36ap -f lD.36b|) identifies here the right- 
hand sides with the angular eigenfunctions '^^'^^f^.-x and ^^'^^fs^-i which thereby become 
represented by 

^'^^fR,-i{^l <P) = -/f c^^^ ■ cos(£pt9) - cot ^ sm{ir^)} (D.42a) 
2 £p ■ e"^* 



(^-)/5,-i(^, 0) = ^J-^-Ji:^^ ■ M^v^) ■ (D-42b) 

But if our ladder hypothesis is correct, these two angular functions must obey the coupled 
first-order system flD.5al) - flD.5b|) which for the present case (4 = 1) appears as 



+ cot ^9 ■ = {iv - 1) ■ ('"^5,-1 (D.43a) 



cot^- (^^)/5,-i = -{ip + 1) ■ • (D.43b) 



Inserting here both functions '^^'^^fn-i and ^^'^^fs-i (]D.42ap - (]D.42bp leads us without any 
contradictions to the relation 

D_(£p,0_) = --i^ .C_(£p,0_) , (D.44) 



which of course is the counterpart of (ID.32P for negative values of iz-, i-e. for the lower 
ends of both ladders. 
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Ladder Formalism for i-p = 3 

In order to become somewhat more acquainted with the pecuharities of the present 
algebra of angular momentum it may be instructive to explicitly write down the action of 
the ascendive and descendive operators on the members of an £p-multiplet. For £-p = 3, 
the members are already known (see the table on p. [84]) . so that by straightforward 
calculation one can determine the constants C+,D+ flD.28ap - (ID.28b|) as well as C7__,D__ 
(ID.37ap - flD.37b|) . Turning first to the "i?-ladder", one finds the following results (see Fig. 



D.I): 

The upper edge (€-p = 3, iz = 3) is given by 

('Vi?,3(^, 0) = a/P e^^^ sin^ ^ , (D.45) 

V OTT 

see equation fllV.lSSp for £p = = 3; and consequently the ascendive operator 
f lD.22aP annihilates this angular state ( ID.45p . which is of course only a special case of the 



general relation ( ID.23p . Next, consider the lower-order state ^^^fR,2{'&, 



and find 



^^Vi? 2(^9, <!)) = \hr- e^"^ sin^ (D.46) 
V 37r 



^-f(%2 = -^(%3 (D.47) 



which then yields for the transition constant C+(3,2) flD.27aP 



C+(3,2) = --|, (D.48) 

as a specialization of the general relation f lD.28ap . Thus the upper edge of the i?-ladder for 
£j> = 3 is found to be in agreement with the former preliminary remarks. This agreement 
does also apply to the lower edge { i-p = 3, = 3). Indeed, the lowest-order non-trivial 
angular function is here 
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R-ladder 





R+ 



;,=+3 (3)/j^_3 



^+(3.2)=-^ 



R+ 



R,2 



C+(3,l)=-4 



C+(3,0+)=-|V3 



R+ 



4 



R4 



=+2 



i^=-l 



C+(3-2)=-V^ 



R+ 



(^)/r,-2 



C_(3,3)=V6 



C_ (3,2)=x/l0 



(^Vr,i 



R- 



C_(3,l)=2^ 



R,0± 



C_(3,0_)=2v^ 



(^Vr,-i 



C_(3,-1)=a/I0 



(^Vr,-2 



C-(3,-2)=0 



R- 





Fig. D.l Ladder Operations for ^^*"^/r,u(i9, 0), Ip = 3 
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and this angular function becomes annihilated under the action of the descendive 
operator _R_ (]D.35a[) . i.e. 

R. (%_2 = , (D.50) 

as demanded by the general relation flD.34ap . cf. also flD.37ap . 

But now a somewhat more subtle point is to be discussed, namely the occurence of 
two different angular states for iz = 0, see the table on p. 11301 

^^Vr,o+ = cos(3^9) = -4:sin^^) cos'd (D.51a) 

(3)/^o = . - sin(3^9) = 3W- sm^{l -^sin^^) . (D.51b) 

\ 71 \ TT 3 

The logical necessity of two angular states for iz = will readily become evident when 
we now look at the action of R± on the neighbouring states £^ = ±1. If we ascend the R- 
ladder left-hand side of Fig. D.I) beginning on the lowest-possible step (i. e. iz = —2) 
by repeatedly applying the ascendive operator R^ flD.22aP we first come to the step with 
£z = —1. The corresponding angular function is 

I 'iO 

= -\hr- e-^-^sin^T^cost? , (D.52) 
V oil 

and the ascent to this step reads (see Fig. D.I) 

^+^%-2 = -y|-(%-i. (D.53) 
Now we once more apply the ascendive operator and thus obtain 

i?+^%-i = -^^%o. (D.54) 

with ^^■'/ij,o- being given by equation flD.51bp . However, in order to further ascend now 
from the step with = to the higher steps £2 = 1, 2, 3 we cannot apply the ascendive 
operator R+ to that state '•^^/i?,o_ (1D.51b|) just reached; but we have to start from the 
other state ^'^•'/ij,o+ (1D.51aP which is also due to iz = 0: 

i?+(%o+ = -^v^-(%i. (D.55) 
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Indeed one is easily convinced by straightforward calculation that the action of -R+ on the 
state ^^^fRfl_ ( ID.Slbp does not lead us to the desired step ^^^/r^i given by 



= y ^ e^-^sin^? (^1 -^sinH^ . (D.56) 



On the other hand, once the right starting state '■^^/i?,o+ is accepted for the further ascent, 
the highest state ^^^fn^s is reached by repeated action of R+ without any inconsistency, 
see Fig. D.I. 

A similar effect is observed for the descent from the highest-order state ^^^fn^s to the 
lower states by means of the descendive operator i?_ (ID.33ap . see Fig. D.I above. Here, 



the first transition occurs from the highest state ^^^/r^^ flD.45p to the next lower one '•^^/ij,2 
flD.46p and looks as follows 

R-^'^fR,3 = V6-^'^fR,2 (D.57) 
and this yields for the transition coefficient C(3, 3) (ID.36aP 



C_(3,3) = V6, (D.58) 

in agreement with the general relation flD.38p . Then we can further descend downwards 
by means of -R_ to the step ^^^/r^i flD.56p . But the next transition downwards to the step 
iz = requires again some attention. By straightforward calculation one finds 

^_(%i = 2v^-(3)/^,o^ . (D.59) 

This says that the descent from above {i^ > 0) by means of R- terminates at that state 
^^^fR,o+ (ID-Slap on the level iz = from which we have to start when we wish to generate 
upwards the states with iz > 0, see equation fID.SSp . 



But for the further descent, one first has to pass over from ''^•*/j?,o+ flD-SlaP to the other 



angular state ^^^fRfl_ (ID.Slbp and then one can continue to descend, i. e. 

i?-^%o_=2y3-(3)/^,_i. (D.60) 
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But the next descent leads us to the last non- trivial step ^^^fR,-2 dEH 

i2_(%_i = v^(%_2, (D.61) 



since this last step becomes then annihilated by the action of see equation (ID.SOp . 



Concerning now the "S'-ladder", the situation is quite analogous to the "i?-ladder", 
cf. Fig. D.I and Fig. D. II. Therefore it may be sufficient to mention here only the most 
striking features. The first of these refers again to the fact that for £^ = one has two 
angular states, i. e. 

^%o+ = sm{3^) = -S^^- sin^ ~ t ^^^^^) (D.62a) 

^%o_ = cos(3^9) = ^cos^ (1 -Asin^^) . (D.62b) 

From here, the next higher state ^^^fs,i{'&, 0) is obtained by applying the ascendive operator 
S+ flD.22bj) to flD.62a|l 



^+^%o+ = -2v^-(%i (D.63) 
(^D+(3,0+) = -2y3) , 

whereas for obtaining the next lower state ^^^fs-i{'&,(f)) one has to depart from ^^^fs,o_ 
flD.62bp : 



^_(%o_ = -^v^-(%-i (D.64) 



Obviously, this is all quite analogous to the case with the i?-ladder; and the annihila- 
tion of the highest {iz = +3) and lowest {i^ = —3) states by means of the actions of 5*+ 
and S-, resp., works also quite similarly, cf. (1D.24P and (lD.34b|) . Clearly, the ladder co- 



efficient -D+(3, —3) is nothing else than the specialization of the general case Z)+(£p, —iv) 
(ID.39P to the present example £-p = 3 

D+{ip,-£p) ^ D+{3,-3) = -V6 , (D.65) 
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with the converse being given by 



D^iir, l-iv)^ D^{3, -2) = --^ , (D.66) 



cf. equation (lD.37bp . Somewhat more interesting are perhaps those relations concerning 



simultaneously both ladders, such as f lD.32p and flD.44p . From the i?-ladder (Fig. D.I) 

we see 

C+(3,0+) = -^v^, (D.67) 

and for the S'-ladder (Fig. D.II) 

L>+(3,0+) = -2^3 , (D.68) 

thus 

D+(3,0+) = ^C+(3,0+) (D.69) 

which validates (lD.32p . Similarly, from the i?-ladder (Fig. D.I) we take 

C_(3,0_) = 2^3 , (D.70) 

and from the S'-ladder (Fig. D.II) one takes 

D43,0_) = -^v^ (D.71) 



thus 



D_(3,0_) = -^C_(3,0_) (D.72) 



which validates the general claim (lD.44p . 



Thus, it should have become obvious that there are many interrelationships between 
the ladder coeffcients; but in place of studying this for special examples it would surely 
be more efficient to deduce their general shape (i. e. for general i-p and i^)- 
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S-ladder 





A 

s+ 



l,=+3 (3)/< 



D+(3,2)=0 



Z?+(3,l)=-V10 



S+ 



D+(3,0+)=-n/12 



D+(3,-l)=-V12 



D+(3-2)=-VT0 



S,-2 



r>+(3 -3)=-V6 



s+ 



S,-3 



S,3 



i3_(3,2)=-iVlO 



=+1 



^^^/s,o± = '-=o = 



5. 



^^-(3,1)=-^ 



S,0± 

D_ (3,0-)=- 1^3 



D_(3,-l)=-|v^ 



(^Vs,-2 



^-(3.-2)=-^ 



iz=-3 



(^Vs,-3 





Fig. D. II Ladder Operations for ^^*^Vs,iz(''^? 0)? = 3 
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Determination of the Ladder Coefficients 



For the determination of the general ladder coefficients C^li-pjiz) and D±{l'p^lz) it is 
very helpful to ffist determine certain products of them. Namely, by applying in subse- 
quent order an ascendive and/or descendive ladder operator one must come back to the 
original state of departure. Thus, one obtains the following eigenvalue equations for the 
operator products, e.g. for the i?-ladder 



R. ■ R+j = + 1) ■ C+(£p,4)J ■ 

and analogously for the S-ladder 



(D.73a) 
(D.73b) 



(D.74a) 
(D.74b) 



The eigenvalues emerging here for the special case i-p = 3 can be read off directly from 
both Fig.s D.I and D.II and may be collected in the following table: 

£-p = 3 





4^ 


-3 


-2 


-1 





+1 


+2 


+3 


R+-R^^ 


C+(3,4-l)-C_(3,4)^ 






-5 


-8 


-9 


-8 


-5 


R.-R+^ 


C_(3,4 + l)-C+(3,4)^ 




-5 


-8 


-9 


-8 


-5 






D+(3,4-l)-I^-(3,4)^ 




5 


8 


9 


8 


5 






D_(3,4 + l)-/^+(3,4) 


5 


8 


9 


8 


5 







With the coefficients of both types of ladders being known, one can now test also those 
"inter- ladder" relations such as flD.32p and flD.44p . For the ffist of these we find for the 
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present case £-p = 3: 



and therefore 



D+(3,0+) = -2^3 (D.75a) 
C+(3,0+) = -^v^ (D.75b) 



D+(3,0+) = ^-C+(3,0+) (D.76) 



which vahdates the general claim f lD.32p . The other assertion flD.44p can be validated by 
an analogous conclusion. 

Of course, it is not very satisfying to determine all the ladder coefficients C±{i'p,iz) 
and D±{l'p,lz) by explicit calculation; this method can only be of heuristic character 
in order to get some feeling of which numbers are to be expected. For the purpose of 
deducing the general result in form of a closed formula for all admissible ^z and one 
first writes down (by explicit calculation) the product operators, i. e. for the i?-ladder 

92 



i?+ ■ = 1 - (1 - cot^ d)-Lz- cot^ ^9.^2 + -— (D.77a) 
R--R+ = ■ 4 - coe^- Ll + ^ (D.77b) 



and analogously for the S'-ladder 

5+ • 5_ = -1- ■ Lz + cot^ ^-Ll-^ (D.78a) 
sm V ov'^ 

^_ ■ 5+ = -1 - (1 - cot^ d)-Lz + cot^ ^.Ll- — . (D.78b) 

Now observe that these operator products are required to act upon the angular functions 

^^^^Jr/z and *-^^^/5/^ where then the angular-momentum operator (ID.lOcp becomes 
replaced by its eigenvalue iz, i-e. we get for the i?- ladder 

^+ ■ i?_ ^ 1 - (1 - cot^ 7?) ■ 4 - cot^ (D.79a) 

R--R+^ ■ 4 - coe^- + ^ (D.79b) 
sm tf ov'^ 
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and analogously for the S'-ladder 

S+-S_^ ■ 4 + cot^^? . il-^ (D.SOa) 

sm 01}^ 

■ S+ ^ -1 - (1 - cot^ ^9) ■ 4 + cot^ i9 . 4 _ (D.80b) 

Furthermore, since the angular eigenfunctions ^^^''fR/^{'&, (j)) and ^^^''fs/^i'^,<P) must obey 
the second-order eigenvalue equations (IIV. ITap - QIV. 17bp . or (]D.19ap - (]D.19bp . resp., one 
can substitute from there the second-order derivatives (with respect to into the present 
operator products flD.79al) - flD.80bl which then finally yields their eigenvalues (ID.73ap - 
flD.74bp in terms of 4 and iv, i. e. for the i?-system flD.73ap and flD.73bD 

C+{ir, 4) ■ C_ i, + l)=il- 4 (D.81) 

and analogously for the S"- system (lD.74ap -( 1D.74bp 



D^iv, 4 - 1) ■ D^iv, 4) =rv-K- (D-82) 

From this result it is evident that the action of the operator products in reverse order 
entails merely a change of iz by on unity {iz =^ ±1), see the table for £p = 3 on p. (?). 
It should also be self-evident that the two independent equations ( ID.Sip and ( 1D.82P are 



not sufficient in order to fix completely the four ladder coefficients C±{£p,£z), i5±(£p,4) 
as functions of the two quantum numbers i-p and iz- At the most, we may conclude from 
these equations that the desired coefficients must be of the general form 

C^ivJ,) = X(lrJ, - 1) • (4 - <f - 1) (D.83b) 
4) = (Cp + 4 + 1) ■ 4 + 1) (D.83c) 

with unknown functions X(£p,4) and y(£p,4)- Surely, by comparison with Fig.s D.I 
and D.II one could perhaps guess how these unknown functions X(£p,4) and F(£p,4) 
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should look like; e. g. the guess 



would not be in conflict with the coefficients C±{ip,iz) of Fig. D.I; but of course it is 
better to rigorously work out the general form of the unknown functions X and Y. 

For this purpose, it is very helpful to first analyze a little bit closer the action of the 
ladder operators R± and S± along those arguments, which led us to the former special 
results flD.32p and flD.44p . Indeed, when one lets act these operators on that product form 
(lD.6a|) - (]D.6bl) of the angular eigenfunctions and furthermore substitutes the derivatives 
with respect to the variable from the first-order equations (lIV.8ap - (]IV.8bp one can 



represent the angular functions due to by those due to ± 1, i.e. for the i?- ladder 

^'^^fR,eM 0) = n -^'"^/^A-il^) (D-85a) 



■\iv-2\ <4 < 



(-|^P-1|<4<^P-1), 
and similarly for the S-ladder: 

^'^^fs,iA^, 0) = {(1 - 24) cot^ ■ + (D.86a) 

(-|^P-1|<4<^P-1) 



- (4 + 1) 
D_(£p,4 + 1) 



(D.86b) 



2). 
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Obviously, these inter-ladder relations do connect the angular eigenf unctions of either 
ladder to those of the other ladder, albeit for a neighbouring step. Therefore, if we wish 
to generate the totality of eigenfunctions for given we need not repeatedly solve the 
eigenvalue system (lIV.8ap - (]IV.8bp for any allowed value of nor apply the ascendive and 
descendive ladder operators but we can obtain the remaining unknown functions from the 
linear combinations of a known minimal subset of eigenfunctions. The coefficients of these 
linear combinations f lD.85aP - flD.86b|) are essentially the ladder coefficients C-|-(£p,£^) and 



D±{iv,^z) which thus yields additional motivation for determining their general form. 
Moreover, the present results (]D.85a|) - (lD.86b|) do also clarify the fact why the same 



angular functions (apart from the different pre-factors and steps) do emerge on both the 
i?-ladder and the S'-ladder (see the table for £-p = 3 on p. \5^ . 

The final task is now to check whether those inter-ladder relations f lD.85a|l -f lD.86b|) 
are compatible with the original eigenvalue equations (IIV.8a|) - (IIV.8bj) . This compatibility 
test is positive and yields one further inter-ladder relation, namely 

X{iv,iz)-Y{iv,iz + l) = l, (D.87) 

whose consistency may be exemplified by means of the ladder coefficients being displayed 
by Fig.s D.I and D.II for £-p = 3. But obviously this result is not sufficient in order to 
completely fix the ladder coefficients! The reason is here that we did not yet use the nor- 
malization condition flIV.15a|) for the angular eigenfunctions. But if this is implemented 



for one step of one ladder, then one can deduce all the other coefficients of both ladders. 
Thus, one obtains for the ascent of the -R-ladder 



C+(£p,4) = -(^P + 4)-^^^^^-^, -{^v-l)<iz<iv-l (D.88) 
and similarly for the descent 



C_(£p,4) = (£p + 1 - 4) ■ \ L t ^ ' (-1^^ - 2| < 4 < iv) (D.89a) 

V — + i 

C_(£p,l-£p) = 0. (D.89b) 
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For the S'-ladder, the ascent is described by the following ladder coefficients 

D+{^vAz) = -v/(^p-4)(^P + 4 + l) , -^r < 4 < (^P - 2) (D.90a) 

- 1) = (D.90b) 

and the descent by 



D^{ivA.) = -{^v-Q\j , -i^v-l)<i.<{iv~l) . (D.91) 

The ladder coefficients displayed by Fig.s D.I and D.II do again exemplify these results 
for i-p = 3. Zhe solutions f lD.88p -f lD.9ip do satisfy both product requirements f lD.8ip and 
f lD.82p and firthermore do also obey the following symmetries for the allowed values of Iz 
and Ip-. 

C+(/p,-/,) = D_(/p,y (D.92a) 
= -D+(/p,g (D.92b) 

Compact Representation 

Up to now we represented the angular-momentum algebra in component form by 
considering separately each component ^^'^^fR,iX'^-i4>) ^'^^ ^^''^^fR/A^t'P) of complex 
two vector (^^¥^^(^9,0) flDJ[3|l . But for a concise survey of the essential algebraic features 
it may be more advantageous to resort to a more compact abstract representation. The 
proper eigenvalue problem has already been reformulated in such an abstract form, cf. 
equations flD.14p - flD.18p . so that this must now merely be completed by the corresponding 
abstract representation of the ladder formalism. 

To this end, we condense the ascendive ladder operators /2+ and 5*+ f lD.22ap -f lD.22bp 
to the two-dimensional abstract object 7+ 



%=\ . (D-93) 
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and similarly for the descendive case 

( R \ 

r = \ 1 . (D.94) 

V 

The ascendive and descendive properties of these newly introduced objects are now ex- 
pressible by their action on the two- vector i.e. 

r+ (^^)f,, = A/V(£p, 4) ■ ('^^£.+1 (D.95a) 
fJ'^\ =M-{lvJ.)-^'^\-i , (D.95b) 
where the matrix-valued normalization constants M± are given by 





+ 4 + 1 



?P + 




?p + 



£p + + 1 

'p - 1 



+ 1 



(D.96a) 



(D.96b) 



£p - 

The equations (lD.95ap - (lD.95bp say that the action of 7+ lets increase the quantum num- 
ber iz of the eigenvector ^^''■*ff^ by one unit, and 71 lets it decrease by one unit. Indeed it 
ts easy to see that this behaviour is an immediate consequence of the following relations 



-T- 



(D.97a) 
(D.97b) 



which themselves are easily deducible from the corresponding commutation relations for 
R± and S±. 



L^, R± 
Lz, S± 



±R± 
±S± . 



Thus, we arrive at 



Czf- ^'^K = (4 - 1) ■ r_ ^'^^i^ 



(D.98a) 
(D.98b) 

(D.99a) 
(D.99b) 
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which express the actions of 7± without any reference to the normahzation constants Af±. 

With respect to the other quantum number £p angular momentum it is also instructive 
to consider the commutator of these abstract ladder operators T±, i.e. 



1-2Q, . 



(D.lOO) 



This result for the ladder operators may be exemplified by their immediate actions 
flD.95al) - flD.95b[) on the eigenvectors ^^'^^e^ where one can make use of the commutative 
relations of the (matrix- valued) normalization constants J\f±{ip,iz)'- 



U- 4) ■ M+{(.v. 4 - 1) - M+{lv. 4) • A4 4 + i) = i - 24a, 



(D.lOl) 



But of course, the claim fID.lOOp does hold quite generally and independently of the actions 
f lD.95ap - flD.95bp on the eigenstates of angular momentum. This can easily be verified by 
reference to the component form 



R-\-i R- 



1-2Q, 
1 + 2Q, . 



(D.102a) 
(D.102b) 



The result f ID.lOOp says that the sequence of ascent ('r+) and descent (f_) (or vice versa) 
does one not lead back to the original angular state ^^'^'^i^ (not even up to an overall 
pre-factor) . 
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